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Abstract 



With the introduction of computed tomography (CT) 
scanning in 19^2, advancements in diagnostic radiology and 
related disciplines have soared. This rapid growth has led 
not only to more advanced diagnostic technologies, but also 
to more advanced CT applications. A reference - display 
system has been developed to superimpose reconstructed CT 
images on the operative field, providing the surgeon with 
information that permits greater precision during operative 
procedures than presently possible. 

The CT images, reconstructed from conventional 
transverse oriented CT scans to match the surgical 
perspective, are displayed on a miniature CRT and 
superimposed at the focal plane of the operating microscope 
by a beam splitting assembly. The spatial registration of 
the CT data, operating . microscope and patient are 
accomplished by the adaptation of an ultrasonic digitizer, 
incorporating an array of three small spark gaps mounted on 
the microscope and used in conjunction with an array of 
ceiling mounted microphones. Glass fiducials, temporarily 
attached -to the patient allow the CT scans, patient and 
microscope to be optically referenced. 

The total reference system error was found to be 
approximately one millimeter. 
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Chaoter 1 



With the introduction of computed tocography (CTJ 
scanning in 19(2, advancements in diagnostic radiology and 
related disciplines have soared C 20 J . This rapid growth has 
led not only to more advanced diagnostic technologies, but 
also to more advanced CT applications. Through the use of 
computers, high resolution three-dimensional visualization 
of previously unseen structures is now available [20]. 

Integration of this technology with operative 
neurosurgery remains in its infancy. Stereotactic surgery 
utilizes a fixed reference frame attached to the patient's 
head to place E£3 electrodes, take tissue biopsies, etc., 
with better than 1 millimeter precision. This technique has 
the capability to incorporate the available computer 
technology with the operating room environment and many 
neurosurgical procedures. Concurrent with these 
developments, the field of neurosurgery in general has 
evolved to a higher art through the development of the 
operating microscope. We propose to integrate CT and 
computer technology with stereotactic principles and the 
operating microscope to develop a computer-based optical 
system for use during raicroneurosur-gical procedures. 

At present, a neurosurgeon's ability to perform 
intracranial procedures for tumor, vascular disease or 
functional disorder is dependent upon his mental integration 
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of the visualized operative field with his knowledge of 
neuro-anatomy and the available radiologic studies such as 
CT scans. Available technology could greatly improve that 
mental process and achieve a far superior degree of 
operative precision and safety. 

Conventional CT scans are oriented transversely to the 
body axis and as the operative approach is rarely along the 
axis of conventional scanning, the ability to reconstruct a__ 
CT scan to match the surgical perspective is highly 
appealing. ine major oojective of this project is zo 
develop a system that will superimpose reconstructed CT 
images over the field of view of the operating microscope. 
The neurosurgeon would then see, for example, the outline of 
a tumor superimposed on the operative field. 

There are a number of advantages of this 
reconstruction/projection system: 

1) There would be no dependence on the surgeon's mental 
reorientation of CT scan information. 

2) The information would be displayed such that it 
would not interfere with the neurosurgeon's procedure or 
require the reading of x-rays off a light screen some 
distance away. 

j) A computer-based anatomical atlas could be developed 
• that would superimpose on the operative field important but 
otherwise unseen structures, such as normal neuronal 
pathways and nuclei and major vascular structures [2J. 



4) The neurosurgeon could use the superimposed image as 
a map accurately guiding operative procedures with greater 
precision than presently possible. 

1 . 1 Recent Research 

Recent work in CT scanning and reconstruction has been 
well established and involves primarily image processing 
software and computer graphics [1,3]. This project will 
require original application of presently developed 
techniques . 

Adaptation of stereotactic technique to CT technology 
has been approached in a number of ways. One useful 
technique is that developed by Leksell at the Karolinska 
Institute in Stockholm, utilizing an adapted metal frame 
fixed to the patient's head at the time of scanning [26,27], 
(e.g. Figure 1.1). Stereotactic coordinates, relating the 
target position of CT demonstrated pathology to the 
stereotactic instrument, are generated directly from the 
scans and the patient is then transported to the operating 
room [27,30]. Other techniques are adequate but often more 
cumbersome [4-7,9,12-17,21,22,28,29,31-37,39-42,48]. All of 
these enable stereotactic procedures generally characterized 
by "blind" insertion of needle-like instruments through 
small openings utilizing previously obtained CT-determined 
landmarks.. This has been a vital development and a 
foundation for this project. They have not generally been 
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Figure 1.1 - Stereotactic Frar.ie 

amenable to "open" procedures such as craniotomy for tu:nor 
or vascular disease and, as previously noted, do not allow 
access to CT data after selection of a target. The CT 
information utilized is limited to the coordinates of a 
point. All instruments select a target and set the 
instrument accordingly; the proposed system, operating in 
reverse, allows free positioning of the microscope with 
subsequent stereotactic "positioning" of the CT data. 

The operating microscope has been incorporated into CT 
stereotactic work by Kelly at the State University of New 
York at Buffalo [2j,24j. This development has also employed 
surgical laser instrumentation and shown the feasibility of 



achieving a synthesis of technologies and the prospects of 
refining neurosurgical operative technique. His technique 
of linking the operating microscope and the stereotactic 
reference system requires utilization of a large 
stereotactic frame. The proposed system will eliminate the 
encumbrance of such a frame and in doing so permit a 
potentially wide-spread applicability to general 
neurosurgery. Ally's system has not employed a projection 
system as proposed in this project. 

1.2 Problem Statement 

The' goals of this project are to design a system that 
will reconstruct CT images from diagnostic CT scans and 
superimpose the reconstructed image on the magnified field 
of the operating microscope. 

The overall project can be divided into four parts. 

1) The data from the CT scanner/computer must be 
transferred from the computer's memory into another form 
that can be manipulated in the operating room. 

2) Computer algorithms and programs must be developed 
to calculate the new reconstructed CT image based on the 
position of the operating microscope with respect to the 
patient's neuro- anatomy , as presented in the CT scans. 

i) A system must be developed to display the new 
reconstructed CT image within the operating microscope. 



U) Ihe operating microscope must oe linked or 
referenced to the coordinate system of the CT scan and the 
patient. 

Our initial design will be to reconstruct only certain 
CT boundaries or contours (su.ch as the outline of a brain 
tumor) and if this technique proves useful, future 
developments will reconstruct grey scale CT information/ 

1 .3 Thesis Objectives 

The specific objectives of this thesis are to address 
the last two points of the Problem Statement and to 
construct the hardware necessary for displaying referenced 
CT scans within the operating microscope. In order to more 
clearly define the reference - display system, the 
reconstruction technique must be outlined. It has been 
determined that manipulation and reconstruction of CT images 
would be best handled by the Treatment Planning Computer (a 
Data General Eclipse SljO based system manufactured by 
General Electric) located in the Norris Cotton Cancer Center 
at the Mary Hitchcock Hospital, since this system has the 
capability to read the formatted CT scans stored on magnetic 
tape by the G.E. CT Scanner. Figure 1.2 shows a block 
diagram representing the overall system. 

The objectives of this thesis are: 
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Fiourg 1.2 - Sustem Block Diagram 



1) to provide the Treatment Planning Computer with 
enough information . concerning the position of the operating 
microscope with respect to the patient and CT scans to 
reconstruct the appropriate image; 

2) to superimpose -the reconstructed image onto the 
field of view of the operating microscope; 

3) to desi'gn and build all the hardware necessary for 
the reference - display system and 

4) to evaluate the system's performance, 

1 ,3.1 Specifications 

There is only one specification concerning the design 
of the reference - display system: the system should 
display a reconstructed CT image of the correct structural 
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anatomy within 1 millimeter (comparable to the stereotactic 
frame) or prove that constraints prevent this accuracy from 
being achieved. 

I.j.2 General Constraints 

This system must meet the constraints of the operating 
room environment and not impair or reduce the effectiveness 
of the operative procedure. The system must also be simple 
to use, allow modifications for new applications, and remain 
within budget restraints of "J 1 *, 500 . Specific constraints 
are indicated in each chapter. 

1.4 Thesis Organization 

This Thesis is divided into twelve chapters and seven 
appendices. In general, each chapter describes a different 
aspect of the project and includes specifications, 
constraints, alternatives and a "description of the selected 
alternative. Chapter 2 describes the operating room 
environment and the constraints associated with 
neurosurgical procedures. Chapter 3 describes the optical 
display system and, more specifically, the choice of beam 
splitter and image display. 

The remaining nine chapters address the problem of 
referencing the reconstructed CT image to the position of 
the operating microscope. Chapter 4 outlines proposed 
alternatives for the reference system. The selected 



registration technique, an ultrasonic linkage, is discussed 
in Chapter 5. Chapter 5 explains the various coordinate 
systems and Chapter 7 describes the theoretical error 
analysis. Chapters 3 through 10 detail the specific 
hardware design of the spark gap multiplexer, reference 
system software, and mechanical designs of the spark gap 
holder and microphone mounts, respectively. The system 
evaluation and experimental results "are described in Chap^el 
11, and Chapter 12 concludes the Thesis with an outline of 
future work and a future design. 

The appendices include schematics, test statistics, 
software, mechanical designs, experimental results and a 
digitizer operator's manual. 



Chapter 2 - Operating Hooai Inv ironment 

The operating room environment imposes three general 
constraints on any procedure: 

1) physical constraints on the positioning and movement 
of equipment; 

2) aseptic constraints on the cleanliness of equipment 
in the operating room and within the operative field; 

3) safety constraints on electrical equipment used 
within 15 feet of the patient. 

2.1 Physical 3^ri3tr3int5 

For the purpose of this project, the operating room can 
be broken down into six basic components that place 
restrictions on the positioning and movement of equipment: 
the patient, the surgeons, the scrub nurse and surgical 
instruments, the anesthesiologist and anesthetic equipment, 
the operating microscope," and the ceiling-mounted operating 
room lights. £ach of these components are in different 
positions within the operating room depending on the 
surgical procedure. There are three neurosurgical operative 
setups depending on the type of surgical procedure and the 
craniotomy site. Procedures involving operations with the 
patient in the sitting or supine positions for frontal, or 
right temporal craniotomies are set up as shown in Figure 
2.1 and 2.2. Procedures involving left temporal 

10 
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Supine Le-Ft Temporal Cramotomu 

craniotomies are set up as shown in Figure 2.3. These 
setups allow sterile draping such that the anesthesiologist 
has proper access to the patient. The operating lights are 
attached to ceiling-mounted light tracks, between which is 
an air output vent. See Figure 2.*). 
2.2 Aseptic Constraints 

Every instrument, person, and piece of equipment in the 
operating room must be either clean or sterile. All items 
must be clean - i.e. free of dirt, oils and other 
contaminants conducive to the growth of potentially 
infective organisms. Equipment brought into the operating 
room must withstand a general cleaning procedure. Items 
directly and indirectly involved with the operative field 
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must be sterile, i.e. not only clean, but ideally free of 
all microbial organisms. This includes anything touching 
the operative site. 

Wot all equipment within the sterile field must itself 
be sterile, but it must be covered with a sterile barrier. 
Examples of items that must be sterile or contained within 
sterile drapes include surgical instruments, the surgeon's 
hands, arms and front body and anything hanging over the 
operative field that has the potential to come in contact 
with sterile areas. 

There are two different sterilization procedures 
depending on the material to be sterilized. Stainless steel 
and other metals that will not oxidize if exposed to steam 
or aqueous solutions must withstand a temperature of 243°F 
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for 20 minutes. Teflon and plastics must undergo ethylene 
oxide 2as sterilization at 1j0°f for 24 hours [44]. 

2.3 Safety Constraints 

Any electrical equipment used within 15 feet of the 
patient must meet the National Fire Protection Association 
(NFPA) safety standards [44]. These standards include: 
proper three conductor a-c power line cord with 
nospi cai/oper acing room grade plug, properly connected 
grounding pin with no more than 100 microamps of leakage 
current, and a resistance between the grounding pin and 
chassis of less than 0.2 ohms. 

The most difficult of the standards to maintain is the 
leakage current limit. The leakage current can be reduced 
by powering the equipment through an isolation transformer. 

2.4 Operating Microscope 

The operating microscope used for neurosurgical 
procedures at the t-Iary Hitchcock Hospital is a Zeiss 0PMI-1H 
with a Contraves floor stand, (Figure 2.5). The microscope 
is used to magnify the operative field approximately five to 
thirty times. The optics incorporate either a 250 or 300 
millimeter objective lens with five parafocal lens settings, 
allowing five magnifications with little focal length 



Figure 2.5 - Zeiss Operating Microscope 



adjustment. The range of the magnified operative field is 

9.5 to 1.45 centimeters in diameter. See Table 2.1 below. 

Table 2.1 - Optical Specifications 

Magnification Magnification Q 1 9 J Field Depth of 

Setting Diameter Field 

in cm. in mm. 

2.5 13. Ox 2 

1.6 8.3x 2.30 4 
1.0 5.2x 3.70 7 
0.6 3.1x 6.10 19 
0.4 2.1x 9.50 >19 

The internal ocular- to-objective optics are all focused to 

infinity, allowing for the addition of other optical 
equipment, i.e. irises and beam splitters, without changing 
the optical characteristics of the microscope. The 
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aicroscope also incorporates a through- the-objective 
illumination system with two degrees' of brightness. 

The aicroscope mount is a six jointed, six degree of 
freedom Contraves floor stand with pistol grip, and a foot 
pedal or mouth switch activated servo-aechanisas to lock the 
microscope in place, ^hen the microscope is properly 
counter-balanced, the surgeon can use the mouth switch to 
release and reposition the "floating" aicroscope. 

The microscope is also equipped with a beam splitting 
device for the attachment of a stereo-observation tube =nd 
either a color video camera or a 35 millimeter camera. 

To insure proper aseptic technique, specially designed 
sterile transparent plastic bags are draped over the 
microscope and most of the floor stand, (Figure 2.6). 




Figure 2.6 - Draped Operating Microscope 



Chapter 3 - Optical Display System 

This chapter describes the constraints and alternatives 
of the optical display system and the video connections 
necessary for displaying the. reconstructed CT image. 

3.1 Constraints 

. To display a reconstructed CT image at the focal plane 
of the operating microscope, a beam splitting device is 
introduced into the optical path of the microscope. This 
device can be used to superimpose or "add" two images. See 
Figure. 3-1. The constraints on the beam splitting assembly 
ar e : 
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Figure 3. 1 - Beam SpMttar 
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1) The beau splitter must have a reflection - 
transmission coefficient ratio such that superimposed images 
are not "washed out" by the illuminated operative field. 

2) The prism optics must not introduce any visually 
disturbing glare or aberrations. 

3) fhe additional optical path length must be less than 
25 centimeters to maintain an acceptable surgeon eye-to-hand 
distance of less than 45 centimeters. 

4) The weight, of the beam splitter must be less than 4 
kilograms dnd the size must be less than 10x20x10 
centimeters to remain within limits for balancing the 
microscope and not interfere with the surgical procedure. 

The image source of the reconstructed CT scan must be 
introduced into the beam splitting assembly in order to be 
superimposed on the focal plane of the microscope. Since CT 
scans are computer generated images displayed on a video 
screen (unlike conventional X-rays that are produced on 
translucent, plastic X-ray film), this CRT (Cathode Ray 
Tube) image would provide an appropriate and convenient 
image source for the beam splitter. The constraints 
governing the image display are: 

1) The display screen must be bright enough to 
effectively superimpose an image on the illuminated 
operative field. 
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2) The CRT must be easily mounted on the microscope, 
weighing less than 1 kilogram for proper balancing, and not 
interfere with the surgical procedure. 

3) The CRT must possess enough resolution to display a 
crisp image. 

j,.2 Zeiss Drawing Tube 

" The alternatives for a beam splitting assembly are 
either to purchase a manufactured assembly, design our own, 
or combine the best features of both. Since Zeiss 
manufactures the microscope and assumes liability for its 
operation (any piece of non-Zeiss equipment attached to the 
microscope would void the Hospital's liability contract), 
manufactured assemblies are preferable and were investigated 
first. A beam splitting assembly manufactured by Zeiss is 
optically compatible with the OPMI-1H operating microscope. 
It is sold as a drawing tube to allow biological 
illustrators to superimpose their sketch pad on the 
magnified field for more accurate illustrations. The model 
that is optically compatible with the microscope is ifM7M622 . 
See Figure 3.2. It consists of a 50/50 binocular beam 
splitter (that will mount between the objective lens of the 
microscope and the observation tube/video catoera beam 
splitter) and a projection tube approximately 20 centimeters 
long and U centimeters in diameter. There are two movable 
lenses in the tube; one is for focusing the superimposed 



Firure 3.2 - Zeiss Drawing Tube 



image and the other is for magnif ication . A prism at the 
end of the tube directs the drawing surface through the 
projection tube to the beam splitting prism. Since the 
0PMI-1H microscope is an older model, the drawing tube 
presently available was mechanically modified to fit the 
microscope. This modification, done by Fred Schleipman, 
involved reducing the size of the beam splitter housing by 
1.27 centimeters on one side and replacing a thumbscrew with 
a recessed set screw. The drawing tube adds approximately 4 
centimeters to the optical path length and, for a proper 
mechanical fit, a 2.5 centimeter spacer must be added to 
lift the beam splitter above the focusing knob. 
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3.3 CRT Display 

There are -lany CRT's that will display composite video 
images from the Treatment Planning Computer, but few that 
will meet the constraints of size and weight. A screen size 
of 5 centimeters square or less is preferable. Most 
portable video cameras have electronic viewfinders 
consisting of miniature CRT displays which accept composite 
video inputs and would be compatible with the beam splitter. 
An exectronic viewi'inaer unodel ,/vi-190u; was purcnasea fro a 
J.V.C. that has a 38 millimeter diagonal black and white CRT 
and all the display driving electronics, accepts a standard 
composite video signal, weighs approximately 0.5 kilograms, 
is small enough to not interfere with any surgical 
procedure, has brightness and contrast control, and has very 
good image resolution. See Figure J. 3. 

Several simple modifications were required to make the 
viewfinder compatible with the Zeiss beam splitter. 3UC 
connectors were added for the video input signal and +12 
volt dc power supply. The field lens and flexible rubber 
eye gasket were removed on the viewfinder and the mirror 
prism assembly was removed from the drawing tube to allow 
the plastic housing of the viewfinder to close around the 
end of the projection tube. To ensure a rigid fit, a hard 
rubber coupling was contact cemented to the projection tube 
and a small metal alignment pin was machined to fit the lens 



Figure 3.3 - Miniature CPJ Display 



guide, preventing twisting and drooping of the display while 
attached to the microscope. The angle of the display and 
drawing tube was chosen not to interfere with the surgeon's 
hands. See Figure j.H. This type of attachment was made as 
an initial connection. Should the surgeon later- find this 
cumbersome, or a more permanent structure becomes necessary 
(i.e. such that the projection tube be shortened .or just the 
CRT itself be attached to the projection tube) this design 
can be easily modified. Figure 3-5 shows the attachment of 
the beam splitter - CRT assembly to the operating 
microscope . 



Figure 3.4 - Beam Splitter - CRT Assembly 




Figure 3.5 - ilicroscope - Display System 
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3.4 Interact N Connections 

The work station for the Treatment Planning Computer 
has a composite video signal output. The reconstructed CT 
image that can be displayed on the work station monitor must 
be displayed on the miniature CRT in the operating room. 
Therefore, a video connection is necessary. Since the work 
station is in the basement of the Hospital and the operating 
suite is on the second floor, a simple, quick video 
connection «ao unlikely. The alternatives for making this 
connection were to string a video cable from tne work 
station to the operating room or to use some existing 
connections. Since the Hospital has an audio-visual system, 
Interact T'/, that can provide this type of connection, the 
last alternative is viable and was selected. Radiation 
Therapy also has audio-visual lines connecting various rooms 
in the department, including the computer room. Upon 
consulting Ray Culig at Interact TV, a patch design was 
made. It involved connecting the work station output to the 
Radiation Therapy network, which terminates in the Hospital 
sub-basement, and connecting the sub-basement signal to the 
Interact T7 lines that run up to the fifth floor where a 
video patch connects the line to the Interact Control Room 
on the first floor, from there the signal is patched to the 
operating suite - operating room 10 (which conveniently is 
used for neurosurgery). See Figure 3.6. These video lines 
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Figure 3, 5 - Interact TU Uideo Connections 



were tested and there is no need for signal . amplification or 
frequency compensation (high end boost). .This dedicated 
video line has been rented until March, 1935. 

One other modification was required on the miniature 
CRT. The work station was manufactured in Europe and 
consequently has a 50 Hz video " display scan rate which is 
incompatible with the miniature CRT. The horizontal hold on 
the CRT did not have enough range to display a stable image, 
so the horizontal hold potentiometer and a resistor were 
replaced to allow enough range for a stable image. See 
Appendix A for modified schematics. 



This chapter describes constraints (4.1) and 
alternatives (4.2) for the design of a reference technique 
that will coordinate the registration of the reconstructed 
CT image with the physical anatomy as observed through the 
operating microscope. The specific requirement of such a 
reference system is that the information from the diagnostic 
Ci scans be reconstructed to match the focal plane from the 
operating microscope and superimposed in the microscope 
within 1 millimeter of the .correct anatomy. Therefore, 
enough information 'rjst be available to calculate the 
equation of the focal plane and its orientation with respect 
to the microscope. 

4 . 1 Constraints 

The constraints on the reference system are: 

1) It must meet the physical, aseptic and safety 
constraints of the operating room. 

2) It must allow clear visualization of the magnified 
operative field and, if needed, only require a minimum of 
alterations to the microscope. 
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4.2 Alternatives 

4.2.1 Stereotactic Linkage 

Since the stereotactic frame (described in Chapter 1) 
determines its own coordinate system, mechanically linking 
the microscope to the frame after the patient has been CT 
scanned would provide a direct technique for determining the 
position of the focal plane. Such a linkage is shown in 
Figures 4.1 and 4.2. unfortunately, tne scereocacuie frame 
not only limits the positioning of the microscope, but it 
also restricts the operative field and craniotomy site, 
making this system applicable to few procedures. This 
design was, therefore, not considered. 

4.2.2 Cranial Linkage 

To avoid the problem of the stereotactic linkage, a 
mechanical connection could be made from a small plate 
.rigidly attached to the patient's skull (away from the 
craniotomy site) to the microscope. This cranial linkage 
could include the minimum number of joints (monitored by 
transducers) to allow the microscope enough degrees of 
freedom for proper motion. A coordinate system could be 
established by CT scanning the patient with the skull plate 
in place, detecting at least three known points on the CT 
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Figure -4.1 
Stereotactic Linkage - Front View 
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Figure 4.2 
Stereotactic Linkage 
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Figure 4.3 
Cranial Linkage 



scans. The microscope's position could then be monitored 
with respect to this coordinate system. See Figure 4.3- 
U.2.3 Microscope Stand Redesign 

The microscope stand could be redesigned to translate 
and rotate with respect to a fixed coordinate system. The 
coordinates of the focal point and the equation of the focal 
plane could be determined directly from the scales on the 
microscope stand or by computer-monitored transducers. 




In order to relate the equation of the focal plane in a 
fixed coordinate system to the- CI scans, the relative 
location of the Cf coordinate system must be known. This 
can be achieved by determining the transformation matrix 
(see section 4.2.4.2) to convert coordinates in the fixed 
coordinate system to coordinates in the CT coordinate 
system. Transformation matrices can be defined by 
determining three linear translation constants and three 
rotational angles - a total of six pieces of information. 
Since angular information cannot be directly derived from 
the relative positions of the fixed and CT coordinate 
systems, the coordinates of at least three non-coll inear 
points (nine pieces of information), common to both 
coordinate systems, must be determined. The transformation 
matrix can then be calculated with six of the nine 
coordinates. However, using only two of the three points 
will not uniquely relate the positions of the coordinate 
systems about the line through the two -points. 

The common points can be three CT-detectable markers or 
fiducials placed on the patient's head before CT scanning. 
The fiducials will appear in the CT scans and can be 
assigned coordinates by the scanning computer based on the 
CT coordinate system. When the patient is brought into the 
operating room the microscope can be focused on each of the 
fiducials and their coordinates determined in the fixed 
coordinate system. Once the transformation matrix is 
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calculated the equation of the focal plane can be determined 

in CT coordinates. 

This alternative might be the best design once the CT 
scan - operating microscope technique has been proven 
useful. It would be too time consuming and expensive to 
consider at this point. 

4.2.4 Microscope Position Tracking 

By mounting position transducers on the six joints of 
the microscope stand (design I) or connecting the microscope 
with a multiple joint linkage to a fixed point with respect 
to the patient (design II), the microscope's position could 
be monitored with respect to a fixed coordinate system. 
Relating the fixed coordinate system of the microscope to 
the CT scans would be accomplished as described in section 
4.2.3 above. This design would allow almost all procedures 
to be carried out normally. There would be no mechanical 
connection between the patient and the microscope to 
restrict the surgeon's procedure, and the microscope could 
be repositioned anywhere once a fixed coordinate system has 
been defined. The angular sensitivity required of the 
transducers for a less than 1 millimeter accuracy at the 
focal plane might require expensive transducers and 
sophisticated driving electronics. Also, applying the 
transducers would involve mechanical modifications of the 
Microscope stand,, which was to be avoided if possible. 
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Since this design appeared practical, it was 
investigated further. 

; 

' H.2.H.1 Theoretical Linkage Analysi s 

In theory, at least , six pieces of information must be 
* known about a rigid body (microscope/focal plane) to locate 

it precisely and without redundancy in a fixed coordinate 
space, see Figure 4.H. Therefore, any position locating * 
technique must determine or monitor six degrees of freedom 
[13]. 

From the mobility formula Uutzbach criteria) for 
■ f spatial mechanisms (used commonly for closed kinematic 

chains, but whicn can also give degree of freedom 
information for open chains): 
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- Pi end Sodu fin al us is 



a = 6(n-1) - 5 j, - 4j 2 -jj 3 - 2 j 4 - j. 

wnere m = * of degrees of freedom, n = ,/ of links (rigid 
connections between joints) and j ] to j g = if of joints with 
1 to 5 degrees of freedom [42]. 

Six degrees of freedom are required, so m = 6 , and 
since the basic single-degree-of-f reedotn , revolute joint is 
the :nost common and easily monitored by transducers, j^ to 
j- = 0 and j 1 = 6. To determine the number of required 
links , ii : 

6 = 6(n-1 ) - 5(6) , n = 7 links. 
See Figure 4,5. 

4.2.4.2 Coordinate System Transformation 

Transformation matrices are used to transform 
coordinates from one coordinate system to another. The 
position of the focal plane can be transformed into a fixed 
coordinate system by knowing the appropriate translations 
and rotations. Rotations can be characterized by a 3xj 
matrix, indicating the axis of rotation, and a vector of the 
initial coordinates (x,y,z) [25]. See Figure 4.6. 

Note: All rotations are counter-clockwise about the 
rotational axis while looking toward the origin. 
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Multiplying the coordinates of a point in one 
coordinate system represented by a 1x3 column vector by a 
rotational matrix produces the coordinates of the point in a 
rotated coordinate system. Translations can be represented 
by a 4x4 matrix, where T x , T y and T z are the magnitudes of 
the vectors along the x, y, and z axes, as shown below. 

10 0 0 

T x 1 0 0 

T y 0 1 o 

0 0 1 

Combining the two matrices yields a 4x4 matrix with the 
general form: 



1 

T x 

Z 



0 


0 


0 


p2 1 


*12 

R 22 


R 1 3 
r 2 3 


R 31 


R 32 


R 33 


is the translation 


vector 


along 



z axes respectively, while R is the. rotational matrix. 
Written in this form, the translation takes place before the 
rotation [38]. Point(s) in one coordinate system can be 
transformed into another by multiplying the appropriate 
transformation matrices together. Six degrees of freedom 
requires six transformation matrices. Since matrix 
multiplication is not commutative, the multiplicative order 
of the transformation matrices is important. 




Angular Measurement Alternatives 

3oth designs I and II depend on the accurate 
measurement of rotational angles. These angles can be 
measured electronically (providing a more precise 
processor-usable input) by using position transducers that 
convert a mechanical displacement into an electrical signal. 
Position transducers can give two types of output: analog 
and digital. Since the position of tne microscope and the- 

latter output is preferred. 

Transducers presenting digital output for an angular 
displacement are called shaft encoders, and their output 
sensitivity (resolution) is measured in numbers of bits 
(angular fraction of j60°) . There are two different kinds 
of shaft encoders: incremental and absolute. Incremental 
encoders give a pulse output such that by counting the 
number of pulses, the relative angular displacement can De 
found with respect to a user-defined zero value. Absolute 
encoders have a fixed binary output for each position of the 
shaft. Absolute encoders are preferred for this application 
since the electrical signal can be directly sent to a 
processor without pulse-counting electronics. Table 4.1 
gives price ranges for various resolutions (ill of bits). 



Table 4.1 - Absolute Shaft Encoders 



Resolution 



Price Range 



9 bits 

10 bits 

11 bits 

12 bits 

13 bits 

14 bits 



$400. 

$100. 
$2300. 
$2300. 
$2400. 
$3200. 



$2500. 
$3000. 
$3000. 
$3000. 
$3200. 
$4000. 



The most accurate analog output transducers are rotary 
variable differential transformers (RVDT). The accuracy of 
RVUT's is measured in percent nonlinearity - how well the 
voltage output versus displacement curve conforms to a 
linear relationship. This percentage increases with an 
increase in angular range. Table 4.2 shows the prices for 
various resolutions (converting percent nonlinearity to 
angular error in degrees) for Pickering RVDT's. 



Table 4.2 - Pickering RVDT ' s 



RVDT Model J 



Range 



Resolution 



Price 



23501 



23511 



23300 



23380 



0- 
0- 
0- 
0- 
0 
0 
0 
0 
0 
0 
0 
0 
0 




+0.100° 
TO. 275° 
TO. 650° 
TO. 875° 
TO. 025° 
TO. 120° 
TO. 525° 
TO. 010° 
TO. 100° 
TO. 875° 
TO. 005° 
TO. 060° 
TO. 700° 



$143.35 



$210.45 



$247.95 



$27f.20 



J3 

4.2.4.4 Microscope Stand Analysis 

In order to estimate the accuracy and costs associated 
with tracking the position of the microscope, the microscope 
stand design was further analyzed, figure 4.7 represents 
the microscope stand. The links were measured as accurately 
as possible (+ 5 millimeters) and a transformation matrix 
was determined for each link-joint combination. The 
starting point vector provides the coordinates of the focal 
poms wi:n respect to an assignee Cartesian coordinate space 
whose origin is at joint 1. The angular range for each 
joint is indicated. See Figure 4.3. 

A computer program was written to evaluate the 
transformation matrices and determine the tolerable 
rotational error for each joint, given the maximum focal 
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point displacement error of 1 millimeter. Table 4.3 below 
shows the rotational error in degrees, the number of parts 
of a circle (fraction of 350°) and the number of bits of 
resolution required for each joint. 

fable 4.3 - Rotational i£rr or 



Joint Rotational Error $ Farts Bits 

1 +0.392° 920 10 

2 +0.031° 4445 13 

3 ±3.070° 5145 lj 

4 +0.054° 5670 13 

5 . +0.057° 5j20 13 

— ^ . , o . - - 



The minimum referencing cost for the microscope stand 
linkage using shaft encoders is approximately 313,200 and 
using RVDT's, approximately $1,500. Analog- to-digital 
conversion circuitry would have to be built for the R/DT's, 
so the cost involved would be greater, but less than that 
for shaft encoders. This evaluation was based on the 
assumption that the total system error is due to angular 
microscope error and that errors in other parts of the 
display - reference system were zero. The actual design 
must be much more accurate. 

4.2.5 Ultrasonic Linkage 

Another technique for locating the microscope/focal 
plane would be to use ultrasonic principles to determine the 
distances between the microscope and fixed ultrasonic 
transducers. This would require a transmitting - receiving 
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network to send, receive and interpret the ultrasound. 
Science Accessories Corporation (S.A.C) manufactures such a 
system, called a three dimensional ultrasonic digitizer. 
The system works by accurately measuring the transit time of 
a sound wave from a small spark gap source to a point 
microphone. Multiplying the transit time by the speed of 
sound in air, compensated for changes in temperature, 
produces the distance from the source to the receivers. Tne 
distances from the source to at least three receivers 
(referred to 25 slant ranges) with known orientation can 
then be triangulated to produce the coordinates of the 
source in the receiver coordinate system. Since we are 
trying to monitor the relative positions of two rigid 
bodies, a minimum of three sources and three receivers are 
required. 

More specifically, the system works by triggering a 
spark gap which emits an audible wide band frequency 
pressure wave ("click"). At the same time a start signal is 
generated which turns on high frequency counters associated 
with each microphone. The microphones "listen" for a signal 
between 55 and.Go' kilohertz and send a stop signal to the 
counters when the sound is received. A temperature 
sensitive transistor, mounted on a preamp box next to one of 
the microphones, monitors the ambient temperature and 
adjusts the number of clock cycles per centimeter 
accordingly. The only major constraint on this system is 
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that the spark gap source must have a clear line of sight to 
at least three microphones to determine its relative 
position in space. 

Science accessories offers two models of the three 
dimensional digitizer. Model 1 ($10,700) outputs only slant 
range distances which are accurate to +0 . 1 J of the measured 
slant range distance and are precise to within +0.01 
centimeters. Model 2 ($12,000) converts the slant range- 
values into C3rtesian coordinates and requires the 
microphones to be accurately positioned in a specific array. 
3oth models come with four point microphones to allow system 
redundancy or to make sure at least three microphones meet 
the line of sight constraint. 

4.j Conclusions 

Given the above alternatives for a reference system, 
microscope position tracking (Section 4.2.4) and the 
ultrasonic linkage (Section 4.2.5), and based on an estimate 
of performance and constraints it was determined that the 
ultrasonic linkage appeared to be the best choice because it 
meets all the physical constraints on the reference system 
and will not interfere with' the surgical procedure, it will 
cost less than most of the other designs, and it 
incorporates a more elegant design. A tfodel 1 (slant range 
output) ultrasonic digitizer was purchased from S.A.C. with 
three spark gaps and an RS-232 serial communications board 
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to transmit the data to a host computer, an ISM PC XT. The 
important specifications are listed below. 

1) The spark gap - microphone distance (slant range) is 
accurate to within +0.1* of the measured distance. 

2) < The resolution or precision of the slant ranges are 
within +0.01 centimeters. 

3) Trie maximum slant range is approximately 250 
centimeters. 

4) The maximum sampling rate (transmitting data at 9500 
baud) is approximately jO points per second. 

The accuracy of the digitizer was investigated further 
"to determine its limitations. 



This chapter describes the fundamental constraints on 
the digitizer based on the physics cf sound propagation in 
air (5.1), the sensitivity of the digitizer to temperature 
variations (5.2) and counter error (5.3), and an ultrasonic 
linkage description (5.4). 

5.1 Physics of Sound Propagation in Air 

The purpose of this section is to .determine the 
sensitivity of" the speed of sound in air to various 
parameters such as temperature and relative humidity. 

It is assumed that air behaves as an ideal gas for 
normal atmospheric pressures (valid up to several million 
.iewtons per square meter) and that for sound propagation its 
compression as a pressure wave is isentropic (reversible and 
adiabatic). Applying the laws of thermodynamics to sound 
propagation in the operating room environment' with relative 
humidity 20 to 60S and temperature of 19 to 24°C, gives the 
results in Table 5.1 below [11 J. 



Table 5.1 


- Speed of 


Sound in Air 


T (°C) 


* 


c ( m/ s ) 


19 


20'JS 


342.788 


50* 


J 4 3. 257 


24 


20* 


346.150 


60* 


346.616 
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This Table indicates that the speed of sound in air, c, 
is much more sensitive to variations in temperature tnan 
variations in relative humidity. 

5.2 Sensitivity to Temperature 

The sonic digitizer is based on the simple principle, d 
= ct, where d is the distance calculated in meters, c is the 
speed of sound in meters/second, and t is the transit time ^ 
between firing the spark gap and sensing the signal at the 
microphone. Since the speed of sound in air is aost 
directly a function of temperature and the temperature is 
only monitored at one of the microphone preamps, there must 
be some fundamental relationship between the variations in 
temperature over the path and the accuracy of the digitizer. 
Letting A = kR/mw and c = AVT (from Section 5.1) and 
applying the chain rule of differential calculus to d = ct, 
where k is the ratio of specifid heats of air at constant 
temperature to that at constant volume, R = 831H.3 
J/kg-mole-°rC is the universal gas constant, raw is the 
molecular weight in grams/mole , and T is the absolute 

o,, 

temperature in 

d = AtVf, 
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Since we are interested in the sensitivity of d with respect 
to changes in temperature within the active volume only, let 
At = 0 . 

Ad = f^AT = AtAT/2VT 
relative % change in d = Ad/d 
Ad/d = AtAT/2dVT = AtAT/2AtVTvT 
Ad/d = AT/2T 

Since the accuracy of the system is + 0.1* of the measured 
distance for slant range calculations, 

0.001 = AT/2T, 
AT = 0.002T(°.<) , 
AT. n . n = 0.002(283. 15°0 = 0.5Y°C, 
AT v = 0. 002 (297. 59°«0 = 0.60°C, 
AT = 0.6°C. 

Therefore, AT is an estimate of the a:nount the temperature 
can vary over the sonic path and deviate from the 
temperature determined at the microphones. This indicates 
that the percent change in the slant range values 
(specification 1 from Section 4.3) can tolerate a 
temperature variation of up to 0.5°C before exceeding +0.11. 
This assumes that temperature is the only factor affecting 
the speed of sound in air, which was concluded in Section 
5.1. 
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5.J Sensitivity to Counter Error 

The basic counter clock cycle frequency, 3.6406 MHz, is 
adjusted to 100 counts per centimeter or 0.1 ram per count, 
based on the temperature reading at the microphones. 
Therefore, since the counters can be in error by up to one 
clock cycle due to the counter resolution, the counter error 
is + 0.01 centimeters. 

5.4 General Ultrasonic Linkage Design 

The reference system must determine the position of the 
microscope and its focal plane with respect to the patient's 
head and CT scans given the distances between spark gaps and 
microphones." The purpose of this section is to determine 
the appropriate placement and orientation of the spark gaps 
given the constraints of the ultrasonic digitizer. 

As mentioned in Section 4.2.4.1, at least six pieces of 
information must be known about a rigid body (microscope) to 
locate it in a fixed coordinate system (CT). In Section 
4.2.3 it was determined that the coordinates of three points 
(nine pieces of information) are easier to determine than 
relative angles and sufficient to locate the position of the 
microscope. Therefore, three spark gaps are needed in 
addition to at least three microphones. 

Physical and aseptic constraints in the operating room 
would prevent mounting either the microphones or the sparks 



i. 
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gaps on the patient's head to determine the relative 
rnicroscope-CT scan positions. These constraints also 
indicate that counting the microphones (connected to 3x5x13 
centimeter pre-amp boxes) on the operating room ceiling 
between the light tracks and air vent, and mounting the 
spark gaps on the microscope would probably be the most 
acceptable design. 

The orientation and separation of the spark gaps are 
important as errors in the slant ranges can produce errors 
at the focal point. The optimal spatial distribution of 
three points would be an equilateral triangle since random 
errors at each point are more likely to cancel and reduce 
the error at the geometric center than the errors resulting 
from other distributions [3j. For example, three points in 
an isosceles triangle would magnify the error along the line 
from the apex bisecting the base. The optimal spark gap 
separation (5) and distance (D) from the geometric center of 
the spark gaps to the focal point (f) can be determined for 
an error at each point (e). If the tnree spark gaps lie on 
an equilateral triangle in the x-y plane (See Figure 5.1) 
and the focal point lies in the z direction, the maximum 
error at F in the x-z plane (e xz >j given a maximum error at 
each point (e), can be determined by similar triangles. See 
Figure 5.2. Point f' is the new position of the focal point 
due to e. 
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Since angle 9 can be considered very small, these error 
vectors can act perpendicularly to each other. Therefore, 
their magnitude is given by: 

e = Ve 2 + e 2 , 
max * xz vz 

e = Jno 2 e 2 /S 2 * Ue^CD* + S^/12)/3S 2 , 
max T 

e = e-IWD/S) 2 + 1/3. 
max ' 

This relationship shows that the minimum error at the 
focal point occurs when S is a maximum and D is a minimum. 
This makes sense, as minor microscope movements will produce 
greater spark gap displacements for larger values of D and, 
therefore, produce greater effective resolution in 
determining the location of the focal point. The spark gap 
separation, however, is limited by the required mobility of 
the microscope. The sterile draping and physical 
constraints of the operating room limit the spark gap 
separation to approximately 30" centimeters. Minimizing the 
lever arm, D , to zero also cannot be achieved because of 
physical operating constraints. Depending on the position 
of the patient, two different distances are practical.- With 
the patient in a sitting position the lever arm can be no 
less than .approximately 45 centimeters and with the patient 
in a supine position no less than 25 centimeters. 

The general registration procedure involves CT scanning 
the patient with at least three fiducials attached to the 
head. Aiih the spark gaps mounted on the microscope and 
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microphones mounted on the ceiling, the patient 1 s head will 
be placed in a rigid clamp which anchors the position of the 
fiducials with respect to the microphones. See Figure 5.4. 
The positions of the microphones with respect to the CI 
• scans are then determined by focusing on each fiducial (see 
Section 4.2.J) and will be discussed in detail in Chapter 6, 




Fi pure 5. 4 
Fiducial - Microphone Relationship 



Chapter 6 - Reference Frames 

This chapter can be divided into two major parts. The 
first part describes three possible reference frames 
considered for the calculations of the reconstructed CT 
iraage at the focal plane - the CT, digitizer, and microscope 
coordinate systems. Each of these offer advantages and 
disadvantages which are outlined in sections 6.1 - 5.3. The 
second part (section 6.4) details the reference system 
design including the equations and algorithm references 
necessary for determining the spark gap - focal plane 
relationship (5.4.1) and the conversion of focal points to 
CT coordinates (6.4.2). Also included in section 5.4.1 is 
the description of a fourth reference frame, an oblique 
spark gap coordinate system. Section 5.4.3 describes the 
general reference system procedure. 

6.1 CT Coordinate System 

The CT slices are stored in a CT Cartesian coordinate 
system established by the CT scanner. The z axis 
corresponds to the transverse body axis for conventional 
scans where the CT gantry is zero degrees from the vertical, 
The CT scans are slices in the'x-y plane. See Figure 6.1. 

The advantage of performing reference and 
reconstruction calculations in the CT coordinate system is 
that all the information (CT data) is in that reference 
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Figure 6.1 - CT Coordina te System 



frame, requiring few transformations between coordinate 
frames, which can take a lot of processor time and introduce 
error to the system. The disadvantage of working in CT 
coordinates is that it restricts the reference system to one . 
that, in practice, may be awkward or undesirable, further 
investigation will determine if this is true. 
6.2 Digitizer Coordinate System 

A digitizer coordinate system can be established to • 
determine the coordinates of the spark gaps with respect to 
the microphones. As mentioned in Section 4.2.5, 
triangulation of the slant range data given three 
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orthogonally placed microphones whose separation is known 
determines the Cartesian coordinates of the spark gaps with 
respect to the microphones. See Figure 6.2, The 
microphones do not necessarily have to be orthogonal; an 
oblique coordinate system can also be established. 

The spark gaps (therefore, the microscope and focal 
plane) are repositioned within the digitizer coordinate 
system established by the microphones on the operating room, 
ceiling. The patient's head and CT scans qre 3lso within 
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this space, which may make calculations in the digitizer 
reference frame easier to understand. The disadvantages of 
working in digitizer coordinates are that all the CT data 
must be transferred to the digitizer system and the exact 
microphone separations >nust be known. 

6.j Microscope Coordinate System 

The reconstructed CT scan must eventually be displayed 
in the microscope as a two dimensional CRT image. This 
involves converting the reconstructed slice from a matrix of 
three coordinates to one of two coordinates, (x,y). A 
microscope coordinate system could represent the focal plane 
as x and y, normal to the optical axis with. the origin at 
the focal point. See Figure 5.3. This technique requires a 
transformation of coordinates because the microscope 
coordinate system will be constantly changing with respect 
to the location of the microphones and CT scans as the 
surgeon moves the microscope. Regardless of the reference 
frame used for reconstructing the slice, in order to display 
the proper image, the slice must be transformed into 
microscope coordinates. The advantage of transforming all 
the data and performing the calculations in the microscope 
coordinate system is that if the surgeon only moves the 
microscope slightly or along the optical axis, the 
reconstruction calculations could be greatly reduced and 
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Figure 5.3 
Microscnpp coordinate Sustem 



allow the surgeon to quickly call up successive new slices 
for viewing. 

5.M Reference System Design 

This section describes an optimal reference system 
design and the goals of such a system. 

The number of CT slices and whether or not grey scale 
information will be processed (only contours will be 
reconstructed initially) will determine the appropriate use 
of che CT coordinate system. Since the Cf data will be 
manipulated by the Treatment Planning Computer, it is 
probably best to perform the reconstructions in the CT 
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coordinate system. Since we know that the reconstructed 
slice must be displayed in the microscope, the microscope 
coordinate system must be established for transforming the 
new slice at each microscope position. A digitizer 
coordinate system is unnecessary because the positions of 
the microscope, sparks gaps and microphones can be 
determined in CT coordinates. 

The general reference system design will involve 
determining the position of the focal plane in CT 
coordinates. Since this is an initial design the 
information needed by the Treatment Planning Computer to 
reconstruct and display a CT image- can be divided into three 
parts'. 

1) The equation of the focal plane in CT coordinates is 
necessary to determine the appropriate reconstructed image. 

2) The CT coordinates of the focal point must be 
determined such that the center. of the reconstructed • slice 
will be properly displayed. 

1) Since the equation of the focal plane and the CT 
coordinates of the focal point do not uniquely define the 
orientation of the focal plane about the optical axis, more 
information is needed. The three direction cosines of the 
Y M axis of the microscope coordinate system with respect to 
the X--. Y„ T and Z rr axes will provide the necessary angular 
information to orient the reconstructed Cf image. See 
Figure 5.4. 




6.4.1 Spark Gap - Focal Plane Relationship 

The patient is anchored to the operating table 
generally in one of two positions - sitting or supine (see 
section 5.4). The physical and sterile constraints of the 
operative procedure restrict the placement of the three 
spark gaps such that two spark gap holder brackets are 
necessary, fixing the position of the spark gaps with 
respect to the microscope and focal plane in each case. 
Therefore, the first step in the reference procedure is to 
detertaine the relative positions of the spark gaps and the 
focal plane for both the sitting and supine cases. 
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The information needed to calculate the equation of the 
focal plane, the CT coordinates of the focal point, and the 
three direction cosines is explained below. 

Three simple ways to define a plane include 
determining.: 1) the coordinates of three non-colinear 
points in the plane, 2) a normal vector to the plane and one 
known point in the plane, or 3) the coordinates of two 
points along a normal line, one of which is in the plane. 
In order to uniquely define the focal plane, we also need to 
know the orientation of the pline about the optical axis 
(focal point). Therefore, the coordinates of at least three 
points are needed to define and orient the focal plane in CT 
coordinates (X rT Z rT ) : the focal point, a point along 
the optical axis (normal point) and a point in the focal 
plane (orientation point). See Figure 5.4. 

6.4.1.1 Focal Point 

As the first, step this section introduces an oblique 
spark gap coordinate system and describes the equations 
necessary -to determine the focal point with respect to the 
spark gaps. 

One technique is to calculate the coordinates of the 
focal, normal and orientation points in an oblique 
coordinate system defined by the vertices of the equilateral 
spark gap triangle (explained in section 5.4). See Figure 
6.5. The locations of the focal, normal and orientation 
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points ;nust be determined with respect to the fixed spark 
gaps, independent of the microscope's position. First we 
will define the oblique spark gap coordinate system. The 
origin is at spark gap 1 (SI),- the X s axis is along the line 
connecting S1 to spark gap 3 (S3), the * s axis is along the 
line connecting SI to spark gap 2 (32), and the Z s axis- is 
along the line from the SI perpendicular to the spark gap . 
plane. The obliquity angle, w,- can be calculated by 
applying the Law of cosines to the measured spark gap 
separations S 12 , S 2 ji and S^: 

cos(w) = (S 12 2 + S 31 2 - S 23 2 )/2S 12 S jl (6.1). 
As will be shown below the oblique spark gap 
coordinates of the focal point can be calculated by focusing 
the microscope on the tip of an additional spark gap. With 
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the microphones on the ceiling, the digitizer can determine 
the slant range distances from the four spark gaps to the 
three microphones. The first task is to determine the 
location of the microphones in the oblique spark gap 
coordinate system. This can be done based on the slant 
range distances from the three microscope mounted spark gaps 
to the three microphones. Once the microphone coordinates 
in the spark gap coordinate space are known, the focal point 
can be determined. The appropriate equations are found as 
follows. 

Figure 6.6 shows the relationships between the slant 
range distances (SL^ and the spark gap coordinates -p(x s , 
y , z ) for one microphone. The perpendicular projections 

s s 

^ or <-u~ v y avp<* (P P ) are determined by 

of SL -j on the and axes u x , r y' 

applying the Law of cosines: 

p = 5L. cos( a 1 ), (5.2) 
x 1 * ' 



x 

x = w 31 T ^1 

Similarly, 



P v = (S 2 «■ SL 2 "- SL 2 )/2 Sj1> (6.3) 



Py = (S 12 2 ♦ SL, 2 - SL 2 2 )/2S 12 . (6.4) 
Coordinates x c and can be calculated from: 

S w 

cosCw) = (P x - x s )/y s = (P y - y s ) /x S ' (6 - 5) 



therefore , 

x = 
s 

(p - x ) /cos( w) • (6.7) 



x = (P - y ) /cost w) , (6.5) 
s y s 



Figurs 5. S - Spark Gap Coord inatos 



Substituting Co. 7) into (5.6): 

x = (P - P cos(w))/sin (w) . (5.3) 
s x y 

Alternately, substituting (5.6) into (5.D yields: 

y = (p - P cos(w) )/sin 2 (w) . (6.9) 

z' can be found by applying the Law of cosines: 

z' 2 : x 2 + y 2 - 2x y cos( n - w) (6.10). 
u s 3 s s s 

The coordinate z, can then be found by applying the 

s 

Pythagorean theorem: 

z 2 = SL 2 - z' 2 (6.11), 

5 I 

z s = VsL 2 - x 2 - y 2 . - 2x y,cos(w) (5.12) 
Therefore we have determined the coordinates of one 
microphone in the spark gap coordinate system. ihe 
coordinates of the other two microphones are found 



similarly. In summary, the equations needed to calculate 
the oblique spark gap coordinates of the three microphones 
are : 

1) apply 6.1 to determine w, 

2) apply 6.3 and 6.4 to determine -? x and P y , 

3) apply 6.3 and 6.9 to determine x s and y s , 

4) apply 6.12 to determine z g . 

Once the oblique spark gap coordinates of the three . 
a,icrop,io»ea are deter-ninea, the cooruiflaces of one fouai 
point (spark gap 4) can be calculated by solving three 
nonlinear equations for the distances between each 
microphone and the focal point. See Figure 5.7. M 1 , \, 
and .< represent microphones 1, 2, and 3; F is the focal 
point whose oblique coordinates Cx a ,y s ,z s > are unknown; and 
D , D 2 and D 3 are the slant range distances between F and 
M lf M 2 and A y The equation relating the distances, D., 
between F and the microphones in oblique spark gap 
coordinates (for i = 1 to 3) is: 

D i 2 s ( V x s. )2 + (y s- y s. )2 + (z s- z s. )2 

+2(x -x )(y,-y, )cos(w). (6.13) 
s s i s s i 

Solving these three equations by an iterative Newton's 
.method (see Chapter 9) for x s ,y s ,z s will give the oblique 
coordinates of the focal point in the spark gap coordinate 
system . 
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5.7 - Focal Point Ca I cu I at.i ons 



6 . M ♦ 1 .2 formal Point 

This section describes the equations and procedures 
needed to determine the. oblique spark gap coordinates of the 
normal point along the optical axis . The first -step 
involves establishing a Cartesian coordinate system 
(indicated by subscript F) that is fixed with respect to the 
microphone, by focusing the microscope on three points on a 
piece of plexiglas where the coordinate's of the three points 
are known. This will allow both the normal and orientation 
points (see section 6.4.1.3) to be determined independent of 
the oblique spark gap coordinate system. Once this is 
known, the fixed coordinates of the focal point can be 
calculated for any microscope position since its location is 
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known with respect to the spark gaps (oblique spark gap 
coordinates) and the location of the spark gaps are 

known with respect to the microphones (digitizer slant range 
distance). Assuming that the focusing knob of the 
microscope accurately moves the focal point along the 
optical axis, the fixed coordinates of a" normal point can be 
determined and then its coordinates in the oblique spark gap 
system referred to the original focal point. 

To determine the coordinates of the three microphones 
in the fixed coordinate space we rely on the fact that the 
distances between the focal point and microphones can now be 
calculated. Refer to Figure 6.3. The oblique spark gap 
coordinates of the focal point ( x fp ,y f p ,* f p ) are known and 
fixed for each spark gap bracket .(sitting or supine), and by 
applying the steps outlined in 6.4.1.1 (Equations 6.3, 6.9 
and 5.12) the oblique spark gap coordinates of each 
microphone Cx^.zJ can be calculated given the three 
slant range distances SL, , SL 2 and SL3. The distances 
between each microphone and the focal point can be 
determined' by applying the oblique distance formula 
(Equation 6.14) 

° 2 = ( vV 2 * (y nfV 2 * c vV 2 + 

2(x m- x fp )(y m- y fp )c0;3Cw) (6 - 14) ' 

3y focusing on three points (P 1 ,P 2 ,P J in Figure 5.9) 
whose coordinates are known in the fixed coordinate system 
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Figure 5.3 - Focal Point - Microohnng Distar.cs 



and calculating each point - microphone distances O n - 
D 3j ) using Equation 6.14, the fixed coordinates of the 
microphones (M,,^,,^) can be determined by solving three 
nonlinear equations for each microphone. J^-M-j represent 
the microphones, P r P 3 represents the points in a fixed 
Cartesian coordinate system, and u^-D^ represent the 
distances calculated between microphone i ( x m 'ym^m^ and 
P^P The fixed coordinates of microphone i can be 
calculated by solving the three (j = 1 - 3) nonlinear 
equations' belbw using iJewton's method given the coordinates 
of Fj and distances D^: 



63 




2 <x 2)5 j a ,2y> 



.'Via ^li^Ma 



Fiaura 6.9 ~ Micraphona Coordinates 



<v^ )2 *'v y j )2+( v i J )2 " ° ij2 (5 ' 15) " 

Once the fixed coordinates .of each microphone are 
known, the fixed coordinates of the focal point can be 
determined independent of the microscope's position. In 
order to calculate the fixed coordinates of a normal point 
the oblique coordinates of the three microphones 



z ) must be determined by applying Equations 



(x ms. ,y ms i , 'ms i 
6.8, 6.9 and 6.12 with the microscope in a fixed position. 

Then by adjusting the focusing knob and moving the 

microscope along the optical axis, as shown in Figure 6.10, 

the spark gaps will now be at points ST, S2« and S3' and 
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the new focal point CF 1 ) will be at the normal point, N. 

The fixed coordinates of N (x ,y ,z ) can be found by 

n n n 

solving Equation 5.14 for the normal point - microphone 

distances, D i , and then solving the three (j=1 - 3) 

nonlinear equations by Wewton's method for x , y^ and z 
^ J n J n n 

given D. , and x m , y m and z„ : 
i i i i 

(X -X m ) 2 +(y -y m ) 2 +( z - z ) 2 = D 2 (5.16). 

n m^ n zik n l 

jp to here we know the following: 

1) the r i y o j coordinates of the three i - r c p h - n e s — 

(x m. ,y :n. ' z m^ 9 
ill 

2) the oblique coordinates of the three microphones 
with respect zo the spark gaps (and microscope) at SI, 52 
and S3, Figure 5.10 - (x,,,^,,,^^), 

3) the fixed coordinates of the normal point, M - 

(x.y n ,z ) . 
n rf n n 

The oblique coordinates of 'A ^ x ns i y ns » z ns ^ with the 

spark gaps at SI, S2 and Sj of Figure. 6.10, can be 

calculated by solving the three nonlinear equations by 

Newton's method for x ns> y ns and z ns , given the normal point 

- microphone distances D^^ (from Equation 6,15) and the 

oblique microphone coordinates (x ,y ,z ): 

cos • ms. ins. 
Ill 

(x -x ) 2 + (y -y ) 2 +(z -z ) 2 + 
ns ms. ; ns ; ais. ns ms. 

i l 2 i 

2(x n<r x ,nc )(y„-y fl . )cos(w)=D * 16. H) 
ns ms. 'ns'ms. i 
l i 
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in summary, the steps needed to calculate the oblique 
spark gap coordinates of the normal point are: 

1) determine the coordinates of the three microphones 
in a fixed coordinate system by focusing the microscope on 
three points whose fixed coordinates are known, and applying 
Equations 6.3, 6.9, 6.12, 6.14 and 5.15, 

2) determine the fixed coordinates of point W by moving 
the microscope along the optical axis ( focusing "knob) and^ 
applying Equations 6.14 and 6.16, 

3) calculate the oblique coordinates of point A by 
solving Equation 6.17 for x ns , y ns and z ns . See Figure 
6.10. 

6.4.1.3 Orientation Point 

This section describes the steps and equations 
necessary to calculate the oblique coordinates of an 
orientation point, 0, along the Y.< axis of the microscope 
coordinate system (see Figure 6.4). This point, in 
combination with the focal and normal points will determine 
the orientation of the focal plane about the optical axis by 
indicating the "top" of the displayed image as seen through 
the microscope. The ocular crosshairs of the microscope are 
instrumental in focusing the microscope, and determining the 
microscope coordinate system. See Figure 6.11. 

The oblique coordinates of the orientation point can be 
found by following essentially the same steps in determining 
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Figure 5.13 ~ Normal Point 




Figurg B. 11 - Orientation Point 
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the normal point. The normal point is only used to 
determine the equation of the focal plane and its distance 
from the focal point •( although limited by the range of the 
focusing knob) is arbitrary as long as it falls along- the 
optical axis. Likewise, the distance between the focal and 
orientation point is arbitrary as long as the line 
connecting them is 90° to the optical axis of the microscope 

coord mate system . 

To take advantage of the fixed coordinate system 

a u 1 ? we will assume the fixed 
established in section o.M.l.^, we 

coordinates of the microphones are known C^.,^..^. 5 and 
that the oblique coordinates of the microphones 

r x y z ) have been determined such that the fixed 

ms ins ' hi s 

coordinates of\he focal point (F) and a point (P) located 
along the f, axis of the microscope coordinate system are 
also known. See figure 6.12. As with the normal point, the 
fixed coordinates of an orientation point must be determined 

first . 

This can be accomplished by projecting the known point 
P (x ,y ,z ) onto the focal plane and determining its new 
fixed' coordinates. If the microscope is focused with 
minimum parallax error such that the Y„ ocular crosshair is 
closely aligned along the line between F and P, then the na 
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projected point can be considered the orientation point. 

See Figure 6.13. 

in order to determine the fixed coordinates of 0, 

(» y i ) the equation of the focal plane must be 
established. The normal vector alon g the optical axis can 
be represented in the fixed coordinate system by the 

equation : 

N = At f «■ B} f ♦ Cfc f (6.18) 

where * fi J f and ft f are the unit vectors along the fixed X, 
t and Z axes. The coefficients A, 3 and C can be determined 
fro. the fixed coordinates of the noraal (x n , y n and z fl ) and 
focal (x fp , y fp , z fp ) poi^s. 

N = ( V x f P }1 f + Cy n" y fP ,J f * ^"VV 

(6.19) 

Since N is a vector nor-.l to the focal plane, then for any 
vector T = (X f ,i f ,Z f ) in the plane: 

N-CT - V f ) = 0 (6.20) 

where T f is the focal point vector Cx fpl y fp .* f p > ■ 
Therefore, 

(v»fp )a f- x fP ) ^ y h- y fp )af " yfp><Un "' fpMZf ' Ifp) ' ° 

(5.21) 



to 



is the equation of the focal plane in the fixed coordinate 
system, in order to calculate the fixed coordinates 
(v y o , Zo ) of the projected point P(x p .y p ,* p >, the equation 
of 0 the line through P normal to the focal plane must be 
determined. It is given by: 

where A, 3 and C are the normal vector coefficients. If we 

let this common ratio equal t: 



x = tA + x , (5.23) 
o P 

y = ta + y_, (6.24) 

J o P 

z = tC + z , C5.25) 
o P 



which are the parametric equations of the line. 

intn the eauation of the focal 
Substituting x Q , y Q and z Q into tne equ* 

plane where D is a constant: 

Ax + 3y + Cz = D, 
A(tA+x p ) + 3U3+y p > +■ C(tOz p ) = D, 
t= (D-Ax p -By p -Cz p )/(A 2 + 3 2 + C 2 ). (6.26)' 

Substituting Equation 6.25 into 6.23, 6.24 and 5.25 yields 
the fixed coordinates of the projected orientation point 
(x o ,y of z 0 ) CU53. The oblique coordinates of 0 with the 
spIrkVps at SI, S2 and S3 of Figure 5. TO, can be 
calculated by solving tne three nonlinear equations by 

* - „ and z ziven the orientation 

Newton's method for x QS , y os ana z os S1 
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point - microphone distances (from Equation 5.15) and the 

oblique microohone coordinates (x ,y ,z ): 

•n s ^ 7 ' m s ^ 1 m s ^ 

(x os- x r flSi )2+( yos-ym Sl )2 + (z o S - z tn Si )2 

1 1 

In summary, the steps needed to calculate th-e oblique 
spark gap coordinates of the orientation point are: 

1) determine the fixed coordinates of point 0 by 
projecting point ? onto the focal plane using Equations 
5.22-5.26, 

2) calculate the oblique coordinates of point 0 by 
solving equations 6.27 for x , y • and z_ . 

OS OS OS 

The relationship between the oriented focal plane and 

the spark gaps can now be defined by the oblique spark gap 

coordinates of the focal, normal and orientation points, 
independent of the position of . the microscope. 

6.4.2 Focal Points in CT Coordinates 

Now that the relationship between the spark gaps and 
the focal points has been established, the general reference 
system procedure must be described. 

6.4,2.1 Registration Procedure 

The next step is to relate the CT scans and CT 
coordinate system with respect to the microphones and, 



therefore, the focal, normal and orientation points. As 
mentioned in Section 4.2.3, three CT-detectable fiducials 
can be. attached to the patient's head and their coordinates 
determined in the Cartesian CT coordinate system. 3y 
applying the principles of Section 6.4.1. 2, the fiducials 
would represent P,, ? 2 and ? 3 and the CT axes would 
represent the fixed axes of Figure 5.9. The coordinates of 
the microphones can be found in the CT coordinate system by 
focusing on each of the fiducials and solving Equations' 
5.15. At this stage we assume the patient's head will reualn 
fixed during the rest of the procedure, e.g. that the 
patient is anesthetized and anchored in the head clamp 
(fiducials fixed with respect to the microphones). 

6.4.2.2 Focal Plane Determination 

With the microphone coordinates known, the neurosurgeon 
can- perform the craniotomy and bring in the microscope when 
needed. Once the microscope is focused, the digitizer and 
an IBM PC XI computer will determine the slant range 
distances and the coordinates of focal, normal and 
orientation points calculated in the CT coordinate system. 
This is achieved by solving (i = 1 to 3): 

(x -x ) 2 + (y -y ) 2 *Cz p -z, ) 2 = Dp 2 ( 5 - 23) 
U P V ; y P y m i p m i *i 

for x p) y p , z p , the focal, normal and orientation points, 
using'! Newton's method, given D p (the distances from each 
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point to the jicrophones; , which is determined by the 
oblique distance formula. 

The equation of the focal plane can be determined by 
applying Equation 6.21, and the coordinates of the focal 
point are known such that the reconstructed CT image is 
appropriately centered. The orientation of the focal plane 
about the focal point must now be calculated. The three 
direction cosines of the vector from the. focal point to_the 
orient^-ion point will p-oviie e-ou-h infomat ion . See 
Figure 5.14. The vector in the CT coordinate space is: 

. t = (x 0 -x fp )£ + Cy 0 -y fp )J+ <* 0 - z fp )ft (5 ' 29) 
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Figure 5. 14 - Direction Cosines 
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Step 2 CT Scanning 

CT scan the patient ,1th three CT-detectable fiducials i» 
pla ce and determine their coordinates in the CT coordinate 
system. 

Step 3 Registratio n Procedure 

Anesthetize or sedate the patient and anchor his head in the 
head clamp. With the soar* gaps in place on the microscope, 
focus on each of the fiducial, and determine the coorui.ates 
of the microphones (Section 5.4.2.1). 

Step 4 Focal Plane Determ ination 

Remove the fiducials if necessary and perform the 
craniotomy. Focus the microscope and calculate the 
coordinates of the focal, normal and orientation point. 
Determine the equation of the focal plane and the direction 
cosines . 



rha pter 7 - F rrnr Analysis 

This chapter reviews all the major sources of error for 
the reference - display system. The purpose of such an 
analysis is to determine the theoretical worst case 
magnitude of error that could be expected and to see if this 
error can be tolerated or reduced if found unacceptable. 
The techniques used to reference the focal plane and CT 
scans must be determined, such that they propagate the least 



error . 



The questions this chapter must answer are: What 
techniques should be used and what parameters should be 
optimized to obtain a maximum error preferably less than or 

equal to 1 millimeter? 

The purpose of the reference system is to determine the 
location of the focal plane with respect to the CT 

There are four main. sources of error for 
coordinate system. There are 

c^m. ^) error in the slant range distances 
the reference system: U error 

t-r -n error in focusing the 

from the digitizer (7.1), 2) error 

,7 ^ V CT scan and display error (7.3), and 4) 
microscope (7.2), 3) scan 

(i m Miscellaneous 
propagation or computation error C7.«U. 

errors are discussed in section 7.5. The limits on each of 
these errors will be determined in each section, and the 
total error evaluated in section 7.6. 
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7.1 Digitizer Error 

The errors in determining the slant range distances can 
be categorized as either errors in accuracy or precision. 
According to S.A.C., the accuracy of the slant range is 
within +0.1* (bias error) of the measured distance and the 
precision (repeatability) is within +3.01 centimeters. 

The parameters that affect the speed of sound in air, 
and, therefore, the bias error were discussed in Chapter 5_. 
The t;;ree aajsr factors are: temperature variatlar.J CAT), 
relative humidity (<£) and air motions represented by a. 
Temperature variations and relative humidity cannot be 
easily controlled, but air disturoances that will distort 
the sonic wavefront and alter the digitizer reading can be 
minimized by not opening the operating room door while 
digitizing, etc. The following equation describes this 
relationship: 

T(d, AT, i, a) = k(AT, <f>, m)d 

where e* is the slant range distance vector error (in cm.) 
a 

due to the accuracy limitations of the digitizer, K is the 
bias error percentage (worst case +0.1*), and d is the slant 
range distance (in cm.). T a is represented as a vector 
because this error acts along the line between the spark gap 
and microphone. 
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oased on the fundamental analysis of Chapter 5, the 
speed of sound is affected most by temperature variations. 
If the bias error percentage, k, given by S.A.C., is based 
only on temperature, a temperature variation, AT, of 0.6°C 
corresponds to a k of +0.1>. Temperature differences from 
the operating room ceiling to the operating table were 
measured during a typical surgical procedure and the mean AT 
was experimentally found to be 0.2°C. Since k is a linear 
function of the temperature variation, experimentally k = 
+0.0343 for the operating room environment. The discrepancy 
between our experimental value of k and the worst case value 
of 0.1* indicates that the operating room is a well 
controlled atmosphere for the operation of the digitizer and 
may produce smaller errors than assumed by S.A.C. 

The slant range distance error in precision is a 
function of the digitizer's clock cycle frequency and is 
fixed at +0.01 centimeters. 

e* = 0 .01 cm. , 
P 

where eT p is the precision slant range distance error. 

The worst case slant range distance error for. each 
spark gap - microphone path (e^ in cm.) can be expressed as 
the sum of the worst case errors in accuracy and precision: 

e* d = k(AT, m)d + F p , 
= 0 . 1 id + 0 . 01 cm . 

I 
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The worst case digitizer error (e" d ) , again, is a vector 
because it acts along the line between the spark gap and 
microphone. Its impact on the magnitude of the total 
digitizer error (e td ), therefore, depends on the direction 
of the spark gap - microphone line. An estimate of this 
worst case error for three slant ranges acting at typical 
angles is 

= 2(0. Ud + 0.01 ) cm. 

td 

Reducing d will reduce e t(J . 

One technique to reduce digitizer errors due to random 
disturbances, such as air motions, is to average multiple 
samples. Since the digitizer can determine approximately 30 
slant range distances per second (transmitting data at 9600 
baud), averaging several seconds of data could reduce the 
error that would result with one sample by a factor of 1/ n, 
where n is the number of samples for a normal distribution. 
[13] This technique is explained in the software of Chapter 
9. 

7.2 Focusing Error 

There will be a certain amount of error in focusing the 
microscope, partly due to optical limitations and partly due 
to human error. A quantitative description of the 
microscope's optical properties were presented in Table 2.1. 



There are two types of focusing error: errors in depth of 
field and transverse error in the focal plane. The larger 
error normally comes from the depth of field. From Table 
2.1 it is observed that an increase in magnification results 
in a decrease in depth of field. The smaller the depth of 
field, the smaller the focusing error. Therefore, the 
highest magnification should be used when registering the 
microscope. 

At the. highest magnification (microscope scale factor 
2.5) the depth of field is two millimeters and the optical 
field size is 1.45 centimeters. The transverse error is 
determined by . the ability to focus non- intersecting 
crosshairs on the focal point within the focal plane. At 
the highest magnification the distance between the 
crosshairs is approximately one millimeter. See figure 
At best, an untrained observer can focus the microscope 
within a volume of a cylinder with a transverse diameter of 
one millimeter and a height (depth of field) of two 
millimeters. The surgeon, however, can probably reduce this 
error through experience with the microscope to be within 
one millimeter in all directions. 

Two other factors affect the focusing error - parallax 
error and microscope balancing. Since the crosshairs are in 
only one ocular, there will be a small amount of parallax 
error in the visual interpretation of the focal point. If 
the microscope is not. balanced properly in all planes, it 
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Figure 7. 1 - Microsco pe Crosshairs 

will be difficult to focus on a specific point because the 
microscope will drift after the locking mechanism is 
engaged . 

The focusing error is, therefore, a function of 
microscope magnification (depth of field), surgeon 
experience with the microscope, parallax error and errors 
due to an improperly balanced microscope. 

Based on an estimate of the above errors, the worst 
case focusing error is on the order of one millimeter with 
properly balanced microscope and an experienced surgeon. 

T f = 1.0 mm , 

where e" f is the worst case focusing error written as a 
vector because it can act in any direction. 
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f. j CT Scan/Display Error 

Error can also result from the CT scanning procedure, 
the fiducial composition, size, shape and placement and the 
reconstruction/display calculations and procedure. These 
sources of error are described in this section. 

J . 1 . 1 CI Resolut ion 

The head CI scanning field covers a 25 x 20 centimeter 
rectangular area and, since tne scanner resolution xs * 
320 voxels, the measured resolution in the center of the 
scan in the x-direction is 200mm/320 voxel s or 0.5 mm/voxel 
and in the y-direction 250mm/ j20vox el s or 0.3 mm/voxel. The 
minimum voxel size in the z-direction (body axis) is the 
minimum slice thickness, 1.5 mm. Unfortunately, practical 
constraints on the use of the CT scanner prevent such 
step-by-step scanning at 1.5 millimeters per slice. A 
maximum of 20 slices are usually taken at least 5 mm thick 
with no overlap, so the voxel z dimension is typically 5 mm. 
The worst case CT resolution error occurs when a fiducial is 
calculated to be at the center of a voxel and it is actually 
at one end. 

|Tj = V(0.6mm/2) 2 + (0.8mm/2) 2 + (5mtn/2) 2 = 2.5 mm. 

This vector error, from a worst case point of view, prevents 
the specification of total system error less than 1 
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millimeter, but it can be reduced as explained in the next 
section . 

7 .'5.2 Fiducial Error 

Since the CT detectable fiducials will link the 
microscope position to the CT scans, they can also 
contribute to the error of the reference system. Their 
careful design is, therefore, important and will be 
discussed in this section, 

7.3.2.1 fiducial Composition, Siie, and Shape 

There are several constraints on the fiducials that 
will determine their composition and size, their ability to 
be CT scanned and, therefore, the error that they will 
impose on the reference system. Since a CT slice is 
actually a density map, the fiducial composition must be 
dense enough to appear clearly in the CT scans. Metals meet 
the density criteria, but their repetitive molecular 
structure causes x-ray scattering that distorts the CT scan. 
However, amorphous materials with irregular molecular 
structures such as glass will not cause scattering and they 
can be doped with lead to increase the density. Glass beads 
of 2 and 5 millimeters in diameter were CT scanned, but only 
the 5 millimeter sample was clearly seen. Glass appears to 
be a suitable fiducial material. 



fhe shape and size of -trie fiducial is important for two 
reasons. It must not be large and uncomfortable since it 
might be attached to the patient's head for up to twenty 
four hours. It also must be large enough to be visible to 
the surgeon, as- well as small enough to define a unique 
point in the CT scans. 

Since the thinnest practical CT slices are 3.5 
centimeters thick, any fiducial smaller than that could be 
detected in the scan, but its position within the slice 
would be unknown. The shape of the fiducial ::rjst x ze 
irregular to overcome this problem without resorting to an 
uncomfortable .0 . 5 centimeter diameter glass bead. Two- 
different designs are proposed: 

1) a glass wedge whose CT scan density would vary 
depending on the position of the slice. See Figure 1.2. 

2) a glass "W'-shape that is used with the 
stereotactic frame and defines the scan position by the 
distance between the elliptic intersection of the CT slice 
and the fiducial. See Figure [.'1. 

Since the glass wedge would be- difficult to make 
precisely, the M N"-shaped glass fiducial was chosen to be 
made approximately 2 centimeters square out of 2 millimeter 
diameter glass rod. The "N" must be oriented such that the 
CT scan intersects both the parallel posts and the 
crosspiece. Initial work, however, will be conducted with 
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the 5 millimeter glass beads," and the construction of an ,! N M 
shaped fiducial left for future work. 



7.J.2. 2 fiducial Placement 

The spatial distribution of the required three 
fiducials, like that for the spark gaps, is important in 
reducing error. The spatial distribution contributing the 
least error would be to place the fiducials at the vertices 
of an equilateral triangle (see Section 5.4) [31. Since it 
might be uncomfortable to attach three fiducials in such an 
orientation on the patient's head, the proposed sites would 
ue one temples and fjreneau. 

f.3.2.j Total Fiducial Error 

The ability of the fiducials to define the CT 
coordinate system based on CT scans and register the 
position of the microscope based on their visibility is 
primarily a function of two variables, as indicated below: 

^Vid tr r r r } 

where eTf^ is the vector error due to the fiducials, e^ is- 

the focusing error, and e* is the CT resolution error. 

r 

Since the M N ,! -shaped fiducial will determine the center of 
the 0.5 centimeter CT slice, the CT resolution is also a 
function of the fiducial shape and is, therefore, much less 
than the worst case value of 5.1 millimeters given in 
Section 7.J.1. The worst case fiducial error must be 
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determined experimentally before a quantitative value can be 
assijned, and will be left as future work. 

7 . j . 'j Display/Reconstruction Error 

There is one major source of error that comes under 

this category - errors due to reconstructing or 

interpolating the CT slice (e^) or contour. It is based on 

the number of CI slices (n ), the slice separation (sep), 

s 

the angle between the reconstructed slice and the serial 
slices (9) (hence, the number of points common to both the 
reconstructed and serial CT slices) and the type of 
interpolating algorithm used (A). 

The most important factor in interpolation is the 
number of data points - the more points the smoother and 
more accurate the reconstructed curve. The number of data 
points depends on the slice separation and the number of 
slices. Therefore, the error between a CT image 
reconstructed from many CT scans and a CT slice produced by 
the CT scanner directly can be represented as: 

e\U, 9(n s ,sep)) , 

where e^ is the interpolation error. Again, the effect of 
this error must be determined experimentally and will be 
left as future work. 



93 



f.4 Computation irror 

The errors due to digitizer calculations, focusing, and 
CT scan reconstruction and display can be magnified as they 
propagate . through the calculations necessary to reference 
the focal plane and CT scans. The effect of this 
propagation depends on the calculations performed and the 
orientation of the vector errors, and will be discussed in 
detail in Chapter 9. This-error will be referred to as 

U COillp" 

1.5 Miscellaneous Errors 

There are three smaller and probably negligible sources 
of error that should be mentioned. 1) There is a 
possibility that the operating table is disturbed after the 
registration procedure, which will relocate the position of 
the CT scans with' respect to the digitizer microphones and 
contribute to the error of the reference system. 2) Since . 
the fiducials will be attached to the patient's skin, there 
may be some relative skin - brain motion between CT scanning 
and the registration procedure, j) There may be computation 
error in transforming the coordinates of the reconstructed 
scan to the microscope coordinate system, which would result 
in some display error without careful programming using 
floating point arithmatic. 



7.6 Conclusions 
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The main purpose of this chapter was to determine the 
sources of error for the reference - display system and 
determine the values of certain parameters to produce the 
least amount of error. 

Although the probability of a worst case error 

situation occuring is small, an analysis of this type was 

performed instead of a statistical analysis to simplify the 

calculations ana get a general idea of the feasibility of 

the reference - display system. The design specification is 

that the total system error must be less than or equal to 

1.0 millimeters. This analysis indicates that there is only 

one adjustable parameter, the microphone - spark gap 

distance, which is actually not that adjustable given the 

operating room constraints. 

The computation error (e~ ) is a function of the 

comp 

digitizer error (e^) f focusing error (T f ), and, most 
importantly, the algorithms used to reference the system. 
As mentioned in section 7.4, this error can propagate 
through the various registration steps and calculations. 
Therefore, to determine this error, the reference algorithms 
must be written and evaluated then adjusted appropriately. 
This work is described in Chapter 9. 



Finally, the error due to interpolation (e^) will be 
estimated. and adjusted based on the CT scan reconstruction 
procedure not covered in this thesis. 

Unfortunately, since the CT scan resolution, focusing 
error and fiducial error cannot be improved by adjusting 
available parameters and the magnitude of T f is at least one 
millimeter, the total worst case error cannot be less than 
or equal to one millimeter. 

Since these errors can act in different directions and 
this is a worst case analysis, the concept of develocir.3 2 
practical system has not been disproven. 



Chapter 3 - Spark Gap ,-lul tiplexer D esign 

Since three spark gaps are needed to determine three 
points on the microscope and, therefore, its position, in 
space, and the digitizer we purchased has only one spark 
line, a device is needed to switch the single input spark 
line to the three spark gaps. Science Accessories sells a 
multiplexer that can switch up to eight spark gaps for 
33,500, but since only three sparks gaps need to be switched 
and the switching rate can be slow, it was less expensive to 
design our own three spark gap multiplexer. This chapter 
describes the specifications, constraints and problems 
associated with the building and design of the multiplexer. 

3.1 Design Specifications/Constraints 

The design specifications/constraints were to build a 
switching circuit that could direct the incoming spark gap 
signal from the digitizer to one of three user selected 
spark gaps. More specifically:" 1) The switching circuit 
must handle the 350 volt spark signal and not allow it to" 
arc within the switch. 2) The circuit must be properly 
shielded to prevent inductive interference with the IBM PC 
and the digitizer. 3) The control circuit must allow both 
manual and computer control of the selected spark gaps. 4) 
The final cost must be less than $3,500, including labor for 
development . 
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6.2 Circuit Description 

The circuit can be divided into five parts: control 
circuitry, switching circuitry, display circuitry, power 
supply circuitry and shielding considerations. Refer to 
Appendix S for schematics and circuit layouts. 

3.2.1 Control Circuit 

In order to make the multiplexer computer or manually 

mounted on the front panel. In manual mode, the user can 
select the appropriate spark gap by two toggle switch*:; 
which turn on and off two 5 volt lines (TTL logic-) . In 
computer mode, the input to the control circuitry is tne : i 
volt TTL parallel port output from the TdA sent to the 
circuit by a back panel mounted D3-25 connector. The switch 
select inputs are processed by a f4lj9 2-4 channel decoder 
which outputs three high signals and one low signal for the 
unique two bit input signal. Since three spark gaps are 
used only three of the four decoder outputs are connected. 
This signal is then inverted by a f405 open collector quad 
inverter so that the switch selected is high and the others 
are low C 47 3 . This signal then goes to the power 
transistors . 
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A back panel mounted jack, switched on and off by a 
?RAk 1 AOS reed relay and controlled by the ISM parallel 
port, was added for accessories. 

3.2.2 Switching Circuit 

The main part of the switching circuit are the single 
pole-single throw mercury wetted reed relays R1 , R2, and R3. 
A 5 volt, 91 milliamp signal is needed to close the normally 
open switch. Mercury wetted relays were purchased that have 
a dielectric breakdown greater than 1000 volts dc so the 
signal will not arc in either of the two open switches. 
Mercury wetted reed relays are position sensitive so they 
had to be mounted vertically. Diodes were placed across 
each relay control input to prevent possible reverse current 
surges from damaging the relay and 0.1 microfarad capacitors 
were placed from 5 volts to ground to reduce power line 
noise. Since standard TTL control logic does not carry 
enough current to close the relays, TIP29A power transistor 
stages were , added to draw enough current through the relay 
when the base signal goes low, 

8.2.3 Display Circuit 

The display circuitry consists of three 2NJ906 bipolar 
transistors that sense the power transistor collector state 
and open and close their collector-emitter junction when the 
base voltage is +5 volts and ground, respectively. rfhen the 
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circuit closes, the current passes through the appropriate 
L£D indicating the multiplexer status. 

8.2.4 Power Supply Circuit 

The power supply circuit consists of a modular Polytron 
120 volt ac to 5 volt dc, 500 rnilliamp power supply. A 100 
microfarad electrolytic capacitor and a 0.1 microfarad 
ceramic capacitor were placed between 5 volts and ground fc to 
reduce power supply noise. The power supply is protected by 
a slow blow fuse and has a back panel mounted on-off switch. 

3.2.5 Shielding and Grounding 

The main problem with designing a switching network to 
handle a 350 volt spike signal is that the sharp voltage 
transition induces current in nearby wires causing the 
digitizer to make false readings or other spark gaps to 
fire. To eliminate this problem, all connections going into 
and out of the multiplexer were made with shielded coaxial 
cable whose integrity was internally maintained by coax 
connections and copper ground plane switch circuitry. The 
switch board was separate from the control and power circuit 
boards for further isolation. 

To prevent ground loops internally and between the 
digitizer and multiplexer all circuit components were 
soldered point-to-point with a common ground connection. 
See Figures 3. I, 8.2, and 3.3. 



Chapter j - Software and Computation irro r 

This chapter describes the software written in BASIC 
for the 13.1 ?C XT with the 3037 coprocessor (sections 9.1 - 
9.4) and includes a quantitative error analysis for each 
step in the reference procedure (section 9.5 - 9.6). 

The software is divided into two files. File" FPL 
perforins the steps in section 5. U.I and establishes the 
relationship between the relative position of the spark gaps 
and the focal plane by calculating the oblique spark gap . 
coordinates of the focal (5.4.1.1), normal (6.4.1,2) and 
orientation (5.4.1.3) points. This program needs to be run 
once for each spark gap bracket (sitting or supine) and only 
occasionally thereafter to reestablish the spark gap focal 
plane relationsnip . 

File REG performs the steps in sections 6.4.2 and 
6.4.J, and is run for each operation in the operating room 
once the CT coordinates of the fiducials have been 
determined by the Cf scanner/ computer . The first step is 
the registration procedure (6.4.2) to calculate the CT 
coordinates of the microphones by focusing on each fiducial. 
The second step (5.4.3) calculates the equation of the focal 
plane, the coordinates of the focal point and the three 
direction cosines of the Y.^ microscope axis 'in CT 
coordinates for each microscope position. The software and 
documentation are in Appendix C. 
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9.1 Digitizer Control 

The subroutine SRINP, used by both files FPL and REG, 
was written to control the digitizer and the spark gap 
multiplexer and to store the distances in an array. The 
first part writes a binary code (30, 01, or 10) to the 
parallel port (address 33C hex) on the IBM monochrome 
adapter card which sends a signal to the spark g.ap 
multiplexer and sets the appropriate relay. Then the RS-232 
serial port of the I3rt is opened to receive the slant range 
data from the digitizer. This data is sent as a 26 
character string as shown below. 

11 11 1S22222S33333S444U4CRLF 

1, 2, 3 and 4 represent the five digit metric slant range 
distances for each microphone, S represents .a space 
character, and CRLF is a carriage return and line feed. The 
string is searched until an ASCII line feed character is 
found marking the end of a set of four slant ranges. The 
next 26 character string is from the digitizer's "last 
sample" buffer and is ignored. The data is then received 26 
characters at a time until the desired number of samples has 
been obtained. Then decimal points are inserted to read the 
slant ranges in centimeters (XXX. XX) and then stored in 
array SLANT. While the data is received, the IBM 
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communications input buffer is .nonitored, interrupting and 
resuming the data flow depending on available buffer space. 

Since the operating room equipment and patient are 
generally in two different setups (depending on the 
craniotomy site) blocking a different microphone each ti.ne, 
this subroutine throws out the slant range data from the 
blocked microphone, which is determined when the user 
identifies the procedure (left temporal or right frontal 
craniotomy) . The I3M serial port is then closed by writing 
a binary c=Ja (1010) to the port's address (3F: hex) 
disabling digitizer communication and turning off the spark 
gaps. The means and standard deviations of the slant ranges- 
are then calculated. If the data from any microphone has 
more than five bad data points (distances > 250 centimeters, 
indicating the signal was blocked or never received) or the 
standard deviations are found to be greater than a 
predetermined limit (SDLIMID set in the software (presently 
set at 0.05 cm.), the subroutine will prompt the user to try 
again . 

After the mean slant ranges are stored in array SRMEAM , 
the spark gap multiplexer switches to the next spark gap and 
repeats the process until all three spark gaps have been 
fired. 
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9.2 Spark Gap - focal Plane Relationship 

This section discusses the subroutine FPLAIJS (in file 
FPL), which executes the procedure described in Section 
6.4.1 and calculates the oblique spark gap coordinates of 
the focal, normal and orientation points for the sitting or 
supine spark gap brackets. File FPL requests the type of' 
spark gap bracket, sitting or supine, and sets FFLAG such 
that the oblique coordinates are stored in the proper file". 

The user is then instructed to anchor the fourth spark 
gap under the microphones , being careful not to block the 
spark gap - microphone line of sight with the microscope., 
SFLAG is then set to indicate to subroutine SRI.JP that only 
one spark gap should be fired. SRINP is called and the 
square of the slant range distances from the fourth spark 
gap are stored in array DIST. The user is then asked to 
focus on the tip of the fourth .spark gap and call SRINP, 
firing the three spark gaps on the microscope. The oblique 
spark gap coordinates of the three microphones are then 
calculated by following Equations 6.3, 6.4, 6.S, 6.9, and 
6.12. Once these are determined and stored in array FX, 
subroutine NEWTOM is called with NFLAG equal to one, 
indicating that the cosine of the angle w (Equation 6.1) 
will be represented by CW. 

NEWTOM, with tNFLAG equal to one, specifically solves 
three nonlinear equations determined by the oblique distance 
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formula (Equation 5.13) [10j. Preset negative initial 
values are used since the z coordinate of the focal point is 
negative and positive initial values could cause the 
solutions to converge to a geometrically valid mirror result 
above the x-y plane. N'ewton 1 s method solves equations in 
the form f^x) = 0 where F^x) is the vector of i equations. 
The solution vector can be found by calculating: 

x.<*> = x^- 1 ' - J .u i (1 <- 1 >rV 1 u 1 ll <-") (9.i)' ■ 

where k is the iteration number (X^ 1 ^ are the initial 
values), and J^X^^" 1 ^) is the Jacobian matrix of partial 
derivatives of each equation i with respect to the unknowns 
X i evaluated at X.^" 1 * [10]. Since i = 1 to 3 , the 
Jacobian matrix is 3x3: 
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NEKTON calculates the Jacobian matrix and calls subroutine 
INVERSE which inverts the Jacobian matrix. Equation 9.1 is 
then solved for each iteration, k, until either 30 
iterations are reached or until; 



(x (k). x (k-1) )/x (k) < r0L 

1 1 1 - 

where fOL is a set tolerance value of 0.001 centimeter. The 
user is notified if the solutions have not converged in less 
than *j0 iterations. 

Subroutine INVERSE uses the Crout algorithm (a 
variation of the standard Gaussian Elimination technique 
only with maximal pivoting strategies) to invert the 3x3 
Jacobian matrix [10]. The elements of the inverted matrix- 
are solutions co the equations: 

vr 1 = 1 

where I is the 3x3 identity matrix. 

The solutions to Equation 9.1 can be calculated 
directly by Gaussian Elimination by letting Z = J^ 1 ^ and > 
therefore, J-^Z = F i and then solving for Z. 

At this point, the solutions to Equation 6.13, the 
oblique spark gap coordinates of the focal point, are stored 
in array FOCAL and section 6. 4,1.1 is complete. 

After solving for the focal point, FPL AN E instructs the 
user to focus the microscope on Points 1, 2 and 3 of a fixed 
grid coordinate system with the crosshairs of a small 
protractor device at the origin. This is to establish a 
fixed coordinate system such that the fixed coordinates of 
the microphones can be determined. See Figure 9.1. 
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POINT 1 




Figure 9.1 - Grid Coordinate System 

rfhen focusing on Point 3, the user is asked to align 
the microscope crosshairs with the protractor crosshairs and 
focus as normal to the protractor as possible. Then the 
user is instructed to move the microscope away from the 
protractor by turning only the focusing knob of the 
microscope. The slant range distances for each of these 
points is stored in array SDATA. 

Subroutine DISTANCE converts the slant range data for 
Points 1, 2 and 3, in SDATA, into the oblique spark gap 
coordinates of the three microphones (Equations 6.8, 6.9 and 
6.12) and calculates the distances between the microphones 
and focal points (Points 1, 2 and J) by the oblique distance 
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formula (Equation 5.14). The fixed coordinates of the three 
microphones are then determined by solving three sets of 
three simultaneous nonlinear equations (Equation 6.15). 
Then MFLAG is set to zero to let cos(w)=0 and allow 
subroutine NEWTON to solve the three sets (one for each 
microphone) of three nonlinear equations (Equation 6.15). 

At this point, we know the fixed coordinates of each 
microphone relative to the grid coordinate system. These 
coordinates are stored in array MIKE. 

The next step is to calculate the fixed coordinates of 
the normal point (5.4.1.2) based on the slant range data in 
SDATA. Subroutine DISTANCE is called again and the 
distances between the microphones and normal point are 
calculated and stored in array DIST. NEWTON is called again 
with NFLAG equal to zero and solves the three nonlinear 
equations o-f Equation 6.16 for the fixed coordinates of the 
normal point. 

The oblique spark gap coordinates of the normal point 
are then calculated by first determining the normal point - 
microphone distances by applying the standard Cartesian 
distance formula of Equation 6.15. NEWTON is then called 
again (NFLAG=1 ) to solve Equation 6.17 for the oblique spark 
gap coordinates of the normal point. These oblique 
coordinates are then added to array FOCAL. 

The next step is to follow the procedure of section 
6.4.1.3 and calculate the oblique spark gap coordinates of 
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the orientation point. The first step is to calculate the 
coefficients of the focal plane at Point 3 of the grid 
coordinate system by solving Ration 5.21. Once the focal 
plane is known, Point 1 is projected onto the focal plane 
and the fixed coordinates of the orientation point is 
determined by Equations 6.23, 6.24, 6.25 and 5.26. The 
distances between the orientation point and the microphones 
are then determined by applying Equation 5.15 and the 
oblique spark gap coordinates of the orientation point - 
calculated by calling subroutine NEvfTON WFLAG.1 ) to solve 
Equation 6.27. The oblique coordinates are then stored in 
array FOCAL . 

The oblique spark gap coordinates of the focal, normal 

vv, a n t-rsnsferred from FOCAL and 
and orientation points are then transterrea 

stored in files SIT or SUP depending on the spark gap 
position and surgical procedure. 

Registrati on Procedure 
File REG is the driver for the registration procedure 
and focal plane determination. Once the patienf s head has 
been fixed in the head clamp, the surgeon is asked to 
•indicate which type of procedure will be performed, left 
temporal craniotomy or right frontal craniotomy (see Section 
2 „ This will indicate which of the four microphones 
(slant range values) will not be used. Should another, less 
common, operational setup be required, the microphone that 



1 1 1 



will not be used can be indicated in the software prior to 
the procedure by changing M.-iOT . The surgeon is also asked 
the position of the patient for the procedure (sitting or 
supine) , which will determine the appropriate spark gap 
bracket (see Chapter 10) and, therefore, relative positions 
of the spark gaps and microscope. Depending on the 
position, file SIT or SJP will be read and the oblique spark 
gap coordinates of tne focal, normal and orientation points 
stored in array FOCAL. 

Subroutine REGISTER is called and the surgeon is then 
instructed to focus on each of the three fiducials, pausing 
between each one for SRINP to fire the spark gaps and 
determine the slant range distances. REGISTER then asks the 
surgeon to enter the CT coordinates of each fiducial. 
Subroutine DISTANCE determines the distances from the 
microphones to the fiducials and NEWTOrJ solves the three 
sets of simultaneous nonlinear equations for the coordinates 
of the microphones in CT coordinates and stores the data in 
array MIKE (Equation 6.15). 

The CT coordinates of the array MIKE are then written 
to file MIKE . Therefore, if the I3M fails, the registration 
procedure need not be repeated. The user is asked at the 
beginning of the procedure whether or not the computer 
failed after the registration procedure. The microphone 
coordinates are then recalled, if necessary, and the 
registration procedure is skipped. 
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9.M Focal Plane Determinat ion 

With the microphone coordinates determined in CT 
coordinates, the surgeon is then instructed to remove the 
fiducials, if necessary, and continue with preparing the 
operative site and perform the craniotomy. If the 
craniotomy site has not been determined, the reference - 
display system can be used to reconstruct and display images 
that might be useful in planning where to "open". REG then 
calls subroutine SLICE, which instructs the surgeon to 
indicate wnen the microscope is focused so SRIWP will be 
called to fire the spark gaps. 

The CT coordinates of the focal, normal and orientation 
points are then determined by subroutines DISTANCE and 
WEWTOrl, which specifically solves Equation 6.23 as explained 
in Section 6.4.2.2. The equation of the focal plane is 
determined by Equation 6.21* and the direction cosines 
calculated by Equations 6.30, 6.31 and 6.32. The focal 
plane coefficients and direction cosines are displayed oh 
the screen along with the coordinates of the focal point. 
Once communication between the ISM PC and the Treatment 
Planning Computer has been established, this data will be 
transmitted directly to the Data General Eclipse for the 
image reconstruction and display. 
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9.5 Error Analysis 

The purpose of tnis section is to quantitatively 

determine the error due to the computat ions (e ) 

comp 

described in this chapter and outlined in Chapter 7. 

9.5,1 FPL Error 

The computation error due to establishing' the spark gap 
- microscope relationship with the program FPL ( e C0[n p_crpL ) 

Se 3 f''np^ i "n f rV»a H i 71 f i' 7 o r error o -» ^ ^ a r -\.>n« ^ ^ • 
-j - 1 ^ _ - - 0 _ „ ~ ; 1 w ' I -•*'- j 

error, eT f . Vhis error was evaluated by calculating the 

oblique spark gap coordinates of the focal point, adding the 

magnitude of e^ to the slant ranges and the magnitude of e^. 

to focal point coordinates and then observing the magnitude 

of the resulting error as a function of |e*^|. See Figure 9.2. 

The solid line is the error when k = +0.1} (worst case bias 

percentage) and the dashed line is with k = +0.0341, the 

expected worst case based on operating room temperature 

variations. It can be seen that with a focusing error of 1 

millimeter, e_ rDT = 1.8 millimeters with a 0.5 
COiflp-r rL 

millimeter change in e coinp- cpL P er millimeter change in 
e* f . This value is a combination of worst case errors for e* f 
and e^, but for typical slant range distances stored in file 
TSLAdT . 

There is another source of error that should be 
uientioned. The oblique coordinates of the focal, normal and 
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orientation points will be determined prior to the procedure 
and their • accuracy is dependent on the precision of the 
spark gap holder, which will be removed and then replaced 
when the microscope is draped. This will introduce another 
error, e g ^ 0 ^j er i which was experimentally determined (see 
Appendix F) to be 0,9 millimeters and in a worst case 
analysis can be added to the propagation of the focusing 
error. This error might be reduced by redesigning the spark 
gap holder - microscope mount, which may eventually be 
changed should the procedure in the operating room introduce 
more physical constraints. 
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9.5.2 REGISTER Error 

The computation error due to the propagation of 

e , T . u and e*, through the registration 

comp-FPL' f s-holder d 

procedure can also be quantitatively determined. Again, the 
average error in the determination of the microphone 
coordinates is based on typical slant range data (stored in 
file TSLANT) and shown in Figure 9.3. This is also shown as 
a function of the focusing error as this error is amplified 
oy focusing on each fiaucial. Again, the solia line 
represents k = +0.1* and the dashed, k = +0.034*. e comp _RjTG 
at |e* f | = 1 millimeter is approximately 5.5 millimeters where, 
in general, a 0.5 millimeter change in le^l corresponds to a 
0.5 millimeter change in e comp _ REG . 
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9.5.j SLICE £rror 

e coinp-SLICE re P resents 'the magnitude of errors due to 
e coinp-R£G' ®d and ®f as they have P ro P a Ssted through the 
calculations of the focal point for a typical reconstructed 
CT image. It is shown in Figure 9.4 as a function of 
focusing error for k = +0.15 (solid line) and k = +0.0j*ij* 

(dashed line). With |e*J at 1 millimeter, e m CT T „ = 3.4 

i comp-SLICc 

millimeters. 

It id interesting to note that the magnitude oi cae 

final error, e^ m ^ c T T rr , is less than the sum of e _ „ r 
comp-SLxCE 1 C0Jip-r ?L 

and e^ m ^ This is because each of these errors is a 

comp-HEu 

vector and acts in specific directions, partially cancelling 
others out. 
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9.6 Error Conclusions 



This error analysis is a combination of a typical and 
worst case situation with typical slant range values and 
worst case error parameters. The total theoretical worst 
case error for the entire reference - display system is the 
sum of the propagated error, e CO ; flp _3Lic£ » the fiducial 
error, eT^^i and the interpolation error, ef\: 

£ T = e como-3LICE + ' T f id 1 + |T i"' 

^omp-SLICE' for a focusin S error of 1 millimeter, is j.4 
millimeters. Therefore, 

E T = 3.4 + ir f . d l + le^I, 

which is unacceptable for neurosurgical procedures. From 
this analysis, the error due to focusing appears to be the 
most sensitive parameter and unfortunately cannot be easily 
adjusted . 

At this point the total worst case error is greater 
than the specification of one millimeter. This means that a 
more statistical approach should be considered to determine 
the probability of more common situations. However, should 
this large error prove to be the case experimentally, new 
algorithms will have to be developed to try to reduce the 
propagation of fixed errors such as focusing. Attaching 
three additional spark gaps to the three fiducials might 
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reduce the focus ing • error but since the microphone placement 
is restricted to the operating room ceiling between the air 
vent and light tracks (see Figure 2.4) the digitizer line of 
sight constraint might geometrically prevent proper sound 
transmission. The experimental results will indicate the 
next step. 



Chapter 1 0 - Mechanical Design 

This chapter describes the hardware that required 
mechanical design and machining. Spark gap holders were 
designed to mount the spark gaps on the microscope for two 
different surgical positions. Microphone mounts were also 
designed to allow quick and simple placement of the 
inicrophone array on the operating room ceiling. 

10.1 Spark Gap Holder 

The specifications for the design of the spark gap 
holder are: 1) the spark gap separation can be no greater 
than jO centimeters to prevent interference with the patient 
drape and scrub nurse's procedure, 2) the weight of the 
holder and spark gaps must be less than 2 kilograms, 3) the 
spark gap holder must mount on the microscope such as not to 
interfere with the surgical procedure and manipulation of 
the microscope, 4) there must be a clear line of sight 
between all three spark gaps and -at least three microphones, 
5) since the spark gaps and holder will be in the vicinity 
of the operative, field they must withstand aseptic 
sterilization, and 6) the spark gaps must be easily and 
precisely connected to the microscope after the drape bag is 
in place. 

The spark gaps are located at the vertices of an 
equilateral triangle as determined in section 5.4. The 
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weight specification was met by building the holder cue of 
plexiglas. The spark gap holder attaches easily to the 
microscope via a male counterpart to the female dovetail 
notch on the microscope, which was designed for the 
attachment of a surgical laser. fhis meets specification 3. 
See figure 10.1. In order to meet the line of sight 
criteria, two different holder arms were designed for the 
cases where the patient is supine or in a sitting position. 
The spark gap holder must undergo ethylene oxide gas 
sterilization because the spark gap tips are made of 
stainless steel and would oxidize with aqueous 
sterilization.' The precise attachment of the spark gap 




Firure iVx - Soark Gap - Microscope Ilount 



121 



holder to the microscope was done by machining an aluminum 
"tongue" that attaches to the laser mount and a female slot 
with two set screws that can precisely relocate the spark 
gaps to the same position. The microscope drape bag is 
stretched over the metal tongue until the tongue protrudes, 
temporarily breaking the sterile barrier. The spark gap 
holder arm is then slid over the tongue and the set screws 
are tightened, anchoring the spark gap holder. See Figures 
10.2, 10.3 and 10.4. 

10.2 Microphone Mounts 

The only major specification for the design of the 
microphone mounts is that the microphones and preamps be 
ceiling mountable in the operating room within one minute 
each so the operative procedure is not delayed. The mounts 
must not be damaged by general .operating room preparation 
and the microphone "array must be removable, allowing their 
use in other operating rooms. Since the exact orientation 
of the four microphones is not critical, they do not need to 
be mounted on a rigid framework. 

The microphones were attached to a small plexiglas 
holder that was mounted on the preamp boxes described in 
section 5.4. The preamp boxes were attached to an aluminum 
male dovetail section whose female counterpart was mounted 
on the operating room ceiling between the light tracks and 



Fiijure 10.2 - Supine Sparl: Gap Holder 




Finure 10.3 - Sitting Spark Gap HolJer 



123 




Figure 10. 1- - .Microscope Mounted Spark Gap Holder 



the air vents. A set screw was made to allow the secure 
positioning of the preamp boxes. The microphones and preamp 
boxes can be quickly mounted by sliding the dovetail pieces 
together and turning the set screws. See Figures 10.5 and 
10.6. The mechanical drawings for the spark gap holder and 
microphone mounts are in Appendix D. 
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Figure Ij.G - Microphor.e Mount - Viaw 2 



Cha pter 1 1 - System Evaluation 

This chapter describes the reference system evaluation 
including digitizer performance (11.1), phantom design 
(11.2), error analysis (11. J), experimental results and 
conclusions (11.4). 

11.1 Digitizer Evaluation 

The digitizer was tested for accuracy, and precision 
under the following test criteria: 1) single spark gap - 
microphone separation, 2) spark gap at 0° to the microphone 
(see figure 11.1), 3) digitizing rate at 30 points per 
second, 4) sampling times of three seconds for 50 samples, 
5.5 seconds for 100 samples, and 3 seconds for 150 samples, 
5) measured spark gap microphone distances were 53 ± 0.1 cm, 
100 * 0.2 cm, 150 + 0.J cm, and 200 + 0.H cm. The test 
statistics are in Appendix E . The conclusions from this 
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Fi pure 11,1 
Spark Gap - Microphone Orientation 
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test were that the digitizer meets its specifications for 
precision, + 0.01 cm, but it was difficult to determine the 
results of the accuracy test because we do not have a 
measuring device that is more accurate than the digitizer 
itself. The slant ranges do appear to be within + 0.1£ of 
the measured distance. 

11.2 Test Design 

In order to determine the accuracy of the reference 
system, a testing technique had to be developed. Since a 
goal of this Thesis is to prove that the calculated focal 
plane is within a specified error (initially 1.0 mm) of the 
correct focal plane, an evaluation introducing the fewest 
sources of error is preferred. The error that will result 
will be easier to trace back to its source. 

A block of plexiglas was milled square approximately 
the size of a human head (13x10x10 cm.) and small points 
(0.8 mm diameter) were drilled into it with known 
separations. Three of these points represent fiducials and' 
the other five are test points whose coordinates are known 
with respect to the coordinate system defined by the 
fiducials. See Figure 11.2. The coordinates of the 
fiducial points and test points are given in Table 11.1 
below. 
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Table 11.1 - Test Point Coordinates^™) 



Fiducials 



Point 



X 



Y 



Z 



2 
3 



0.0 

-1 .6n 

1 .016 



0.0 
7.62 
15.367 



0.0 
4.854 
1 .75 



Test Points 



2 

3 
4 

5 



-4.719 
2.54 

5.949 
6.096 
1.016 



J .62 
0.0 
7 .62 
3 .24 7 
15.367 



4.354 

0.0 

0.226 



0.79 
0.0 



This test phantom will eliminate the errors related ts 
CT scanning and display/reconstruction. The error 
determined at the focal point will be the distance between 
the real and calculated CT coordinates and the only sources 
of error will be due to slant range, focusing and" 
computation error. A subroutine TEST was written to call 
the required subroutines and evaluate the focal point error 
and angular focal plane orientation error. The error at the 
focal point can be easily determined since the coordinates 
of the test points are known. The error in the focal plane 
calculations can be determined qualitatively by refocusing 
on the grid coordinate system (see Figure 9.1) normal to the 
x-y plane along the z axis and determining if the focal 
plane coefficients represent the plane z=0. The x and y 
direction cosines will indicate angles that add up to 90° if 
the orientation vector is in the first quadrant of the grid 
coordinate system. 
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Figure 11.2 - Test Phantom 

11 .3 Test Error Analysis 

This section describes the total worst case focal point 
error that could be expected by using the test procedure 
outlined above. 

The worse case error for the test phantom is given by: 

E test = € comp-SLICE + |e fid-test'' 

where e -„ is j.4 millimeters. The fiducial error 

comp-^LICt 

for the test phantom, |e f id _ test li is different than that for 
the general system. Each test point and fiducial was marked 
within 0.12f millimeters (J. 005 inches), so ' e fi<j-test' is 
0.127 millimeters . 

£. aeV = 3-4 + 0.1J = J.5j millimeters 



The worst case focal point error is 3-5 millimeters. 

The worst case test focal plane error is hard to 
estimate as the focal plane is defined as four coefficients. 
It will be analyzed as described in Section 11.2, above. 

11. U -Experimental Results 

The raw data for the experiments is in Appendix F with 
a brief description. There are three tests that were 
conducted to characterize the reference system, error: focal 
point accuracy and precision, resting focal point precision 
and focal plane accuracy. 

The focal point accuracy and precision was determined 
by first focusing on each of the three fiducials and 
calculating the coordinates of the microphones in the test 
phantom' coordinate system. fhe microscope was then focused 
on each of the various test points on the test phantom 
multiple times at different angles and the magnitude of the 
error at the focal point was determined. fhe average error 
was 1.1 millimeters, less than the estimated worst case 
error (11.3), with an average standard deviation of 0.52 
millimeters. If the errors due to the fiducials (e* fid ) and 
interpolation (e^) also prove to be less than expected, the 
reference system will have wide range applicability in the 
operating room. The focal point precision or repeatablity 
is determined by the standard deviation. If we assume that 
the magnitude of the focal point error follows a standard 



normal distribution, the standard deviation indicates that 
53* of the focal point values will be within 0.53 and 1.62 
millimeters (one standard deviation), and 95% will be within 
0.06 and 2.14 millimeters (two standard deviations). 

The resting focal point precision was determined by 
focusing on the same point and sampling fifteen times to 
determine the variation in the coordinates due to fixed 
errors such as air motions, etc. The standard deviations in 
centimeters of each coordinate were determined and found to 
be 3.351 Cx), .0.066 (y) and 0.055 (z) with values failing 
evenly on either side of the mean. This indicates that the 
variations in error values at the focal point could be 
partially due to a fixed bias error which cannot be removed 
from the system by changing available parameters. 

The focal plane is defined in the CT coordinate system 
by the four coefficients of the equation of the plane: Ax + 
3y + Cz = D . As mentioned earlier in this chapter, by 
focusing the microscope along thez axis of the CT 
coordinate system (test phantom) at the origin ,. normal to 
the x-y CT plane, the resulting focal plane equation should 
be approximately z=0 (A=0, B=0 , C = 1., °=0). This experiment 
was conducted by refocusing along the z axis of the grid 
coordinate system at the origin (Point 3) fifteen times and 
evaluating the equation of the focal plane. The average 
values for the plane coefficients were A=-0.0"j2, o=0.106, 
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C=l, and 0 = 0. 005. This, indicates that the focal plane has 
been determined without any gross error. 

A more rigorous test of focal plane accuracy would be 
to determine the C'f (fixed) coordinates of three corners of 
a large flat surface by focusing on each potnt. The 
equation of this plane could then be calculated. Then, by 
focusing at the center of the surface, the focal plane 
coefficients could be determined and compared to the. surface 
plane coefficients. This error could then be translated 
into focal plane error at the edge of the microscope field. 

The next step is to evaluate the accuracy of the 
orientation point coordinates, or specifically the direction 
cosines indicating the orientation of the focal plane about 
the optical axis. 3y focusing the microscope along the z 
axis at the x-y plane again and keeping the y axis of the 
microscope coordinate system crosshair in the ocular) in 
the first quadrant of the 'grid coordinate system, the sum of 
the direction cosines should be 1. This corresponds to an 
angular sum of 90° from the x and y grid axes to the 
orientation vector. The average angular sum was 90.3° and, 
again, only indicates that there have not been any major 
errors . 



Chapter 12 - C inclusion 

12.1 Reference - Display System Conclusions 

The reference - display system, under simulated 
operating rootn conditions, was found to be accurate to 
approximately 1 millimeter at the focal point and the 
equation of the focal plane was found to be' qualitatively as 
accurate. The initial specification was to have the entire 
operational system error within 1 millimeter. The errors in 
other parts of the project, CT scanning, fiducials, etc.," 
will certainly increase the resulting reference - -1is?lv 
system error, however, the magnitude of this total error and 
the conceptual success and resulting project applicability 
cannot be determined until the system is actually used with 
a patient in the operating room. 

Up to this point the project has been quite successful 
and hopefully the CT image reconstruction phase will 
completely prove the concept of integrating CT scanning and 
the operating microscope. 

12.2 Future Design 

As mentioned in Chapter 4, redesigning the microscope 
stand might be a better way to determine the relative focal 
plane - CT scan positions. Although it has not been 
determined as being cost effective, this design could take 
into account all the surgeon's needs and constraints while 
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providing a microprocessor-based system to determine, via 
position transducers the focal plane with respect to a fixed 
CT coordinate system. Since the present system is most 
sensitive to focusing error, the microscope optics could be 
redesigned to minimize the depth of field. 

The display system could be improved by incorporating a 
high resolution color CRT with thinner beam splitter and 
reducing the extension of the CRT from the optical axis. - 

The optimal computer system design would directly 
receive the CI scan d2t2 frc: the scanning ccuputar zr.i 
process it right in the operating room. The reconstruction 
and display process could be done with the hardware of the 
present generation CT scanners, reducing the calculation 
time . 

A design of this magnitude and expense could be a 
future goal. 

12.j future Work 

The future work on the current project will primarily 
involve the design of CT scan reconstruction software. 
Although contours will be developed initially, the 
algorithms should be adaptable to grey scale information. 
An anatomic brain atlas could also be developed to indicate 
important neurological structures. 

An important feature that should be added to the 
reference - display system is feedback. If a" 3-0 trackball 
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controlled the position of the CRT. image, the surgeon could 
compensate for obvious reconstructed slice reference errors 
by correcting the image placement according to visible 
anatomic landmarks. The results of focusing error and other 
disturbances, which are significant and more or less fixed, 
could then be eliminated in successive slices by storing the 
offset values from the trackball. 
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Appendix A - CR: Schematic Changes 



This appendix indicates the electronic changes to the 
J.tf.C. model vf-1'jOO viewfinder (miniature CRD used to 
display the reconstructed CT image. 

The horizontal hold potentiometer, R23 was changed from 
4.f k to 10 k, and resistor R29 was changed from 12 k to 10 
k. These changes will now allow more adjustment of the 
horizontal hold and provide a stable image for the Treatment 
Planning Computer display. 
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A ppendix 2 - .-iu ltiple xer 3cne:astic and L ayout 

This appendix contains the schematics and layouts of 
the spark gap multiplexer described in Chapter 3. 
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This appendix contains the software described 
Chapter 9. 
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20 9 APPENDIX C 
"30 " 

40 * THESIS TITLE: REFERENCE - DISPLAY SYSTEM FOR THE INTEGRATION OF CT 
50 » SCANNING AND THE OPERATING MICROSCOPE 
oO ' 

70 9 FILENAME: FPL 

30 9 

*0 ' PROGRAMMER: JOHN F. HATCH 
XOO " DATE: AUGUST. 1984 

110 * 

120 * NOTE i REFER TO CHAPTER 6 IN THESIS FOR EQUATIONS. 
130 ' AND REFERENCE 10 FOR ALGORITHMS. 

140 9 

130 » THIS FILE DETERMINES THE COORDINATES OF THE FOCAL, 

160 ' NORMAL AND ORIENTATION POINTS IN OBLIQUE SPARK GAP COORDINATES. 

170 * FOR BOTH THE SITTING AND SUPINE SPARK GAP BRACKETS. AND STORES 

180 ' THEM IN FILES SIT AND SUP. 

190 ' 

200 OPTION BASE 1 ' SET LOWER ARRAY INDEX TO I FOR CONVENIENCE 

210 ' 

220 ALLOCATE SPACE FOR ARRAYS 

230 9 NOTE: ARRAY NAMES FOLLOWED BY A SIGN ARE DOUBLE PRECISION 

230 DIM J#<3.3) ' JACOBIAN MATRIX COEFFICIENTS 

260 DIM I#'.3,3) ' JACOBIAN MATRIX INVERSE 

270 DIM 5RMEAN(3,3) > MEAN SLANT RANGES 
280 DIM SDATA (4,3,3) 9 SLANT RANGE DATA 

2*0 DIM FOCAL (3.3) 9 OBLIQUE SPARK GAP COORDINATES OF THE FOCAL, NORMAL 

-.00 ' AND ORIENTATION POINTS 

310 DIM NS0LC3,3) 9 COORDINATE SOLUTIONS TO NEWTON 

320 DIM FX (3.3) FIXED GRID COORDINATE SYSTEM 

330 DIM MV#(3> 9 ITERATION VARIABLE <N-1) 

340 DIM M* (100,3) f ITERATION VARIABLES 

330 DIM SL<3> * SLANT RANGE DATA 

3*0 DIM T(3) * ITERATION ERRORS 

370 DIM DIST(3.3> * MICROPHONE-FOCAL POINT DISTANCES 

330 DIM MIKE(3,3) ' MICROPHONE COORDINATES IN FIXED GRID COORDINATE SYSTEM 
3^0 DIM NP0INT(3) 9 NORMAL POINT COORDINATES IN FIXED GRID COORDINATE SYS. 
400 DIM 0P0INT(3) ' ORIENTATION POINT COORDINATES IN FIXED GRID SYSTEM 
410 DIM GRID<3.3) ' COORDINATES OF- MICROPHONE PLANE GRID 
420 DIM SLANT (150.-3) 9 SLANT RANGES 

^30 DIM Ft(4) * ITERATIVE SOLUTION TO NEWTON 

440 " 

430 ' MEASURED DISTANCES BETWEEN SPARK GAPS IN CM. 
460 S12-30.041 ' SPARK GAP 1 - SPARK GAP 2 

470 S23-29.993 ' SPARK GAP 2 - SPARK GAP 3 

480 S3 1-29. 637 9 SPARK GAP 3 - SPARK GAP 1 

490 9 

300 CW- CS12~2*S31~2-S23~2) / (2*S12*S31 > 9 COSINE OF ANGLE W, S2-S1-S3 EQN. 6.1 
310 S2W-1-CW1CW 9 SINE SQUARED OF ANGLE W. S2-S1-S3 

320 9 

330 ' IGNORE SLANT RANGE DATA FROM MICROPHONE MNOT 
340 MNOT-4 
330 ' 

370 PRINT-WHICH SPARK GAP BRACKET WILL BE USEDt SITTING (SIT) OR SUPINE (SUP) » 
380 INPUT BRACK* 
390 ' 

600 9 IF THE SITTING BRACKET IS USED SET FFLAG TO 1. IF THE SUPINE BRACKET 
610 ■ IS USED SET FFLAG TO 0. 

620 * 

^0 IF BRACK*-" SIT- THEN FFLAG- 1: GOTO 630 
640 IF BRACK*»"SUP M THEN FFLAG-0 ELSE GOTO 360 
630 ' 
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*dO PRINT 

670 PR I NT "TURN ON DIGITIZER AND SPARK GAP MULTIPLEXER." 
aSO PRINT 

«*0 PRINT-NOTE: THE DIGITIZER /COMPUTER COMMUNICATIONS MAY FAIL INITIALLY." 
700 PRINT"DO NOT WORRY - THIS IS NORMAL. JUST RETYPE >FPL' WITHOUT TURNING" 
710 PR I NT "OFF THE DIGITIZER AND YOU SHOULD HAVE NO PROBLEM." 
"20 • 

77-0 > CALL SUBROUTINE FPLANE TO DETERMINE THE OBLIQUE SPARK GAP COORDINATES 
740 * OF THE FOCAL, NORMAL AND ORIENTATION POINTS 
730 ' 

7oO G0SU8 1170 ? CALL FPLANE 
770 ' 

780 ' NOW THAT THE RELATIONSHIP BETWEEN THE SPARK GAPS AND THE FOCAL PLANE 
790 » HAS BEEN ESTABLISHED, STOP THE PROGRAM. 

300 *• 
310 STOP 
320 * 

330 7 SUBROUTINE: DISTANCE 
340 ? 

330 * THIS SUBROUTINE CALCULATES THE DISTANCE FROM THE FOCAL POINT IN OBLIQUE 
360 " SPARK GAP COORDINATES AND THE MICROPHONES. 

370 * 

330 f INPUTS: FOCAL < 1 • COORD) - OBLIQUE SPARK GAP COORDINATES OF FOCAL POINT: 

3*0 - 3LCP0INT) - SLANT RANGE DISTANCE TO MICROPHONES: 

900 f CW - COSINE OF ANGLE W. S2-S1-S3 (SPARK GAPS) : 

910 f S2W - SINE SQUARED OF ANGLE W, S2-S1-S3 (SPARK GAPS) 

920 ' OUTPUT: D - DISTANCE FROM FOCAL POINT TO MICROPHONE 

970 ' 

940 * POINT 1 - OBLIQUE MICROPHONE COORDINATES 
930 * 

9©0 ' CALCULATE PROJECTIONS OF SLANT 1 ON X AND Y SPARK AXES 
°70 XPl-(SLU>~2*S3i~2-SL(3)~2)/(2«S31) ' EQN. 6.3 
980 YP1«(SL<1 >~2+S12~2-SL (2)^2)/ (2SS12) ' EQN. 6.4 
990 ' 

1000 * CALCULATE OBLIQUE COORDINATES OF POINT 1 
1040 X 1» < XP1-YP1 *CW) /S2W 9 EQN. 6.8 

1020 V1«(YP1-XPUCW)/S2W 9 EQN. 6.9 

1030 Zi-(SL(l)~2-Xl~2-Yl~2-2«XHYltCW)~.3 ' EQN. 6.12 

1MA0 * 

1030 ' POINT 2 - FOCAL POINT COORDINATES 
10*0 ' 

1070 X2-F0CALU, 1) 
1080 Y2-FQCAH1, 2> 
1090 Z2-F0CAL(l,3> 
1100 » 

1110 ' CALCULATE DISTANCE D USING OBLIQUE DISTANCE FORMULA, EQN. 6.14 
1120 9 

1 1 30 D« ( X2-X 1 ) ~2+ ( Y2-Y 1 ) ~2+ ( Z2-Z 1 ) ~2+2* < X2-X 1 ) * ( Y2-Y 1 ) »CW 

1140 • 

1130 RETURN 

U60-' 

1170 ' SUBROUTINE: FPLANE 
1180 ' 

1190 ' THIS SUBROUTINE DETERMINES THE OBLIQUE SPARK GAP COORDINATES OF THE 
1200 9 FOCAL, NORMAL AND ORIENTATION POINTS FOR THE SITTING OR SUPINE SPARK 
1210 ' BRACKET. 
1220 ' 

1230 9 OUTPUT : FOCAL (POINT, COORD) - OBLIQUE SPARK SAP COORDINATES OF THE FOCAL. 

1240 9 NORMAL AND ORIENTATION POINTS 

1230 9 SUBROUTINES NEEDED: SRINP, NEWTON, STORE, DISTANCE 

1260 ' 

1270 PRINT 

1280 PRINT-ANCHOR THE FOURTH SPARK GAP UNDER THE MICROPHONES SUCH THAT IT CAN" 
1290 PR I NT "BE EASILY FOCUSED UPON. CONNECT THE SPARK GAP TO THE * SPARK GAP 1'" 
1300 PR I NT "OUTPUT ON THE SPARK GAP MULTIPLEXER. MAKE SURE MICROSCOPE IS OUT OF" 
1310 PR I NT "THE WAY. DO NOT MOVE OR DISTURB THIS SPARK GAP WITH RESPECT TO THE" 
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1320 PR I NT "Hi CR0PH0NES ONCE IT HAS BEEN ACTIVATED 'FIRED) BY THE DIGITIZER!" 
1330 9 

1340 * SET SPARK GAP FLAG TO 1 INDICATING ONLY ONE SPARK GAP MUST BE FIRED 
1330 SFLAG* 1 
1360 9 

1370 ' FIRE FOURTH SPARK GAP TO DETERMINE THE DISTANCES BETWEEN THE FOCAL FOIN" 
1380 ' AND THE MICROPHONES. 
1390 9 

1400 GOSUB 4070 ' CALL SRINP 
1410 9 

1420 * STORE THE SQUARED SLANT RANGES IN ARRAY DIST TO BE USED BY NEWTON 

1430 FOR MK-l TO 3 ' STEP THROUGH EACH MICROPHONE 

1440 DIST (MK,1>- (SRMEAN (MK,1>>~2 » MICROPHONE - FOCAL POINT DISTANCES 

1430 NEXT MK * NEXT MICROPHONE 

I4o0 ' 

1470 PRINT ■ _ 

1430 PRINT-CONNECT ALL THREE SPARK GAPS TO THE SPARK GAP MULTIPLEXER AND THEN" 
1490 PPINT"FOCUS THE MICROSCOPE, AT THE HIGHEST MAGNIFICATION (2.3), ON THE" 
1300 PRINT-TIP OF THE FOURTH SPARK GAP, MAKING SURE THE FOURTH SPARK GAP IS" 
1310 PR I NT "NOT DISTURBED." 

1 " 

1330 ' COLLECT SLANT RANGES 'FOR FOCUSED MICROSCOPE 
1340 GOSUB 4070 * CALL SRINP 

I3«y * CALCULATE THE OBLIQUE SPARK GAP COORDINATES OF THE THREE MICROPHONES 
1370 " GIVEN THEIR SLANT RANGES. 
1380 ? 

1390 FOR MK-1 TO 3 ' STEP THROUGH EACH MICROPHONE 
laOO ' 

1610 ' CALCULATE PROJECTIONS OF SLANT RANGE 1 ON THE X AND Y SPARK GAP AXES 
lto20 " 

lo30 XP- ( SRMEAN C MK , 1 ) ~2+S3 1 ~ 2 -SRMEAN ( MK , 3 ) ~2 > / < 2* S3 1 > ' EON. 6.3 
1640 YP- (SRMEAN (MK. 1 ) ~2+S12~2-SRMEAN (MK, 2)~2>/ (2*Sl2) ' EQN. 6.4 
1650 * 

1660 " CALCULATE THE OBLIQUE COORDINATES OF MICROPHONE MK 
1670 FX (MK, 1 ) ■ (XP-YP»CW) /S2W ' EQN. 6.8 
1630 FX iMK, 2) ■ ( YP-XP*CW) /S2W * EQN. 6.9 

1690 FX(MK.3)-(SRMEANtMK,l)-2-FX(MK,l)-2-FX(MK,2)-2-2lFX(MK,l)»FX(MK,^)*CW)«. 

-1700 7 EQN. 6.12 
1710 ' 

172*5 NEXT MK * NEXT MICROPHONE 

1^40 ' CALCULATE COORDINATES OF FOCAL POINT IN SPARK GAP COORDINATES BY SOLVIh 
1730 ' THREE NONLINEAR DISTANCE EQUATIONS FOR. THE THREE COORDINATES. 

1770 » CALL SUBROUTINE NEWTON WITH NFLAG-1 WHICH WILL TAKE THE OBLIQUE 

1780 • SlCRoSoNE COORDINATES STORED IN ARRAY FX, AND THE DISTANCES STORED IN 

1790 ' ARRAY DIST AND SOLVE EQN. 6.13. 

1800 ' 

1810 NFLAG-1 ' SET FLAG TO 1 SUCH THAT NEWTON SOLVES EQUATIONS 6.1^. 
1820 GOSUB 3330 * CALL NEWTON 
1830 ' 

1340 ■ STORE THE OBLIQUE COORDINATES OF THE FOCAL POINT, FROM SUBROUTINE NEWT( 
1830 • IN ARRAY FOCAL. 
1860 * 

1870 FOR COORD-1 TO 3 ' STEP THROUGH X,Y,Z 

1880 FOCAL ( I f COORD) -NSOL ( 1 , COORD) ' STORE COORDINATE 
1390 NEXT COORD ' N E* T COORDINATE 

1900 * 

1910 * COLLECT THE SLANT RANGE DATA FOR THE GRID COORDINATE SYSTEM (FIG 9.1). 
1920 ' THIS DATA IS THE DISTANCES FROM EACH OF THE THREE SPARK GAPS TO 
1930 9 THE THREE MICROPHONES FOR GRID POINTS 1, 2 AND 3 (NORMAL TO THE 
1*40 9 2 AXIS OF THE GRID COORDINATE SYSTEM), AND THE NORMAL POINT ALONG 
1930 ' THE OPTICAL AXIS. 
1960 * 

I Q7Q PR I NT "LOCATE THE GRID COORDINATE SYSTEM WITH THE PROTRACTOR. RESET THE" 
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1980 PRINT "PROTRACTOR TO 0 DEGREES AND ANCHOR IT UNDER THE MICROSCOPE." . 
1990 PRINT 
2<X>'"> " 

2010 ' STORE FIXED CARTESIAN GRID COORDINATES 
2020 GRIDU.D-OJ 6R ID < 1 , 2) -10. 246» GRID 1,- -0 
2030 GRID(2.l)-lO.U4i GRID(2,2)-0i GRIDU,„>-0 
2040 GRID (3. 1 > "Oi GRID(3,2)«Oi GRID<3,*)»0 

UTo FOR PT-1 TO 4 ' STEP THROUGH EACH POINT 

2070 IF PT-1 THEN PRINT-FOCUS ON PT. I" 

51? S J?2 tSn pSiSr-pSSi Z I'nORMAL TO GRID WITH THE CROSSHAIRS" 

2^00 F R3 2£! RlS-IFSt "ICROSCOPE ALIGNED WTHTWWOF THE" 

Tc- bt--«- TUffw pot NT "PROTRACTOR SET AT 0 DEGREES. MAKE 5UK£ TMht mt 
~X IF TO ™!£ PR I NT "FOCUS I NG KNOB IS TURNED SUCH THAT THE MICROSCOPE IS' 
£ J£i ™2J PR INT" AS CLOSE TO THE GRID COORDINATE SYSTEM AS POSSIBLE. 

iUo if Prlt ™ p^nt-m^mSroscope along w^™-™™" GRID " 

2130 IF PT-4 THEN PRINT*BY TURNING THE MICROSCOPE FOCUSING KNOB. 

2160 ' 

2170 ' COLLECT SLANT RANGE DATA 

2 ISO GOSUB 4070 * CALL SRINP 

'fqp mk-1 TO 3 ' STEP THROUGH EACH MICROPHONE 

"i . f6* SPK-l TO 3 * STEP THROUGH EACH SPARK GAF 

"--r, • STORE SLANT RANGE DATA IN ARRAY SDATA 

^1 N |;? T ^ T '" K ' SPK, " SR " EANmK;S Nl! t T SPARK GAP 
-"ir N«T MK ' NEXT MICROPHONE 

s£ nSt p? ' N « T BRID P0INT 

2280 ' CALCULATE DISTANCES FROM GRID POINTS 1, 2 AND 3 TO MICROPHONES 

ZZ™ ' pTm i TO 3 ' STEP THROUGH EACH GRID POINT 

F ?5 R IkIi ?0 3 ' STEP THROUGH EACH MICROPHONE 

=S£ FOR SPK-I TO 3 ' STEP THROUGH EACH SPARK GAP 

2330 ' STORE SLANT RANGES IN VECTOR SL 
2340 SUSPK) -SDATA (FT. MK.SPK) 
2330 NEXT SPK ' NEXT SPARK GAP 



2360 

2380 - 

2390 ' STORED IN VECTOR SL. 

2410 GOSUB 830 ' CALL DISTANCE 

2420 
2430 



STORE SQUARED DISTANCES BETWEEN FOCAL POINT AND MICROPHONES IN ARRAY 
2440 » DIST FOR SUBROUTINE NEWTON. 
2 4 SO ' 

24o0 DIST<PT,MK)«D « mu «K,r 

2W0 NEXT PT ' NEXT 6RID POINT 

2360 f 

2370 ' CALCULATE MICROPHONE COORDINATES 
2330 ' 

2390 PR I NT "CALCULATING MICROPHONE COORDINATES. • . " 

2600 * ^ . ai- T ft * etirw THAT NEWTON WILL SOLVE EQUATIONS 6. IS- 

SSI**/ CA^NeStSn^S^E USTtW GRID MICROPHONE COORDINATE 

2630 ' 
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'640 ' STORE MICROPHONE COORDINATES 

SJX : SJ^S^oTooSSK^ °TO T S?H S'M THREE mCPOPHO.ES =v 
2750 ' CALLING SUBROUTINE DISTANCE. 

I7o ' : ' ' , , T „ - ' STEP THROUGH EACH MICROPHONE 

2770 FOR MK-1 TO * _ , THROUGH EACH SPARK GAP 

2780 FOR SW-l TO j , SLANT RANGE DATA IN VECTOR SL 

N |^ SF s ^ SDATft,4 ' nK ' SPK> 

?S ' ChLL SUBROUTINE DISTANCE TO DETERMINE THE DISTANCES BETWEEN THE 
2S30 ' NORMAL POINT AND. ' THE MICROPHONES. 

2930 GOSUB 330 ' CALL DISTANCE 

' STORE THE SQUARE OF THE DISTANCE IN ARRAY DIST 

"*3S0 " 

~'S°0 DIST(f1K. X)-D 

- Q °° ' „ • NEXT MICROPHONE 

2910 NEXT MK 

2920 ' rflnBnTKJATr , TN A o RAY FX FOR SUBROUTINE NEWTON 

2930 * STORE MICROPHONE COORDINATES IN ARKht run 

2O40 ? , SX£P THROUGH EACH MICROPHONE 

2930 FOR MK-1 TO „ p 2 £p THROUGH X,Y,Z 

2980 NEXT COORD , | MICROPHONE 

2990 NEXT MX 

~0O0 ' _ nf uatcc in THE FIXED GRID COORD. SYSTEM 

=010 ' CALCULATE THE NORMAL POINT COORDINATES IN THE FIXfcii * 

-So PRINT-CALCULATING NORMAL POINT COORDINATES. 

K£ NFL AG* 2 • SET NFLAG TO 2 SUCH THAT NEWTON WILL SOLVE EQUATIONS 6.16 
30*0 GOSUB 5330 ' CALL NEWTON 

TO > STORE NORMAL POINT GRID COORDINATES IN VECTOR NPOINT 

3160 ' POINT.) WITH THESE CO^""™ IN ARRAY ^ EQUATION 6. IS TO THE 

31 70 ' MICROPHONE TO THE «IWAL MINT SpO I NT, SUBROUTINE NEWTON WILL 

COORDINATES IN ARRAY MIKE * E "0R NPOINT^ *u 6.17 FOR THE OBLIQUE 

SOLVE THE THREE NONLINEAR EQUATIONS OF EQUATION e. 
SPARK GAP COORDINATES OF THE NORMAL POINT. 

FOR MK-1 TO 3 ' STEP THROUGH EACH MICROPHONE 

gS ' • 'CALCULATE PROJECTION OF SLANT RANGE 1 ONTO THE X AND Y SPARK AXES. 

§S XP-<SDATA< 3 ,M K ,n T S3i : ^ ™\ fcj 

--> 7Q Y P-<SDATA<3,nK,l)~2*Sl2~2-SDATA<-,MK.2> . -J / i- 

SS CALCULATE THE OBLIQUE SPARK GAP COORDINATES OF EACH MICROPHONE AND 



3190 
3190 
3200 
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ZZVO ' STORE IN ARRAY FX FOR SUBROUTINE NEWTON. 
3310 

3320 FXd-IK. l)»iXP-YptCW) /S2W 9 EQN. 6-3 

^fr.KS:i£SS^««™:.?4«.«.2.-««.«...«««.»«..-.« 

3330 ' EQN. 5.12 

7370 ' CALCULATE THE DISTANCE BETWEEN THE NORMAL POINT AND EACH MICROPHONE IN 
"BO ' THE GRID COORDINATE SYSTEM. 

5S DIST (MK, 1 ) - (MIKE <MK, 1 ) -NPQINT < 1 > > ~2+ (MIKE (MK, 2) -NPOINT (2) ) -2* (MIKE (MK, 3) -NP 
OINT (~) >~2 

3420 NEXT MK ' NEXT MICROPHONE 

3440 • SOLVE FOR THE OBLIQUE COORDINATES OF THE NORMAL POINT BY CALLING 
-430 ' SUBROUTINE NEWTON. 

C.4T?. NFLAG-l * SET NFLAG TO 1 SUCH THAT NEWTON .SOLVES EQUATIONS 6.17. 
t.4S0 GOSUB S530 * CALL NEWTON 

3300 ' STORE OBLIQUE NORMAL POINT COORDINATES IN ARRAY FOCAL. 

-li'r, po* roOPD-1 TO ~ ' STEP THROUGH X.Y.Z 

3330 FOCAL (2. COORD >«NSOLU. COORD) * STORE COORDINATE 
33^0 NEXT COORD ' NEXT COORDINATE 



3360 PRINT-CALCULATING ORIENTATION POINT COORDINATES..." 

3330* CALCULATE THE COORDINATES OF THE PROJECTED GRID POINT 1 ONTOTHE FOCAL 
3390 ' PLANE - WHICH IS THE ORIENTATION POINT, STORED IN VECTOR OPOINT. 
3600 ' 

stis ^sji'.ssssrs^ass'iSrss . . . => ,,, «... , . , » . /« 

3630 ' EQN. 6.26 

zt£Z FOR COORD-1 TO 3 n rnft ^ ' STEP THROUGH X,Y.Z 

3o60 OPOINT (COORD) -NPOINT (COORD) «T*GRID(1, COORD) pnnRD I NATE 

3*70 NEXT COORD NEXT C°° RDIN « T6 

Hit • AS WITH- THE NORMAL POINT, THE OBLIQUE SPARK COORDINATES OF THE 

3700 ' ORIENTATION POINT MUST BE KNOWN. FIRST THE OBLIQUE COORDINATES OF THE 

--Jrt • MICROPHONES MUST BE DETERMINED" WITH THE MICROSCOPE FOCUSED ON THE GRID 

--fi • ORIGIN THEN THE DISTANCES BETWEEN THE MICROPHONES AND THE ORIENTATION 

37^ ■ S Xt'AusTK CALCULATED^ AND THIS «-™™ ««™ ™J 

3740 ' NEWTON TO CALCULATE THE OBLIQUE SPARK GAP COORDINATES OF THE ORIENTATION 

3730 ' POINT. 

3770 FOR HK-1 TO 3 ' STEP THROUGH EACH MICROPHONE 

' CALCULATE THE DISTANCE BETWEEN THE ^^l^^^ t ^J^.^SS!!^S» 
3300 DIST(nK,l)»(MIKE(MK,l)-0P0INT(i))-2*(MIKE<MK,2)-0P0INT(2)) 2* (MIKE (MK. .*) —OP 

0INT(3>)~2 

3310 NEXT MK ' NEXT MICROPHONE 

§830 * SUBROUTINE NEWTON WILL THEN SOLVE EQUATIONS 6.27 FOR THE OBLIQUE 
3340 ' SPARK GAP COORDINATES OF THE ORIENTATION POINT. 

ISm NFLAG-l ' SET NFLAG TO 1 SUCH THAT NEWTON WILL SOLVE EQUATIONS 6.27. 
3370 GOSUB 3530 ' CALL NEWTON 

3890 ' STORE OBLIQUE ORIENTATION POINT COORDINATES IN ARRAY FOCAL 

3??0 FOR COORD-1 TO 3 ' STEP THROUGH -X ; Y P Z 

39-70 FOCAL (3, COORD) -NSOLd. COORD) * STORE COORDINATE 
3930 NEXT- COORD ' NEXT COORDINATE 
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394U 

"950 PRINT 

3960 PR I NT" THE OBLIQUE COORDINATES OF THE FOCAL, NORMAL AND ORIENTATION" 
3970 PRINT"POINTS HAVE BEEN DETERMINED AND STORED IN ARRAY - " ; BRACK* 
3930 PRINT 

3990 PR I NT "READY TO RUN FILEi REG" 

4000 ' STORE FOCAL POINT DATA IN FILE SIT IF FFLAG-1 OR SUP IF FFLAG-0 
4010 ' 

4020 IF FFLAG-1 THEN GOSUB 7230 * STORE SIT 
4030 IF FFLAG-0 THEN GOSUB 7410 ' STORE SUP 
4040 ' 
4030 RETURN 
4060 * 

4070 SUBROUTINE: SRINP 
4080 * 

40*0. > THIS FILE OPENS THE RS-232 PORT FOR DIGITIZER COMMUNICATIONS. 
4100 ? CONTROLS THE SPARK GAP MULTIPLEXER, STORES THE SLANT RANGES 
4110 ■' IN AN ARRAY. AND DETERMINES THE STANDARD DEVIATIONS AND MEANS . 
4120 ' 

4130 ? INPUT: MNOT - IGNORED MICROPHONE 

4140 9 OUTPUT: SRMEAN (MIKE, SKGAP) - AVERAGE SLANT RANGES FOR EACH SPARK GAP 
4130 » 
416U PRINT 

-U"0 PRINT "»tJi»PRESS ANY KEY WHEN FOCUSED/READY* ***t «• 
4 ISO PRINT 
4190 ' 

4200 ' PAUSE UNTIL ANY KEY IS DEPRESSED 
4210 P*»INKEY»: IF PS-"" THEN 4210 
4220 ' 

4230 IF SFLAG-1 THEN SMAX-0 ELSE SMAX-2 ? SET NUMBER OF SPARK GAPS 
4240 * 

4^50 FOR SGAP-0 TO SMAX ' SET MULTIPLEXER THROUGH PARALLEL PORT 
4260 ' CHANGE SGAP TO SGAPl FOR PROPER MULTIPLEXER SEQUENCING 
4270 IF SGAP«0 THEN SGAP 1-0 
4280 IF SGAP-1 THEN SGAP I -2 
4290 IF SGAP-2 THEN SGAPl- 1 
4300 * 

4310 " WRITE THE BINARY CODE FOR SGAPl TO THE PARALLEL PORT ADDRESS 3BC HEX 
4320 OUT &H3BC, SGAPl 
4330 ' 

4340 ' OPEN RS-232 PORT. SET BAUD RATE AND PARITY 
4330 OPEN "C0M1: 9600, 0,7, I" AS #1 
4360 ' - 

4370 FLAGM " ' INITIALIZE CHARACTER FLAG 

4380 * WHILE FLAG 9 SEARCH CHARACTER INPUT FOR ASCII LINE FEED IF FLAG-1 
4390 S*-INPUT«(1, «1) 

4400 IF ASC<S»>-10 THEN FLAG-0 ELSE FLAG-1 

4410 WEND 
4420 » 

4430 W«-INPUT«(26,#1> ' IGNORE FIRST VALUES (FROM LAST SLANT RANGES) 
4440 • 

4430 FOR 1-1 TO 100 9 FIRE EACH SPARK GAP 100 TIMES- 

4460 • 

4470 ' MONITOR COMM. BUFFER, IF >40 CHARACTERS IN BUFFER TURN OFF DIGITIZER 
.4480 ' BY WRITING BINARY CODE 10 TO ADDRESS 3FC HEX - THE RS-232 PORT 
4490 9 

4300 IF L0C(1)>40 THEN OUT &H3FC, 10 ELSE OUT 3<H3FC,ll 
4310 * 

4320 ' INPUT SLANT RANGES AND INSERT DECIMAL POINTS TO READ DISTANCES IN CM. 
4330 * 

4340 A*-INPUT*(26,»1> * INPUT 26 CHARACTER STRING (4 SLANT RANGES) 

4360 C«-MID*(A«,1,6): CL«-LEFT« (CS, 3) 1 CR«-RIGHT* <C», 3) : CC*«CL*+" . •♦CRs 

4370 DS-MID*<At,7,6>: DLS-LEFT* (D*, 3) : DR««RIGHT* (D«, 3) : DD*«DL*+" . "*DRS 

4380 E«-MIDS(A«,13,6>: EL*»LEFT» (ES. 3) : ER*-RIGHT* <E«, 3) : EE*-EL«+" . "+ERs 

4390 F**MIDs(As, 19.6) : FL»-LEFTs (F», 3) : FR*-R I GHT* (Ft , 3 > : FF*-FL*+" . "*FRs 
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4600 ' 

4610 ' STORE THE VALUE OF EACH CHARACTER STRING FOR EACH SLANT RANGE IN Gl-4 
4620 Gl-VALCCCt): 52-VALtDD*): G3-VAL(EES): G4-VAL(FF*) 1 61 

4630 f 

4640 f IGNORE SLANT RANGE DATA FROM MICROPHONE MNOT AND STORE REMAINING 
4650 ' SLANT RANGES IN Pl-3 ^ iNQ 

4660 ' 

4670 IF MNOT-1 THEN P1-G2: P2»G4: P3-G3 
4630 IF MNOT-2 THEN Pl-Gli P2-G3l P3-G4 
4690 IF MNOT- 3 THEN Pl«G4s P2«G2j P3-G1 • 
4700 IF. MN0T«4 THEN P1-G3: P2-GU P3-G2 

4710 ' 

4720 ' STORE SLANT RANGES IN ARRAY SLANT 

4730 SLANT (1,1) "Pi : SLANT < 1 , 2) -P2i SLANT (1,3) «P3 

4740 * 

4730 NEXT I * GET NEXT SET OF FOUR SLANT RANGES 
4760 • 

4770 CLOSE #1 > SUPPRESS DIGITIZER COMMUNICATION 
4730 " 

4790 • CALCULATE MEAN SLANT RANGE VALUES 
4900 ? 

4310 C1«0: C2-0: C3-0 • INITIALIZE SLANT RANGE SUM 

4S20 BAD 1-0: BAD2-0: BAD3-0 ' INITIALIZE BAD DATA COUNTER 

4330 CNT1-0: CNT2*0: CNT3-0 ? INITIALIZE NUMBER OF SLANT RANGES 

4840 " 

43S0 FOR I»l TO 100 ' STEP THROUGH EACH SLANT RANGE 

4360 ' 

4370 " TEST EACH SLANT RANGE. IF LESS THAN MAX VALUE OF 243 CM. ADD TO SUM. 
4880 ' IF A VALUE IS GREATER THAN 243 CM.. INCREMENT BAD COUNTER AND IGNORE 
4890 > THE BAD DATA. IF 3 BAD VALUES ARE DETECTED, THE MICROPHONE IN ERROR 
4900 ? MIGHT BE BLOCKED, AN ERROR MESSAGE IS DISPLAYED AND SRINP IS RECALLED. 
4910 ' 

4920 IF SLANT (1,1) <243 THEN Cl-Cl +SLANT ( I , 1 ) i CNT1-CNT1+1: ELSE BAD 1 "BAD 1 1 
4930 ' 

4940 IF BADl.-=3 THEN PRINT "BAD DATA FROM X MIKE-TRY AGAIN, "i GOTO 4170 
4930 ' 

4960 IF SLANT (I, 2X243 THEN C2-C2+SLANT ( 1 , 2) I CNT2-CNT2+1: ELSE BA02-BAD^1 
4970 " 

4980 IF BAD2>«3 THEN PRINT "BAD DATA FROM ZERO MIKE-TRY AGAIN. M : GOTO 4170 
4990 9 

3000 IF SLANT i 1 ,3X243 THEN C3«C3+SLANT ( I , 3) : CNT3-CNT3+1: ELSE BAD3-BAD3+1 
3010 * 

3020 IF BAD3>-3 THEN PRINT "BAD DATA FROM Y MIKE-TRY AGAIN.": GOTO 4170 
5030 ' 

3040 NEXT I > TEST NEXT SET OF SLANT RNAGES 
3030 ' 

3060 9 CALCULATE. THE AVERAGE SLANT RANGE VALUES 
3070 SMEANI-C1/CNT1: SMEAN2-C2/CNT2: SMEAN3-C3/CNT3 
3080 9 

3090 ' CALCULATE STANDARD DEVIATIONS OF SLANT RANGE VALUES TO DETERMINE IF 
3100 9 THERE IS UNACCEPTABLE VARIATION IN THE DATA, IE. A AIR DISTURBANCE; 
31X0 ' 

3120 Cl-Oi C2-0l C3-0 • INITIALIZE SLANT RANGE SUM 

3130 CNT1-0: CNT2»0j CNT3-0 9 INITIALIZE NUMBER OF SLANT RANGES 

3140 ' 

3130 FOR 1-1 TO 100 9 STEP THROUGH EACH SLANT RANGE 

3160 » 

3170 IF SLANT (1,1) <243 THEN C1-C1+ (SLANT < I , 1 > -SMEAN1 > ~2: CNT1-CNT1 + 1 
3180 IF SLANT (I, 2X243 THEN C2-C2+ (SLANT (1,2) -SMEAN2) ~2: CNT2-CNT2+1 
3190 IF SLANT (I, 3X243 THEN C3-C3+ (SLANT ( 1 , 3) -SMEAN3) ~2: CNT3-CNT3+1 
3200 ' 

3210 NEXT I ■ TEST NEXT SET OF SLANT RNAGES 
3220 ' 

3230 ' LET Nl-3 EQUAL THE TOTAL NUMBER OF VALUES-1 
3240 Nl-CNTl-1: N2-CNT2-1: N3-CNT3-1 
3230 " 
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5160 > CALCULATE STANDARD DEVIATIONS 

5270 SD1 = (Cl/Nl)".5: SD2- (C2/N2) ~. 3: SD3« (C3/N3) ^ 

5290 ' COHPARE STANDARD DEVIATIONS TO SDLIMIT TO DETERMINE SLANT RANGE 
3300 9 ACCEPTABILITY. 
3310 ' 

3320 SDLIMIT-. 1 9 SET STANDARD DEVIATION LIMIT 

sSS V SD1>SDLIMIT THEN PRINT M DATA OUT OF RANGE - TRY AGAIN. -, GOTO 4170 
IF SD^>SDLIMIT THEN PR I NT "DATA OUT OF RANGE - TRY AGAIN. 3 GOTO 41,0 
3360 S SK^!!" THIN pSSnT'DATA OUT OF RANGE - TRY AGAIN. H : GOTO 4170 

539^ ' ADD THE COUNTER DELAY EQUIVALENT OF 4.43 CM TO *™ i™* SLANT RANGE 
S390 SMEAN1-SMEAN1+4.45: SMEAN2-SMEAN2+4. 45: S>lEAN^«SMEAN3*4. 4., 

3410 SP-SGAP+1 f INCREMENT 3RMEAN ARRAY POINTER 

5430 * STORE MEAN VALUES IN ARRAY SRMEAN ( M I KE# « SP ARK# ) 

S«0 SRMEAN ( 1 , SP> -SMEAN1 : SRMEAN (2, SP> -SMEAN2: SRMEAN (3, SP) -SMEAN3 

FIRE NEXT SPARK 

RESET SFLAG TO 0 



3460 

3470 NEXT SGAP ' FIRE NEXT SPARK GAP 



3-T*0 3FLA5=>: 
3300 ' 
3310 RETURN 
33^0 * 

3330 y SUBROUTINE: NEWTON 

«S ' THIS SUBROUTINE SOLVES 3 SETS OF 3 NONLINEAR EQUATIONS FOR 3 UNKNOWNS. 

' USIU^ METHOD. ™J« V ^"5 £ 6 13 

3370 ' NFLAG, CW1 IS SET TO ZERO OR TO CW SUCH THAT EQUATIONS 6. lw OR 6.1- 
--on > CAN BE SOLVED , RESPECTIVELY. 

3390 ' THE OUTPUT IS THE COORDINATES STORED IN ARRAY NSOL. 
3ol0 ? INPUTS: FX (POINT. COORD)- P I ^.COORDINATES^ 



Sft-v, • *DIST(P0INT,3) - DISTANCES BETWEEN POINTS; 

-*"o > CW- COSINE OF ANGLE U. S2-S1-S3 (SPARK GAPS); 

• NFLAG - 1 TO SOLVE EON. 6. 13, 0 TO SOLVE EQN. 6.1. 



J630 ? OUTPUTi NSOL ( PO I NT , COORD ) - COORDINATE SOLUTIONS 
3660 * SUBROUTINE NEEDED i INVERSE 

MaS ' EVALUATE THE STATUS OF NFLAG AND SCT APPROPRIATE VARIABLES. 
3690 > EQNUM - NUMBER OF SETS OF EQUATIONS TO BE SOLVED | 
5700 ' INIT - INITIAL VALUES FOR ITERATIVE SOLUTION. 
3710 * 

3720 IF NFLAG-0 THEN CWl-Ot EQNUM-3i INIT-15 
3730 IF NFLAG- 1 THEN CWl-CWi EQNUM- 1: INIT— ^0 
3T40 IF NFLAG -2 THEN CWl-Oi EQNUM-lJ INIT-lO 
5730 * 

3760 ' SET ITERATIVE TOLERANCE 
5770 E-.001 
*7S0 * 

5790 9 LOOP THROUGH EACH SET OF EQUATIONS UP TO EQNUM 
3800 FOR MK-1 TO EQNUM 
3S10 » 

5820 ' STORE INITIAL SOLUTIONS 

5S30 FOR COORD-i TO 3 ' STEP THROUGH X,Y,Z 

5840 M#(1,C00RD)-ITER » STORE IN I TAL SOLUTION 

5830 NEXT COORD ' NEXT COORDINATE 

5860 ' 

5870 ? LOOP THROUGH UP TO 33 ITERATIONS 
3880 FOR ITER-2 TO 33 
5890 ' 

5900 ' STORE ITER-1 SOLUTION 

5910 FOR COORD- 1 TO 3 ' STEP THROUGH X,Y P Z 
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«S ^ CALCULATE J ACQS I AN MATRIX COEFFICIENTS FRC* THE DERIVATIVES CF THE 

3960 ' EQUATIONS. 
3970 ? 

3990 FOR 1-1 TO 3 

5990 FOR J-l TO 2 

6000 IF J»l THEN Jl"2 

IS? I a.a: J IS(;i:la»-FX(i, J ))* 2 »c W i.cnv ( i( J i)-Fx<i, J i)) 

6030 NEXT J 

a040 J#(I.3)-2» (MV«<3>-FX<I.^> > 

6030 NEXT I 
a060 * 

6070 ' INVERT THE JACOBIAN MATRIX 

6080 GOSUB 6420 ' CAUL INVERSE 
oV90 " 

a 100 * DETERMINE SOLUTION ITER 

" -FX (I. 1) >KMV*(2>-FX(I,2) >lCWl-DIST<I,MK> 
6 130 NEXT I 

6170 NEXT I 

t\l°0 : CALCULATE ITERATIVE ERROR AND SEE IF IT IS -ITHIN THE SET TOLERANCE TOL 

6200 FOR I"l TO 3 , Y«rtr© 

6210 T ( I ) -ABS <m t ITER, I ) -MV# < I > ) /ABS CM ( ITER, I > > 

6220 NEXT I 

: IF ERROR IS WITHIN TOL STORE SOLUTION OR CONTINUE TO NEXT ITERATION 
6230 IF <TUXE> AND CTC2XE) AND (TUXE) THEN 6„*0 
6260 ' 

6270 NEXT ITER ' NEXT ITERATION 

lilt I IF NO SOLUTION IS REACHED AFTER 3= ITERATIONS PRINT ERROR nESSASE 
6300 ' AND STORE LAST ITERATION. tTeB _ 

o310 PRINT "NO SOLUTION FOR POINT •■»»<• ITER«~ 

©~**0 * 

6330 ' STORE SOLUTIONS IN ARRAY NSQL- 

o340 FOR I -I TO 3 

6330 NSOL (MK, I ) "M1M ITER, I ) * 

6360 NEXT I 

So NEXT MK ' NEXT SET OF EQUATIONS 

. 6390 ' 

6400 RETURN 
6410 * 

6420 ' FILENAME: INVERSE 

6430 * MATRIX USING THE CROUT ALGORITHM 

6440 > THIS SUBROUTINE INVERTS THE JACOBIAN MATRIX USirua 

6430 * 

6460 ' INPUT! ' J# (3,3) - JACOBIAN MATRIX 
6470 » OUTPUTi 111(3,3) - INVERTED MATRIX 

tSS ' SET INVERSION CRITEREA BY STORING THE IDENTITY MATRIX SOLUTIONS 
6300 FOR 1-1 TO 3 
6310 FOR J"l TO 3 

6320 IF I-J THEN I#CI,J>-1 ELSE I#(I,J>«0 
6530 NEXT J 
6340 NEXT I 

tilt TOUMlE-09 ' SET ZERO TOLERANCE FOR MATRIX SINGULARITY 
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6370 • 

osao for 1-1 to : ? SEARCH for largest element in a column 

6390 X — l 
6600 * 

6610 FOR K-l TO 3 * SEARCH BELOW MAIN DIAGONAL 
o620 IF A8S'.J#CK. I> K-X THEN 6640 

6630 Q-Kt X»ABS(J#<K.I>> * Q IS ROW OF LARGEST ELEMENT 

6640- NEXT K 9 X IS LARGEST ELEMENT 

6630 ' 

6660 9 TEST FOR SINGULAR MATRIX AND STOP' IF NECESSARY 
6670 IF X>»0 THEN 6690 ELSE PR I NT * S I NGULAR MATRIX" : STOP 
6680 * 

6690 IF I=Q THEN 6790 'INTERCHANGE IF NEEDED 

6700 * 

6710 FOR J«l TO 3 ' NO. SWITCH ROWS I AND Q 

6720 T»J#(I.J>: J# ( I « J ) «J# CQ. J ) : J«<Q«J>-T 
o7-0 NEXT" J 
6740 ' 

6730 FOR 0-1 TO 3 

6760 T«HMI,J>i It<I,J>-I#(Q,J>: I#CQ ? J)-T ' SWITCH RIGHT HAND SIDE 
67T0 NEXT J 
©730 ' 

e7«0 FOR J-l TO 3 * ELIMINATE ON THAT ONE ROW 

6800 ' 

6810 IF KJ THEN Ml-X-l ELSE Ml-J-1 

r820 S-0 * PIND INNER PRODUCT OF ROW I AND COLUMN J - S 

6830 ' 

6340 FOR K-l TO Ml 

6830 S-S+J* ( I , K) * J* (K. J ) 

686U NEXT K 

6370 ? 

6880 J#<I. J)-J#(I, J)*S 

6890 ' STOP HERE IF BELOW MAIN DIAGONAL, CHECK FOR SINGULARITY ELSE NORMALIZE 
6900 IF I>«J THEN 6930 

6920 IF ABS<J#(I ? I>XTOL* THEN PR I NT "SINGULAR MATRIX": STOP 
o930 ' 

fl Q^ J*<I.J>— J»<X, J>/J#U, I) 'ELSE NORMALIZE 
6930 NEXT J ' NEXT COLUMN 
6960 NEXT I ' NEXT COLUMN 

6970 ' 

6980 9 REDUCE RIGHT HAND SIDE 

6990 ' ~ 

70'X, FOR J-l TO 3 ' FOR EACH SET OF CONSTANTS 

7010 FOR I- 1 TO 3 * LOOK DOWN EACH COLUMN 

7020 S-0 
7030 ' 

7040 FOR K-l TO 1-1 

7030 S-S+J«U P K)*HMK,J) ' TAKE PARTIAL INNER PRODUCT 
7060 NEXT K 
7070 ' 

7080 It(I, J) — (I»U, J>+S) /J# (I.I) 
7Q90 * 

7100 NEXT I 
7110 * 

7120 FOR 1-3 TO I STEP -1 ' WORK BACK UP COLUMN 

7130 S-0 ' TAKE PARTIAL INNER PRODUCTS 

7140 9 

7150 FOR K-I + l TO 3 
7160 S-S*J*(I,K> *I#<K, J) 
7170 NEXT K 
7180 * 

7190 lit I, J)— IftCX, J)+S 
7200 * 

7210 NEXT I 
722V NEXT J 



7230 

7240 

7220 

7260 

7270 

7230 

7290 

7300 

7310 

7320 

7330 

734«.) 

7350 

7360 

7370 

-380 

7390 

74(.iy 

7410 

7420 

7430' 

7440 

"430 

74*0 

7470 

-480 

7490 

"*300 

7310 

7320 

7330 

7340 

7330 

7360 



RETURN 
r 

' SUBROUTINE: STORE-SIT 

" THIS SUBROUTINE STORES THE OBLIQUE SPARK GAP COORDINATES OF THE FOCAL. 
» NORMAL AND ORIENTATION POINTS IN FILE "SIT". 



' OPEN FILE SIT TO STORE COORDINATES 

9 STEP THROUGH EACH POINT 

' STEP THROUGH X.Y.Z 

' WRITE THE COORDINATE TO THE FILE 

' NEXT COORDINATE 

" NEXT POINT 

? CLOSE FILE SIT 



OPEN "SIT" FOR OUTPUT AS #1 
FOR PNT«i TO 3 
FOR COORD- I TO 3 

PRINT • 1 • FOCAL ( PNT , COORD ) 
NEXT COORD 
NEXT PNT 

r 

CLOSE ttl 
RETURN 

* SUBROUTINE: STORE-SUP 

? THIS SUBROUTINE STORES THE OBLIQUE SPARK GAP COORDINATES OF THE FOCAL, 
? NORMAL AND ORIENTATION POINTS IN FILE "SUP" . 



OPEN "SUP" FOR OUTPUT AS #1 
FOR PNT-1 TO 3 
FOR C00RD«1 TO 3 

PRINT ttl. FOCAL (PNT, COORD) 
NEXT COORD 
NEXT PNT 

CLOSE »1 

RETURN 



? OPEN FILE SUP TO STORE COORDINATES 

' STEP THROUGH EACH POINT 

• STEP THROUGH X.Y.Z 

• WRITE THE COORDINATE TO THE FILE 

• NEXT COORDINATE 

• NEXT POINT 

' CLOSE FILE SUP 
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10 

20 * 

TO * 

40 » 

30 ' 

60 ' 

70 ' 

30 ' 
or. 

100 
110 
120 
130 
140 
ISO 
loO 
1 TO 
180 
190 
200 
210 
220 
230 
240 
230 
260 
270 
290 
290 
300 
310 
320 
330 
340 
330 
360 
370 
380 
390 
400 
410 
420 
430 
440 
430 
460 
470 
480 
490 
300 
310 
320 
330 
340 
330 
560 
370 
380 
390 
600 
610 
620 
630 
640 
650 



APPENDIX C 

THESIS TITLE: REFERENCE - DISPLAY SYSTEM FOR THE INTEGRATION OF CT 
SCANNING AND THE OPERATING MICROSCOPE 

FILENAME: REG 

PROGRAMMER: JOHN F. HATCH 
f DATE: AUGUST, 1984 

» NOTE: REFER TO CHAPTER 6 IN THESIS FOR EQUATIONS, 

> AND REFERENCE 10. FOR ALGORITHMS. 

> THIS FILE IS "THE DRIVER FOR THE REGISTRATION PROCEDURE. 
OPTION BASE 1 ' SET THE LOWER ARRAY INDEX TO 1 FOR CONVENIENCE 



ALLOCATE SPACE FOR ARRAYS 

NOTE: ARRAY NAMES FOLLOWED BY A 



SIGN ARE DOUBLE PRECISION 



DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 
DIM 



J#(3.3> 
It (3.3) 

SRMEAN<3.3) 

SDATA(3,3,3) 

FOCAL (3. 3) 

NS0L(3.3) 

FX (3,3) 

MVt (3) 

M#(100,3> 

SL(3) 

F#(4) 

DIST(3.3> 

D(3) 

T<3) 

MIKE(3,3) 
FP0INTC3) 
FPLANE ( 4 ) 
SLANT (150. 3) 



JACOBIAN MATRIX COEFFICIENTS 
JACOB I AN MATRIX INVERSE 
MEAN SLANT RANGES 

FOCAL, ^ORMAL AND ORIENTATION POINT COORDINATES 

NEWTON SOLUTIONS 

FIXED COORDINATE SYSTEM 

ITERATION VARIABLE CN-1> 

ITERATION VARIABLES 

SLANT RANGE DATA 

ITERATIVE SOLN TO NEWTON 

DISTANCES USED IN NEWTON 

DISTANCES USED IN DISTANCE 

ITERATION ERRORS 

MICROPHONE COORDINATES 

FOCAL POINT COORDINATES 

FOCAL PLANE COEFFICIENTS 

SLANT RANGES 



' MEASURED DISTANCES BETWEEN SPARK GAPS IN CM. 
S12-30.041 ' SPARK GAP 1 - SPARK GAP 2 

S-*3-29.993 ? SPARK GAP 2 - SPARK. GAP 3 

S3 1-29. 637 • 9 SPARK GAP 3 - SPARK GAP 1 

'• REQUEST OPERATIVE POSITION FOR RIGHT/FRONTAL <R1 OR LEFT <LI CRANIOTOMY 
PR I NT 

PRINT-WHICH OPERATIVE POSITION WILL BE EMPLOYED, " 

PRINT" (R) - RIGHT/FRONTAL" 

PRINT- OR- . 

PRINT- (L) - LEFT CRANIOTOMY"; 

INPUT ORPOS« 
# 

» IF "R" THEN IGNORE DATA FROM BLOCKED MIKE 4 
IF ORPOSt-^R" THEN MNOT-4: GOTO 610 
» IF "L" THEN IGNORE DATA FROM BLOCKED MIKE * 
IF ORPQSt«-L - THEN MN0T«2i ELSE GOTO 300 

' REQUEST PATIENT POSITION SITTING (SIT) OR SUPINE (SUP) 
PRINT 

PRINT-WHAT POSITION WILL THE PATIENT BE IN, 1 
PRINT" (SIT) - SITTING" 
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aoO PRINT" or- 

670 PRINT" iSUP) - SUPINE": 

o30 INPUT BRACKS 

6V0 ' 



5« \?^Z,ff^£^?$S- AND 0B1ENT - P0INT C00RDS FRCn S,TTINS 

'60 PRINT 

^70 PRINT "TURN ON DIGITIZER AND SPARK GAP MULTIPLEXER." 
■ 30 PR I NT 
"90 ' 

??? COMPUTER FAILS AFTER REGISTRATION RECALL MICROPHONE COORDINATE* 

310 PRINT-DID THE COMPUTER FAIL AFTER FOCUSING ON THE FIDUCI ALS !S m - 
320 PR I NT " OR NO <N)7 ANSWER 'NO' TO THIS QUESTION IF THIS IS THE FIRST TIME" 
3-0 PRINT-RUNNING THIS PROGRAM. " ; 5 TIf1£ 

340 INPUT FAILS 
330 ' 

ooO HFLAG^O 

370 IF FAILS- "Y" THEN GOSUB 7720: GOTO 900 
SSV IF FAILS' * "N" THEN GOTO 310 
390 • 

900 CFLAG-0 ' SET CALIBRATION FLAG 
910 * 
*20 PRINT 

940 ISK ?S?S Y0U ^ 10 TEST THE REGISTRATI ° N SYSTEM iY OR N>"; 
9 "50 PRINT 

960 PRINT"PRE3S CTRL BREAK TO STOP AT ANY TIME" 
970 PRINT 
*80 * 

?90 IF TSTS="Y" THEN GOSUB 6630t GOTO 920 ' CALL TEST 
1000 IF TSTSO-N" THEN GOTO 920 
1010 » 

10'0 ' DETERP1INE MICROPHONE COORDINATES IN CT COORDINATE SYSTEM 

1040 IF HFLAG-1 THEN GOTO 1080 
1030 » 

1060 GOSUB 4260 » CALL REGISTER ~ 
1070 * 
1080 PRINT 

1090 PRINT - BEGIN PROCEDURE - BREAK TO STOP" 
1100 PRINT 
1110 ' 

1120 ; STORE MICROPHONE COORDINATES (IN ARRAY MIKE) IN ARRAY FX FOR NEWTON 

A lwO 

J J £2 F £™"t£iL T0 3 * STEP THROUGH EACH MICROPHONE 

llwO FOR COORD- 1 TO 3 ' STEP THROUGH X,Y,Z 

1160 FX<MK, COORD ) -MIKE CMK f COORD) ' STORE COORDINATE 

1170 NEXT COORD > NEXT COORDINATE 

1180 NEXT MK > NEXT MICROPHONE 

1190 ? 

1200 ' LOOP THROUGH UP TO 100 RECONSTRUCTED CT IMAGES 
1210 FOR RECON-1 TO 100 
1220 ' 

J 552 \ MLL SUBROUTINE SLICE TO TAKE THE COORDINATES OF THE MICROPHONES IN rr 
240 > COORDINATES AND THE SLANT RANGE DISTANCES AND CALCULATE THE COORD I NATES 

ISO • WsISis: 0 "' P0INT ' EQUATI0N ° F ™ E FGML PLANE AN ° ^ ^EE°S?REC?IoN 
1270 ' 

1280 GOSUB 3260 > CALL SLICE 
1290 * 

J-?? > ^! P b AY THE COORDINATES OF THE FOCAL POINT. THE EQUATION OF THE FOCAL 
1^10 ' PLANE AND THE THREE DIRECTION COSINES. 



1 320 7 
1330 PRINT 

1340 PR I NT "FOCAL POINT COORDINATES CX,Y,Z)" 
1330 PRINT 

1360 PRINT FPOINTU) ,FP0INT(2) ,FP0INT<3) 
1370 PRINT 

1380 PRINT "FOCAL PLANE COEFFICIENTS <A,B,C,D)" 
1390 PRINT 

1400 PRINT FPLANE(l) , FPLANE £2) , FPLANE ( 3 ) . FPL ANE ( 4 ) 
1410 PRINT 

1420 PRINT W X DIRECTION COSINE-" ; DCOSX 
1430 PRINT-Y DIRECTION COSINE-" ; DCOSY 
1440 PRINT-2 DIRECTION COSINE-" ; DCOSZ 
1430 * 

1460 NEXT RECON » REFOCUS FOR NEXT SLICE 

1470 * 

1430 ' STOP AFTER 100 RECONSTRUCTED CT in AGES 
1490 * 
1300 STOP 
1310 1 

1320 ' SUBROUTINE: SRINP 
1330 ' 

1340 r THIS FILE OPENS THE RS-232 PORT FOR DIGITIZER COMMUNICATIONS, 
1330 * CONTROLS THE SPARK! GAP MULTIPLEXER, STORES THE SLANT RANGES 
1360 ' IN AN ARRAY, AND DETERMINES THE STANDARD DEVIATIONS AND MEANS, 
1370 * 

1380 * INPUT: MNOT - IGNORED MICROPHONE 

1390 ' OUTPUT: SRMEAN (MIKE, SKGAP) - AVERAGE SLANT RANGES FOR EACH SPARK GAP 
loOO 

1610 PRINT 

1620 PRINT "XXXXXPRESS ANY KEY WHEN FOCUSED / READY I x x x x H 
1630 PRINT 
1640 ? 

1630 » PAUSE UNTIL ANY KEY IS DEPRESSED 
1660 P«-INKEY«: IF Ps-"" THEN 1660 
1670 ' 

1680 FOR SGAP-0 TO 2 ' SET MULTIPLEXER THROUGH PARALLEL PORT 
1690 ' CHANGE SGAP TO SGAP1 FOR PROPER MULTIPLEXER SEQUENCING 
1700 IF SGAP-0 THEN SGAP 1-0 
1710 IF SGAP-1 THEN SGAP 1-2 
1720 IF SGAP- 2 THEN SGAPl-1 
1730 • 

1740 9 WRITE THE BINARY CODE FOR SGAP1 TO THE PARALLEL PORT ADDRESS 3BC HEX 
1730 OUT &H38C , 5GAP 1 
1760 • 

1770 ' OPEN RS-232 PORT, SET BAUD RATE AND PARITY 
1780 OPEN "C0M1 i 9600,0, 7, 1" AS ftl 
1790 " 

1800 FLAG-1 ' INITIALIZE CHARACTER FLAG 

1810 WHILE FLAG ' SEARCH CHARACTER INPUT FOR ASCII LINE FEED IF FLAG-1 
1 820 St- INPUT* < 1 , # 1 ) 

1830 IF ASC'<SS>-10 THEN FLAG-0 ELSE FLAG-1 

1 840 WEND 
1830 \ 

I860 W«-INPUT«<26,#1) ' IGNORE FIRST VALUES (FROM LAST SLANT RANGES) 
1870 * 

1880 FOR I- 1 -TO 30 9 FIRE EACH SPARK GAP 30 TIMES 

1890 9 

1900 9 MONITOR COMM. BUFFER , IF >40 CHARACTERS IN BUFFER TURN OFF DIGITIZER 
1910 9 BY WRITING BINARY CODE 10 TO ADDRESS 3FC HEX - THE RS-232 PORT 
1920 * 

1930 IF L0C<1)>40 THEN OUT 5cH3FC, 10 ELSE OUT ScH3FC,ll 
1940 ' 

1930 ' INPUT SLANT RANGES AND INSERT DECIMAL POINTS TO READ DISTANCES IN CM. 
1960 9 

1970 A$-INPUT*<26,#1> * INPUT 26 CHARACTER STRING (4 SLANT RANGES) 
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1980 ' 

19*0 C*=MIDs<As,l,6): CL«»LEFT«iC*,3)i CRS-RIGHTS (Cs, 3) : CCs-CLs*" . " +CRS 
2000 D*«MID«CAS,7.6>: DLS-LEFT* (DS, 3) : DRS-RIGHT* (D«, 3) : DD*-DL*+ M . "+DRS 
2010 E*=MID«iA«, 13,6) : EL*-LEFT* (E*. 3> : ER«-RIGHT* <E«, 3) i EE*-EL*+" . "+ERS 
2020 FS«MID*'AS, 19,6) : FL*«LEFT* (F«, 3) : FRs-RIGHTS <F«, 3) : FF*-FU«+". "*FRs 

2030 ? 

2040 ' STORE THE VALUE OF EACH CHARACTER STRING FOR' EACH SLANT RANGE IN Gl-4 
2030 Gl-VAL(CCS): G2«VAL(DD*>: G3-VAL(EE*>: G4«VAL(FF«> 

2060 ' 

2070 ? IGNORE SLANT RANGE DATA FROM MICROPHONE MNOT AND STORE REMAINING 
2080 * SLANT RANGES IN PI -3 

20 Q 0 * 

2100 IF MNOT » 1 THEN P1-G2: P2-G4: P3-G3 
2110 IF MNOT-2 THEN Pl-Gl: P2-G3: P3-G4 
2120 IF MNOT- 3 THEN P1-G4: P2-G2: P3-G1 
2130 IF MNOT -4 THEN P1-G3: P2-G1 2 P3-G2 

2140 ' 

2130 ' STORE SLANT RANGES IN ARRAY SLANT 

2160 SLANT (I, l)«Pl: SLANT ( I , 2> -P2i SLANT (1,3) «P3 

2170 ? 

2180 NEXT I * GET NEXT SET OF FOUR SLANT RANGES 
2 1 ^O " 

2200 CLOSE #1 * SUPPRESS DIGITIZER COMMUNICATION 
22 1 0 * 

2220 ' CALCULATE MEAN SLANT RANGE VALUES 
^"^30 * 

2240 Cl-O: C^*0: C3-0 ' INITIALIZE SLANT RANGE SUM 

2230 BADl-Oi BAD2«0i BAD3-0 * INITIALIZE BAD DATA COUNTER 

2260 CNT1-0: CNT2-0: CNT3-0 ' INITIALIZE NUMBER OF SLANT RANGES 

2230 FOR I- 1 TO 30 ' STEP THROUGH EACH SLANT RANGE 

2290 * 

23<V> ' TEST EACH SLANT RANGE. IF LESS THAN MAX VALUE OF 243 CM. ADD TO SUM. 
-3K» • IF A VALUE IS GREATER THAN 243 CM. , INCREMENT BAD COUNTER AND IGNORE 
-3?0 ? THE BAD DATA. IF 3 BAD VALUES ARE DETECTED, THE MICROPHONE IN ERROR 
2330 ' MIGHT BE BLOCKED • AN ERROR MESSAGE IS DISPLAYED AND SRINP IS RECALLED. 

2330 IF SLANT (I. 1X243 THEN C 1 "C 1 +SLANT (1,1): CNT1-CNT1+1: ELSE BADl-BADl+1 

2360 ■ 

23T0 IF BADl>-3 THEN PRINT "BAD DATA FROM X MIKE-TRY AGAIN.-: GOTO 1620 
2330 ' 

2390 IF SLANT (1,2X243 THEN C2-C2+SLANT ( I , 2) I CNT2-CNT2+H ELSE BAD2-BAD2+1 
2400 * 

2410 IF BAD2>-3 THEN PRINT **BAD DATA FROM ZERO MIKE-TRY AGAIN.": SOTO 1620 
2420 * 

2430 IF SLANT (1,3X243 THEN C3-C3+SLANT ( 1 , 3) I CNT3-CNT3+U ELSE BAD3-BAD3+1 
2440 * 

2430 IF BAD3>-3 THEN PRINT "BAD DATA FROM Y MIKE-TRY AGAIN.": GOTO 1620* 
2460 * 

2470 NEXT I ' TEST NEXT SET OF SLANT RNAGES 
2480 ' 

2490 » CALCULATE THE AVERAGE SLANT RANGE VALUES 
2300 SMEAN1-C1/CNT1: SMEAN2-C2/CNT2: SMEAN3-C3/CNT3 

2320 ' CALCULATE STANDARD DEVIATIONS OF SLANT RANGE VALUES TO DETERMINE IF 
2330 * THERE IS UNACCEPTABLE VARIATION IN THE DATA, IE. A AIR DISTURBANCE. 
2340 * 

~>ZZC Cl-O: C2-0: C3-0 ' INITIALIZE SLANT RANGE SUM 

?360 CNT1-0: CNT2-0: CNT3-0 ' INITIALIZE NUMBER OF SLANT RANGES 

2380 FOR 1-1 TO 30 F STEP THROUGH EACH SLANT RANGE 

2390 ' 

2600 IF SLANT (I, 1X243 THEN C1-C1+ (SLANT < 
2ol0 IF SLANT (I, 2X243 THEN C2-C2+ (SLANT < 
2620 IF SLANT (1,3X243 THEN C3-C3* (SLANT 1 
2630 ' 
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2640 NEXT I ? TEST NEXT SET OF SLANT RNAGES 
2630 > 

2660 " LET Nl-3 EQUAL THE TOTAL NUMBER OF VALUES- 1 
2670 N1-CNT1-1: N2»CNT2-1: N3-CNT3-1 

2bS0 * 

2e«0 r CALCULATE STANDARD DEVIATIONS 

2700 SD1-<C1/Nl>~.3i SD2- (C2/N2) ~. 3i SD3- (C3/N3) ^. 3 

2710 ' 

2720 * COMPARE STANDARD DEVIATIONS TO SDLIMIT TO DETERMINE SLANT RANGE 
2730 * ACCEPTABILITY. 
2740 f 

2730 SDLIMIT-. 1 ' SET STANDARD DEVIATION LIMIT 
2760 ' 

2770 IF SD1>SDLIMIT THEN PRINT W DATA OUT OF RANGE - TRY AGAIN. " 5 GOTO 1620 
278* IF SD2>SDLIMIT THEN PRINTWDATA OUT OF RANGE - TRY AGAIN.": GOTO 1620 
27^0. IF SD3>SDLIMIT THEN PRINT"DATA OUT OF RANGE - TRY AGAIN.": GOTO 1620 
2300 ? 

2S10 ADD THE COUNTER DELAY EQUIVALENT OF 4.43 CM. TO EACH MEAN SLANT RANGE 
2S20 SMEAN1-SMEAN1+4.43: SMEAN2-SMEAN2+4 . 43: SMEAN3-SMEAN3+4. 43 
2330 ? ■ - 

2940 SP-3GAP+1 ? INCREMENT SRMEAN ARRAY POINTER 

2930 f 

23o0 * STORE MEAN VALUES IN ARRAY SRMEAN (MIKE*. SPARK* ) 
2370 " 

2330 SRMEAN ( 1 , SP ) =SMEAN 1 : SRMEAN ( 2 , SP ) -SMEAN2 : SRMEAN ( 3 , SP ) -SME AN3 
2890 ' 

2900 NEXT SGAP ' FIRE NEXT SPARK GAP 

2910 ' 
2920 RETURN 
2930 " 

2 Q 40 5 SUBROUTINE: DISTANCE 

2960 ' THIS SUBROUTINE DETERMINES THE DISTANCE BETWEEN TWO POINTS IN OBLIQUE 
2970 * SPARK GAP COORDINATES, 

2990 ' INPUTS: FOCAL (POINT, COORD) - OBLIQUE SPARK GAP COORDINATES OF FOCAL, 
3000 ? NORMAL AND ORIENTATION POINTS; 

ZOIO ' SLCPOINT) - SLANT RANGE DISTANCES FOR EACH POINT; 

3020 ' CW - COSINE OF ANGLE W, S2-S1-S3 (SPARK GAPS) ; 

3030 " S2W - SINE OF ANGLE W SQUARED 

3040 ' PFLAG - FLAG TO INDICATE REGISTRATION PROCEDURE 

5030 f OUTPUT: D (POINT) - DISTANCES FROM FOCAL, NORMAL AND ORIENTATION POINTS 

3060 7 

3070 ' POINT 1 

3080 ' CALCULATE PROJECTIONS OF SLANT 1 ON X AND Y SPARK AXES 
3090 » , . 

3100 XP1-(SL<1)^2+S31~2-SL<3)~2)/ <2tS31> ' EQN. h.* 
3110 YPI- (SL( 1) ~2+Sl2~2-SL (2) ~2> / (21S12) ' EQN. 6.4 
3120 7 

3130 ' CALCULATE OBLIQUE COORDINATES OF POINT 1 
3140 ' 

3130 X1-(XP1-YP1»CW)/S2W ' EQN. 6.8 

3160 Yl-(YPi-XPi«CW)/S2W * EON. 6.9 

3170 21- (SL ( 1 ) ~2-X 1~2-Y1~2-2*X I *Y 1 *CW) *** • 3 ' EQN. 6.12 
3180 ' 

3190 ' POINT 2 
3200 ' 

3210 ' FOR REGISTRATION PROCEDURE, ONLY THE FOCAL POINT IS NEEDED SO 
3220 9 PTS-l. THIS IS DETERMINED BY PFLAG. 

3230 ' 

3240 IF PFLAG- 1 THEN PTS-l ELSE PTS-3 
3^30 * 

3260 FOR 1-1 TO PTS p STEP THROUGH FOCAL, NORMAL AND ORIENT. POINT COORDS. 
3270 ' 

3280 ? RECALL FOCAL, NORMAL AND ORIENTATION POINT COORDINATES 
3290 ' 
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3300 X2-F0CALU. 1) 
3310 Y2-F0CAL(I ,2) 
3320 22=F0CAL<I,3> 
3330 ' 

3340 ' CALCULATE DISTANCE D USING OBLIQUE DISTANCE FORMULA , EQN. 6.14 
3330 > 

3360 D ( I ) « ( X2-X 1 ) ^2+ ( Y2-Y 1 ) -2* ( 22-2 1 ) ^2+2* iX2-Xl>«< Y2-Y 1 ) «CW 

3370 r 

3330 NEXT I 

3390 ? 

3400 RETURN 

3410 * 

3420 ' FILENAME: INVERSE 

3430 ' 

-440 ? THIS SUBROUTINE INVERTS THE JACOBIAN MATRIX USING THE CROUT ALGORITHM 
3430 ? 

3460 ' INPUT: J# (3,3) - JACOBIAN MATRIX 
3470 * OUTPUT: . I# (3. 3) - INVERTED MATRIX 
3430 ' 

3490 • SET INVERSION CRITEREA BY STORING THE IDENTITY MATRIX SOLUTIONS 
3300 FOR I-l TO 3 
3310 FOR J-i TO 3 

3320 IF I-J THEN I#<I.J>-1 ELSE I#(I,J>-0 
3330 NEXT J 
3340 NEXT I 
3330 ' 

3360 T0L#«lE-09 ' SET 2ER0 TOLERANCE FOR MATRIX SINGULARITY 

3370 ' 

3330 FOR I»i TO 3 * SEARCH FOR LARGEST ELEMENT IN A COLUMN 
3390 X— 1 
3600 * 

3610 FOR K-l TO 3 ' SEARCH BELOW MAIN DIAGONAL 
3620 IF ABS(J#CK, I))^»X THEN 3640 

3630 C-K: X-ABS ( j# <K, I ) ) 9 Q IS • ROW OF LARGEST ELEMENT 

3640 NEXT K ' X IS LARGEST ELEMENT 

3630 ' 

3660 TEST FOR SINGULAR MATRIX AND STOP IF NECESSARY 
3670 IF Xj-0 THEN 3690 ELSE PR I NT "SINGULAR MATRIX" : STOP 
3630 * 

3690 IF I-Q THEN 3790 * INTERCHANGE IF NEEDED 

3700 ' 

7710 FOR J-l TO 3 * NO. SWITCH. ROWS I AND Q 

3720 T-J#(I,J)i J# (X. J>-J#<Q, J) : J#<Q,J>-T 
3730 NEXT J 
3740 * 

3730 FOR J-l TO 3 

3760 T-I#<I.J)i I#(I, J)»I#(Q. J) : I#<Q.J)-T ' SWITCH RIGHT HAND SIDE 
3770 NEXT J 
3780 9 

3790 FOR J-l TO 3 ' ELIMINATE ON THAT ONE ROW 

3800 ■ 

3810 IF KJ THEN Ml-X-1 ELSE Ml-J-1 

3820 S«0 ' ' FIND INNER PRODUCT OF ROW I AND COLUMN J - S 

3830 ' 

3840 FOR K-l TO Mi 
3830 S«S*J# < I , K) * J# <K, J) 
3860 NEXT K 
3870 » 

3830 J# < I , J > «J# < I , J > +S 

3890 * STOP HERE IF BELOW MAIN DIAGONAL, CHECK FOR SINGULARITY ELSE NORMAL I 2E 
3900 IF I>-J THEN' 3930 
3910 * 

3920 IF ABS(J#<I,I))<TOL# THEN PR I NT "SINGULAR MATRIX": STOP 

3930 * 

3940 J#<I, J)— J#<I, J)/Jtt(X. I) 'ELSE NORMAL I 2E 
3930 NEXT J " NEXT COLUMN 
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3960 
3970 
3980 
3990 
4000 
4010 



NEXT I 



FOR J-l TO 3 
FOR I-l TO 3 



* REDUCE RIGHT HAND SIDE 



NEXT COLUHN 



FOR EACH SET OF CONSTANTS 
LOOK DOWN EACH COLUMN 



4020 S-0 

4030 ' 

4040 FOR K- I TO I-l 

4030 S«S+J#iI.K> *I#<K. J) ' TAKE PARTIAL INNER PRODUCT 

4O60 NEXT K 

4070 ? 

4080 I*(2.J)>-(Xtt(I. J)+S> /JO (1,1) 

4090 9 

4100 NEXT I 

4110 ? 

4120 FOR 1-3 TO 1 STEP -1 ? WORK BACK UP COLUMN 

4130 3»0 ' TAKE PARTIAL INNER PRODUCTS 

4140 ' 

4150 FOR K-I + l TO 3 

4160 . S«S+JtCX«K> tXftCK. J) 

4170 NEXT K 

4 i «o * 

4190 XftCX, J)— XftCX, J)«-S 



4270 ' 

4230 * THIS SUBROUTINE CALCULATES THE COORDINATES OF THE THREE MCROPHONES 
4290 ' IN THE CT COORDINATE SYSTEM AND STORES THEM IN FILE MIKE. 

4300 * 

4310 ' INPUT: FOCAL < 1 , COORD ) - OBLIQUE SPARK GAP COORDINATES OF FOCAL POINT 
4320 • OUTPUT: MIKE (MK, AX ) - CT COORDINATES OF MICROPHONES 
4330 ' SUBROUTINES NEEDED: SRINP, NEWTON 
4340 ' 

4330 PFLAG-1 ■ SET FOCAL POINT ONLY FLAG FOR SUBROUTINE DISTANCE 
4360 ' 

4370 * COLLECT FIDUCIAL DATA - - 

4330 9 

4390 FOR FID-1 TO 3 F STEP THROUGH EACH FIDUCIAL 

4400 PRINT 

4410 PRINT "FOCUS ON FIDUCIAL POINT #";FID 
4420 • 

4430 9 CALL SUBROUTINE SRINP TO COLLECT SLANT RANGE DISTANCES 
4440 ' 

4430 GOSUB 1320 9 CALL SRINP 
4460 ' 

4470 9 STORE AVERAGE SLANT RANGE DATA IN ARRAY SDATA 

4430 FOR MK-l TO 3 9 STEP THROUGH EACH MICROPHONE 

4490 FOR SPK-1 TO 3 * STEP THROUGH EACH SPARK GAP 

4300 SDATA (FID, MK, SPK) -SRMEAN (MK, SPK) * STORE SLANT RANGE DATA 

4310 NEXT SPK ' NEXT SPARK GAP 

4320 NEXT MK ' NEXT MICROPHONE 

4330 NEXT FID • NEXT FIDUCIAL 

4340 ' 

4330 PRINT 

4360 PR I NT "CALCULATING MICROPHONE COORDINATES..." 
4370 PRINT 

4580 ' CALCULATE DISTANCES BETWEEN MICROPHONES AND FIDUCIALS 
4390 * 

4600 FOR FID-1 TO 3 ' STEP THROUGH EACH FIDUCIAL 

4610 FOR MK-l TO 3 9 STEP THROUGH EACH MICROPHONE 




RETURN 



NEXT I 
NEXT J 



SUBROUTINE: REGISTER 



4640 
4630 

4q60 



enc -wc-i TO - * STEP THROUGH EACH SPARK SAP 

CALCULATE THE DISTANCES BETWEEN THE FIDUCIALS (FOCAL POINT) AND THE 
4670 > MICROPHONES. 
4650 ' 

4690 GOSUB 2940 * CALL DISTANCE 

SlO ' STORE DISTANCES TO FOCAL POINT IN ARRAY DIST 
4720 ' 

4730 DIST(FID,nK)-D(l) 

*ZVl J , F v T nw ' NEXT MICROPHONE 

473«J NEXT MK FIDUCIAL 
4760 NEXT FID . FIDUCIAL 

' IF TESTING THE REGISTRATION PROCEDURE, FIDUCIAL COORDINATES ARE 
4790 " KNOWN SO DO NOT ENTER COORDINATES, 
4800 * 

4810 IF CFLAG-1 THEN 4990 
48-0 * 

4330 5 ENTER FIDUCIAL COORDINATES 
4940 * 

ScO fVihI -ENTER CT COORDINATES OF THE FIDUCIALS" 

till FOP FID-1 TO 3 ' STEP THROUGH EACH FIDUCIAL 

4390 ? 

4<?00 PRINT 

4*10 FftINT"FIDUCIAL •"jFIDj-X, Y, 2 ; 

4«W ' STORE FIDUCIAL COORDINATES IN ARRAY FX 

£04l"l * 

4950 INPUT FX (FID. 1) ,FX (FID.2) .FX(FID.v) 

4?t5 NEXT FID ' NEXT FIDUCIAL 

: CALCULATE MICROPHONE COORDINATES BY SOLVING THREE NONLINEAR EQUATIONS 
' FOR THREE UNKNOWNS USING NEWTONS METHOD. EQN. 6.1-. 

5020 GOSUB 3910 ' CALL NEWTON 

5040 ' STORE MICROPHONE COORDINATES IN ARRAY MIKE 

Z '' Zt> . ,„ - ' STEP THROUGH EACH MICROPHONE 

5060 FOR MK-1 JO - , . STE p THROUGH X,Y,Z 

3S nS?^ 00 "" • S? nicrophone 

"i? • STORE CT COORDINATES OF MCRCFHONES IN FILE NIKE 

5140 OPEN "HIKE" FOR OUTPUT M .. | OPEN ^i^SSSoNE 

5150 FOR Wins , "«J tJbough X. Y AND 2 

S£ ^.St^.S.Sk.COORD, • S5« EACH VALUE TO THE F.LE 

5180 NEXT COORD 

US dS K 

3210 ? , R£SET F0CAL POINT COORDINATE FLAG 

3220 PFLAG-O- R&5c ' 
3230 ' 
3240 RETURN 
3230 * 

3260 * SUBROUTINE: SLICE 
3270 ' 



5230 " THIS SUBROUTINE CALCULATES THE CT COORDINATES OF THE FOCAL NORMAL 'ND 
2290 * ORIENTATION POINTS TO DEFINE THE FOCAL PLANE. IT ALSO DETERMINES 
3300 * THE FOCAL PLANE COEFFICIENTS *AND THE THREE DIRECTION COSINES. 
S3 10 

3320 ■" INPUTS: HIKE tMK, AX) - CT COORDINATES OF MICROPHONES; 

30 * FOCAL (POINT, COORD) - OBLIQUE SPARK GAP COORDINATES OF THE FOCAL 

**0 ' NORMAL AND ORIENTATION POINTS. ' 

3330 * OUTPUTS: FPOINT (COORD) - FOCAL POINT COORDINATES: 
53a0 ? FPLANE (COEFF) - FOCAL PLANE COEFFICIENTS; 

5370 ' DCOSX - X DIRECTION COSINE OF THE FOCAL - ORIENT. POINT VECTOR 

3380 DCOSY - Y DIRECTION COSINE OF THE FOCAL - ORIENT. POINT VECTOR 

2390 ' DCOSZ - 2 DIRECTION COSINE OF THE FOCAL - ORIENT. POINT VECTOR 

5400 * SUBROUTINES NEEDED: SRINP, NEWTON 
3410 * 

5420 " GET SLANT RANGES 
5 A 30 " 

3-*40 GOSUB 1520 * CALL SRINP 
5450 9 

34e0 PQR MK-l TO 3 " STEP THROUGH EACH MICROPHONE 

34T,;, FOR SPK-1 TO 3 ? STEP THROUGH EACH SPARK GAP 

3450 * STORE SLANT RANGE DATA IN VECTOR SL 
3490 SL ( SPK ) -SRMEAN < MK , SPK ) 

350v NEXT SPK ? NEXT SPARK GAP 

3310 * 

3320 GOSUB 2940 * CALL DISTANCE 
3330 ' 

3340 FOR PT-1 TO 3 • STEP THROUGH EACH POINT 

3330 ' STORE DISTANCES IN ARRAY 
5560 DIST<MK.PT)-D(PT) 

3370 NEXT PT * NEXT POINT 

3330 NEXT MK ' NEXT MICROPHONE 

33*0 * 

3600 GOSUB 3920 ? CALL NEWTON - EQN. 6.28 

3610 ' * 

3o20 ? STORE FOCAL POINT COORDINATES 

3630 FOR COORD* I TO 3 » STEP THROUGH X.Y,2 

3640 FPO I NT ( COORD ) -NSOL ( 1 . COORD ) 

3*30 NEXT COORD * NEXT COORDINATE 

5660 ' 

3670 FOR I«t TO 3 ? CALCULATE AND STORE FOCAL PLANE COEFFICIENTS 

3680 FPLANE ( I ) -NSOL (2.1) -FPOINT ( I ) 
56*0 NEXT I 
3700 ? 

3710 FPLANE (4) -0 

3720 FOR 1-1 TO 3 • STEP THROUGH EACH COEFFICIENT 

3730 FPLANE (4) -FPLANE ( 4 ) +FPLANE ( I ) I FPO I NT ( I ) 
3740 NEXT I > NEXT COEFFICIENT 

3730 9 

3760 ' CALCULATE AND STORE X, Y AND 2 DIRECTION COSINES 
3770 9 

3730 NORM-0 * INITIALIZE VECTOR NORM 

3790 • 

3SO0 FOR COORD- 1 TO 3 " STEP THROUGH X,Y,Z 

5310 NORM-NORM* ( NSOL ( 3 • COORD ) -NSOL ( X • COORD ) ) ~2 * EQN . 6 . 29 
3820 NEXT COORD * NEXT COORDINATE 

3330 ' 

3840 NORM-NORM-". 3 
3330 9 

3360 DCOSX- (NSOL (3.1) -NSOL (1,1)) /NORM ' EQN. 6.30 
3870 DCOSY- (NSOL (3. 2) -NSOL (1,2) ) /NORM • EQN. 6.31 
3330 DCOSZ- (NSOL (3,3) -NSOL (1,3)) /NORM * EQN. 6.32 
3890 • 
3900 RETURN 
5910 * 

5920 * SUBROUTINE: NEWTON 
3930 ■ • 
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5340 f THIS SUBROUTINE SOLVES 3 SETS OF 3 NONLINEAR, EQUATIONS FOR 3 UNKNOWNS. 
3930 ' USING AN ITERATIVE NEWTON'S METHOD. THE EQUATIONS ARE IN THE FORM OF 
3*60 ' DIST«((FXl-NS0Ll)-2*<FX2-NS0L2)^2*(FX3-NS0L3)^)-".3, WHERE DIST AND FX 
3970 ' ARE KNOWN AND NSOL IS UNKNOWN. 

^?30 ' 

5*90 * INPUTS: FX (POINT, COORD) - FIXED COORDINATES; 

6000 9 DIST (POINT. 3) - DISTANCES BETWEEN POINTS; 

6010 * CW - COSINE OF. ANGLE W, S2-S1-S3 (SPARK GAPS) 

6020 ' OUTPUT: NSOL (POINT, COORD) - COORDINATE SOLUTIONS 

©030 ? SUBROUTINE NEEDED I INVERSE 

6040 " 

oOSO ? SET CONSTANTS 

o060 E-.OOl * ITERATIVE TOLERANCE 
6070 ' LOOP THROUGH EACH SET OF EQUATIONS 
oU80 FOR MK»1 TO 3 

60«0 * * ■ 

a 100 ' STORE INITIAL SOLUTIONS 

61 IV FOR PT-l TO 3 

O1--0 M#U.PT)-13 

6130 NEXT PT 

6140 * 

6130 ? LOOF THROUGH UP TO 33 ITERATIONS 
a 160 FOR ITEP»2 TO 33 - 
el70 ' 

zlSO ' STORE ITER- I SOLUTION 
6l<?0 FOR PT»1 TO 3 
i2u0 MV# ( PT * -n* ( I TER- 1 , PT ) 
a2l0 NEXT PT 

o230 > CALCULATE JACOBIAN MATRIX COEFFICIENTS FROM THE DERIVATIVES OF THE 
6240 ' EQUATIONS. 
6250 FOR I»l TO 3 
6260 FOR J-l TO 3 

6270 J#U, J)*2» (MV#(J)-FX(I, J) ) 

6230 NEXT J 
6290 NEXT I 
6300 ' 

6310 ' INVERT THE JACOBIAN MATRIX 
6320 GOSUB 3420 ' CALL INVERSE 

6330 ' 

o340 * DETERMINE SOLUTION ITER 

till "^1^2!^ 
• 

o370 NEXT I 
6380 ' 

"JlITER^-MVta^ 
6410 NEXT I ' 

£sO ' CALCULATE ITERATIVE ERROR AND SEE IF IT IS WITHIN THE SET TOLERANCE TOL 
6440 FOR I-i TO 3 wm v \ \ 

6430 T ( I ) -ABS (M# (ITER, I > -MV# ( I ) ) /ABS <M# < ITER. I ) > 
O 460 NEXT I 

6480 'IF ERROR IS WITHIN TOL STORE SOLUTION OR CONTINUE TO NEXT ITERATION 
6490 IF (T(IXE) AND (T(2XE) AND (T(3)<E) THEN 6380 
6300 * 

6310 NEXT ITER ' NEXT ITERATION 

6330 ' IF NO SOLUTION IS REACHED AFTER 33 ITERATIONS PRINT ERROR MESSAGE 
6340 * AND STORE LAST ITERATION. 

o330 PRINT "NO SOLUTION FOR POINT #";MK: ITER-33 
6360 9 

6370 ' STORE SOLUTIONS IN ARRAY NSOL 
638» FOR 1-1 TO 3 



167 



e s*u nsol Km, I ) a n# i iter, i ) 

6e00 next i 

6620 NEXT HK " NEXT SET OF EQUATIONS 

ooZO * 

6640 RETURN 

ceSO 7 

oabO ' SUBROUTINE: TEST 



■SS THIS UMUTM I. — ™ "J SS^SSSfUTS! POINTS ARE ^NOUN 

I.To • S5£ .^^.^Sf^ES^ « SLn. error „ « 



o720 
a 730 

o7e ' 



POINT. 



SUBROUTINES NEEDED: REGISTER. SLICE 



°*-rt if HFLAG-1 GOTO 6920 

SJS:-— « PHANTON S.OCK „™ RESPECT TO THE n.-H- , 

6800 PRINT 



• STORE "CUCMU COORDINATES ^ 



oJ 1 . 1 Stunt r- « - — «- % — f> 

l£S SET CALIBRATION FLAG FOR SUBROUTINE DISTANCE 
~3S0 CFLAG=»1 

SCTUS .-to • REGISTER-DETERMINE H.CROP-ONE COORDINATES 

i? 30 FOR 4-1 TO Z 
6 940 FX'.I,J)-"lKE(I ? Oi 
6=30 NEXT J 

5 c©0 NEXT I 

' LOOP THROUGH UP TO 100 TEST POINTS 
FOR TPOINT-1 TO I 00 

'»WT "ENTER TEST POINT . «.-.. « ««"«• " 

7^20 INPUT TPT* . - 

' LOOK- UP KNOWN COORDINATES OF TKT POINT ^ 7iS0 

r03^ IF TPT.-" I" THEN TX--4 ; 719: TY-7.62. Ti ^ 
7?60 F TPT.-2- THEN TX-2.S4; TY-O: TZ-O^ 6 gqt0 ?l50 

-070 IF TPT.--3" THEN TX-5.949 *~ , ^Z-O. GOTO 7130 

7X30 GOTO 7010 

?5S ' G OSUB =260 ' CALL SLICE TO CALCULATE VALUES 
SI?? » PRINT FOCAL PLANE COEFFICIENTS 
?JK PrSnT-FOCAL PLANE COEFFICIENTS <A,B,C,D>" 
g?S PRINT FPLANE ( 1 > » FPLANE ( 2 ) » FPLANE ( 3 ) ? FPLANE ( 4 ) 
7220 ' 

7230 PRINT _ , v „ 

•7240 PRINT- -DIRECTION COSINES <X,Y,Z> 
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7230 
72oO 
7270 
7230 
7290 
"300 
7310 
* .'2v 
"330 
"34*;. 
7330 
."360 
"370 1 

"30 NE»T TPOINT 
~3?0 * 

"iOO RETURN 



RRINT 

RRINT DCOSX.DCOSY.DCOSI 
PRINT 

' RRINT TRUE AND CALCULATED VALUES OF THE FOCAL POINT AND THE ERROR IN nn. 
PR I NT" FOCAL POINT COORDINATES" 
RRINT 

PRINT " MEASURED- " I FPO I NT \ 1 ) ,FP0INT(2) ,FP0INT<3) 
PRINT " TRUE- M ? TX , TY , TZ 

EF»lO* ( iFPOINT(l)-TX? ^2+ iFPOINT <2) -TY) ~2+ (FPOINT (3) -T2> ~2V% 5 
PRINT " ERROR IN mm.-":ER 
PRINT 
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•430 
"-40 
'450 
~4a0 



SUBROUTINE: RECALL-SIT 

THIS SUBROUTINE RECALLS THE COORDINATES OF THE FOCAL, NQRriAL AND 
ORIENTATION POINTS FOR THE "SITTING" SPARK GAP BRACKET FROM THE 
FILE SIT AND STORES THE COORDINATES IN ARRAY FOCAL (POINT. COORDINATE" 1 



"43-V OPEN "3IT N FOR INPUT AS #1 
"4-0 FOR PNT** TO 3 



120 



FOR COORD- 1 TO 3 

INPUT # I , FOCAL ( PNT . COORD ) 
NEXT COORD 
NEXT PNT 
CLOSE *1 
RETURN 



OPEN FILE -SIT" TO READ IN COORDINATES 
LOOP THROUGH FOCAL. NORMAL AND ORIENT FT. 
LOOP THROUGH X. Y AND 2 COORDINATES 
STORE DATA IN ARRAY FOCAL (POINT, COORD ■ 
NEXT COORDINATE 
NEXT POINT 
CLOSE FILE 



• SUBROUTINE: RECALL-SUP 



7380 
7390 
"aOO 
7ol0 
7620 
7o30 
"640 
"530 
"660 
"670 
7680 
'690 
7700 
7710 
772-J 
'730 
7740 
^730 
7760 
^770 
"790 
7790 
7300 
7810 
7820 
7330 
7B40 
7330 



THIS SUBROUTINE RECALLS THE COORDINATES OF THE FOCAL. NORMAL AND 

ORIENTATION POINTS FOR THE "SUPINE " SPARK GAP BRACKET FROM THE 

FILE SUP AND STORES THE COORDINATES IN ARRAY FOCAL (POINT, COORDINATE) ■ 



OPEN "SUP" FOR INPUT AS #1 
FOR PNT-1 TO 3 
FOR COORD- 1 TO 3 

INPUT * 1 « FOCAL (PNT , COORD ) 
NEXT COORD 
NEXT PNT 
CLOSE *1 
RETURN 

' SUBROUTINE: RECALL-MIKE 

f THIS SUBROUTINE RECALLS THE CT COORDINATES OF THE MICROPHONES THAT WERE 
' STORED IN CASE THE COMPUTER FAILED AFTER THE REGISTRATION PROCEDURE. 



? OPEN FILE "SUP" TO READ IN COORDINATES 

* LOOP THROUGH FOCAL, NORMAL AND ORIENT PT. 
? LOOP THROUGH X, Y AND 2 COORDINATES 

* STORE DATA IN ARRAY FOCAL (POINT, COORD) 
' NEXT COORDINATE 

* NEXT POINT 
" CLOSE FILE 



OPEN "MIKE" FOR INPUT AS #1 
FOR MK-1 TO 3 
FOR COORD- 1 TO 3 

INPUT #1. MIKE (MK, COORD) 
NEXT COORD 
NEXT MK 
CLOSE *1 
HFLAG-I 
RETURN 



OPEN FILE "MIKE" TO READ IN COORDINATES 
LOOP THROUGH THE MICROPHONES • 
LOOP THROUGH X, Y AND 2 COORDINATES 
STORE DATA IN ARRAY MIKE (MK, COORD) 



CLOSE FILE 

SET RECALL FLAG EQUAL TO 



rs.13 appendix contains the mechanical drawing for the 
spar. 4 ap holder, spar, gap holder brackets, spar* sap 
.older - alcrosccpe aount, and the celling microphone 
mounts. 
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Appendix E - Jixitizer Error Statistic s 

Listed below are the results of the digitizer 
evaluation. It was determined that the digitizer meets its 
specifications for precision and that the observed errors in 
accuracy are most likely due to measurement errors in 
determining the correct slant range distance, and not due to 
digitizer error. All values are in centimeters where d is 
the slant range distance, n is the sample number and sd is 
the standard deviation. 
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Value 


Value 




Theoretical 
Experimental 


50 


50 
50 
.130 
150 


50 .000 
49 .339 
49 .591 
49 .900 


0.002T 
0.0059 
0.0028 


50.01 
49.39 
49.90 
49.91 


49.99 
49.33 
49.38 
49.39 


0.02 
0.01 
0.02 
0.02 


Theoretical 
Experimental 


100 


50 
50 
100 
150 


100.00 
99 .904 
99 .915 
99 .929 


0.0069 
0.0053 
0.0071 


100.01 
99.91 
99.93 
99.92 


99.99 
99.39 
99.90 
99.90 


0 .02 
0.02 
0.02 
0.02 


Theoretical 
Experimental 


150 


50 
50 
100 
150 


150.00 
149.63 
149.55 
149.63 


0.0038 
-0.0067 
0.0048 


150.01 
149.64 
149 .66 
149.69 


149.99 
149.62 
149.62 
149.62 


0.02 
0.02 
0.02 
0.02 


Theoretical 
Experimental 


200 


50 
50 
100 
150 


200.00 
199.30 
199.32 
199.33 


0.0061 
0.0045 
0.0063 


200.01 
199.30 
199.33 
199.34 


199.99 
199.23 
199.31 
199.32 


0.02 
0.02 
0.02 
0.02 
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Ap pendix f - f xpe r imencal Jata 

This Appendix is divided into four sections, focal 
point accuracy and precision, resting focal point precision 
focal plane accuracy and spark gap holder precision. The 
experimental criteria is explained with each section. 

focal Point Accuracy 

This purpose of this experiment was to -determine the 
reference system accuracy and precision at the focal point, 
first the test phantom coordinates of the microphones were 
determined by focusing on each of the three test fiducials. 
Then the microscope was focused at ten different angles on 
four test points whose coordinates are known in. the test 
coordinate system. The table below shows the magnitude of 
the focal point errors in millimeters for each sample, n, 
and the mean and standard deviations Csd) for each test 
point. 

Test Point 



1 -fS 2.47 .0-99 

0.53 1.22 1.23 

1 .33 0.9V 0.74 

1.93 0.57 l.^J 

0.66 0.8b 1.0o 

0.45 • 1.29 ■ 1.59 

0.48 0.38 1.60 

0.53 0.48 0.72 

1.43 1.73 2.23 

0.86 0.91 2.W 



1 


0.79 


2 


0.53 


3 


1 .45 


4 


0.13 


5 


0.67 


6 


0.49 


7 


1 .00 


8 


0.57 


9 


1.10 


10 


0.33 
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Mean 0.77 1.01 1.0* 1.40 

~"sd~ 0.556 0.566 0.6j4 0.545 

Average Mean 1.07 

Average sd 0.526 

Resting Focal Point Precision 

The purpose of this experiment was to determine the 
reference system precision at the focal point, in other 
words, how much do che coordinates of the focal point vary 
without refocusing the microscope between samples. After 
registering the microphones, the microscope was focused on 
one test point and fifteen samples taken. The standard 
deviations were then calculated to give an estimate of the 
variations that could be expected during a typical case due 
to variable air motions, etc. Each spark gap was fired _$0 
times. The number of samples is n, the coordinates are in 
centimeters, and the standard deviation is sd. 



Coordinates 



_n 


X 


x 


Z 


1 


-4.65 


7 . 76 


0.73 


2 


-4. 70 


Y .64 


0.69 


3 


-4 . 67 


7.76 


0. 72 


4 


-4.63 


7.75 


0.72 


5 


-4.64 


7.69 


0.63 


6 


-4.53 


7.77 


0.61 


7 


-4.60 


7.30 


0.69 


8 


-4.65 


7.70 


0.72 


9 


-4.66 


7.75 ' 


0.73 


10 


-4.62 


7.74 


0.69 


11 


-4 . 67 


7.78 


0.76 


12 


-4.63 


7.71 


0.72 


1j 


-4.65 


7.63 


0.58 


m 


■ -4.69 


7.74 


0.75 


15 


-4.56 


7.65 


0.62 
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sd J. 057 O.Ooo 3.055 

focal Plane Accuracy 

The purpose of this experiment was to determine the 
accuracy of zhe focal plane calculations. Since the focal 
plane is defined. by four coefficients in the form Ax + 3y + 
Cz = D, it is difficult to determine the plane accuracy. 
Therefore, the microscope was focused on the grid coordinate 
system (figure 9.1) normal to the z axis where the equation 
of the plane is z = 0, and the plane coefficients determined. 
Ideally, A, 3 and D should equal 0 and C equal 1. 

The accuracy of the focal plane orientation was 
determined by calculating the x and y axis direction cosines 
of the focal point - orientation point vector in the first 
quadrant of the grid coordinate system. The sum of these 
angles should be theoretically 90°. 

The focal plane coefficients and sum of direction 
cosine angles in degrees are indicated below for fifteen 
refocused samples (n), along with the means and standard 
deviations (sd) . 

focal Plane Coefficients 



HA BCD Angle 

1 -.044 .101 1 -.042 34.j 

2 -.037 .108 1 .053 90.3 

3 -.025 .090 1 -.060 34.5 

4 .113 1 -.015 91.4 

5 -.030 .121 1 -.053 91.0 
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6 
7 
3 
9 

10 
11 
12 
13 
14 
15 



-.035 
-.031 
-.029 
-.040 
-.041 
-.041 
-.031 
-.027 
-.022 
-.021 



.091 
.113 
.109 
.100 
.900 
.102 
.101 
.105 
.118 
.119 



-.015 
.064 
.021 
.037 
.01 7 

-.112 
.01 7 
.019 
.035 
.018 



91 .2 
91 .2 
92.3 
90.7 
91 .9 
91.0 
90.7 
95.7 
.92.6 
92.6 



Mean 
sd 



-.0321 
.0072 



.1061 
.01 10 



.0053 
.0503 



90.33 
2.39 



Spark Gap Holder Precision 

The purpose of this experiment was to quantitatively 
determine the error due to repositioning the spark gap 
holder between oblique focal, normal and orientation point 
determination and the actual procedure. The focal point 
coordinates (in cm.) were determined fifteen times without 
refocusing the microscope on the spark gap 4 and the average 
values were -3.509 (x), -3.723 (y) and -39.371 (z) with 
standard deviations of 0.035 (x), 0.035 (y) and 0.022 (z). 
The spark gap holder was then removed and the sterile drape 
bag placed over the microscope. The spark gap holder was 
replaced after the drape bag was stretched over the metal 
tongue of the microscope mount. The focal point coordinates 
(in cm.) were determined fifteen times after the microscope 
was refocused on spark gap 4 and the average values were 
-3.430 (x), -3.733 (y) and -39.391 (z) with standard 
deviations of 0.053 (x), 0.053 (y) and 0.025 (z). fhe 



1ci2 

magnitude of error is the square root of the sum of the 
squares of the differences in x, y and z. 

'S-holder 1 =Vo.006j 2 + 0.0030 2 + 0.0004 2 
'^s-holder 1 = °' 9 »iHi«etera. 

This error might be reduced by designing .a better spark 
gap holder - microscope mount once the limitations and 
physical constraints of the operating room have been 
experimentally determined . 



Appendix G - Digitizer Operator's Manual 

This appendix is a step by step guide to the setup and 
operating procedure for the reference-display system. The 
relationship between the spark gaps and the focal, normal 
and orientation points must be determined first by running 
file FPL and following the instructions ' given . Below is a 
description of the procedure in the operating room. 

Display System 

1. Place the microscope in an upright position, 

2. Remove the ocular assembly and camera-observer 
beam splitting assembly. 

3. Attach the display assembly, with the projection . 
tube to the kft (See Figure 3.5), to the microscope. 

4. Replace the camera-observer beam splitting assembly 
and the ocular assembly. 

5. Connect the +12 volt d-c power supply and the 
composite video input cable (from the Interact TV patch 
panel on the operating room ceiling) to the display 
assembly. 

Spark Gaps 

6. Attach the microscope-spark gap link (unsterile) to 
the laser mount on the microscope, securing the set screw. 

I. Place a sterile drape bag over the microscope. 
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3. Pull che sterile drape over the "tongue" of the 
microscope-spark gap link until a hole is made in the drape 
oag. Uote: be careful not to touch the unsterile 
"tongue" , cutting the hole partially open will help.) 

9. Attach the sterile spark gap holder to the 
mlcroscope.5p.rk gap link and tighten the set screws into 
the recessions on the "tongue". 

10. Attach the spark gaps to the spark gap holder, 

tightening each set screw. 

11. Balance the microscope. 

Digitizer - Compu ter System 

12. Attach the four microphones to the microphone 
mounts on the operating room ceiling. 

lj. Connect the ribbon microphone cable from the 
aerophones to the DB-25 connector on the front panel of the 
digitizer. 

14. Connect the three" spark gaps to the front panel of 
the spark gap multiplexer. 

15. Connect the hex connector spark gap input fro- the 
sack panel of the multiplexer to the front panel of the 

digitizer . 

16. Connect the digitizer output (back panel- terminal 
port) to the RS-2J2 serial port on the back of the I3M PC 
XT, via the DU-25 connector ribbon cable. 



W. Attach the parallel port of the IBM to the back 
panel of the multiplexer. 

18. Set the multiplexer to computer mode. 

19. Turn on the IBM with a system disk in drive A. The 
screen will ask for you to enter the date and time which can 
be ignored by pressing "RETdRN". 

20. On the back panel of the digitizer, make sure the 
foot pedal plug (without cable) is in place and the 
"terminal" switch is on. On the front panel, make sure that 
the toggle switch is on "line" and the rate knob turned 
completely clockwise (maximum). 

21. Type "C:" in response to the prompt "A>" and then 
type "FPL" for the focal plane determination, or "REG" for 
the surgical procedure. 

22. Follow the directions of the software for the rest 
of the procedure. 

Note: It is best to test the registration system (File REG) 
after determining the spark gap - focal plane relationship 
(rile FPL) to see if the error at the test points is 
approximately one millimeter. 
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